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The  energy  basis  for  a  coastal  region  at  Franklin 
County  and  Apalachicola  Bay,  Florida  was  evaluated  with 
models,  field  measurements,  energy  calculations,  land-use 
maps  from  aerial  photographs,  and  computer  simulations. 
Energetic  relationships  between  urban  and  coastal  ecosys- 
tems were  studied  to  help  understand  ecosystem  relationships 
on  a  regional  scale  and  predict  consequences  of  adding  ur- 
ban development  and  changing  water  quality  to  an  economy 
based  on  a  marine  fishery.   Simulations  of  a  regional  model 
suggested  a  near  steady  state  pattern  of  urban  activities 
in  the  county  which  was  consistent  with  previous  studies. 
A  slow  decline  in  urban  activities  was  recorded  when  fur- 
ther inflation  was  introduced.   Trends  were  level  tempor- 
arily when  investments  and  government  spending  were  added. 
Temporary  growth  was  observed  when  the  investment  rate  was 
increased  to  twice  the  present  rate.   Growth  could  not  be 
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stimulated  by  increasing  oyster  harvest  effort  without  devel- 
oping special  feedback  to  stimulate  the  oyster  reefs.   Growth 
trends  initiated  by  investment  in  urban  development  were  re- 
versed and  a  decline  observed  if  the  oyster  industry  was 
closed  because  of  a  decrease  of  water  quality.   Reducing  the 
seasonal  pulse  and  increasing  and  decreasing  river  flow  also 
reduced  oyster  biomass  and  the  urban  economy.   The  model  in- 
dicated that  natural  predation  rather  than  harvest  was  the 
major  cause  of  mortality  of  oysters. 

An  aggregated  submodel  of  diurnal  dissolved  oxygen 
dynamics  in  Apalachicola  Bay  was  programmed  to  investigate 
the  possibility  of  low  oxygen  occurring  under  a  variety  of 
conditions.  Community  metabolism  and  organic  carbon  stocks 
in  the  bay  and  river  were  measured  in  the  summer  of  1973  to 
evaluate  models.  Metabolism  ranged  from  3.1  to  21.6  g  O2/ 
m^/day  and  was  generally  comparable  to  other  estuarine  areas. 

Low  predawn   oxygen  levels  were  not  observed.   Total  organic 

3 
carbon  averaged  6.4  and  5.3  g  C/m  in  the  bay  and  river, 

respecti  vely . 

Simulation  of  the  model  with  reduced  diffusion  and  in- 
creased detritus  input  depressed  oxygen  levels  to  3.2  ppm. 
Nutrient  loading  produced  larger  diurnal  changes  but  not 
lower  oxygen  levels.   Evaluation  of  this  model  suggested 
that  phytoplankton  photosynthesis  was  the  major  source  of 
organic  matter  in  the  summer  (75%).     Eighty  percent 
of  the  nitrogen  used  in  photosynthesis  was  from  the  recycle 
of  nitrogen  from  imported  organic  matter. 
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Energy  calculations  indicated  that  natural  energy  in- 
puts accounted  for  84%  of  the  total  energy  input  in  the  coun- 
try.  The  ratio  of  purchased  to  natural,  free  energy  inputs 
(investment  ratio)  for  the  county  was  low  compared  to  Flori- 
da and  the  nation.   On  a  county  basis,  full  development  of 
the  barrier  island  for  tourism  changed  the  ratio  from  0.19 
to  3.1,  exceeding  the  2.5  national  average.   On  a  local  basis 
calculations  indicated  that  the  planned  development  may  be 
too  dense  for  best  fit  of  natural  and  purchased  resources. 

The  oyster  industry  accounted  for  60%  of  the  total 
county  income  in  1970.   Within  the  county,  natural  energies 
predominated  and  accounted  for  84"i  of  the  total  energy  used 
in  oystering  as  measured  at  dockside.   When  traced  to  final 
consumption  the  natural  energies  involved  in  producing  oy- 
sters were  matched  elsewhere  in  the  national  economy  with 
the  work  of  fossil  fuels  in  ratio  of  1:2.2.   Preliminary 
calculations  suggested  that  the  bay  could  support  twice  the 
current  oyster  population  without  adding  extensive  fossil 
fuel  subsidies. 
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INTRODUCTION 

This  is  a  study  of  principles  organizing  man  and 
nature  in  the  coastal  zone.   The  study  considers  purchased 
and  natural  energy  flows  that  together  support  a  coastal 
fishing  economy.   Questions  are  considered  for  fitting  man's 
increasing  interaction  with  coastal  waters  into  patterns 
which  maximize  regional  value.   Included  here  are  land- 
use  maps  developed  from  aerial  photographs,  an  aggregated 
model  of  urban-coastal  relationships,  energy  analyses  of 
development  alternatives,  simulation  models  of  subsystems, 
and  examination  of  ways  that  changes  in  international  energy 
regimes  control  county  activity.   This  study  was  made  to 
develop,  apply,  and  extend  ecological  principles  cf  organ- 
ization and  function  to  a  regional  scale. 

Over  the  world,  during  the  period  of  accelerating 
urban  growth  based  on  fossil  fuels,  there  was  an  increasing 
harvest  drawn  from  both  recreational  and  commercial  marine 
fisheries  at  the  same  time  that  coastal  developments  were 
removing  habitat  and  interfering  with  the  necessary  food 
chain  base  of  coastal  fisheries.   In  many  areas  there  was 
a  burst  of  economic  activity  as  boats  utilized  a  virgin 
fishery  followed  by  periods  of  overfishing,  coastal 
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disturbance,  loss  of  the  fishery,  and  replacement  with  ur- 
ban development  that  may  have  unnecessarily  decimated  the 
fishery  resource.   We  need  to  test,  clarify,  and  recommend 
ways  to  retain  the  value  of  marine  resources,  while  adding 
fossil  fuel-based  developments. 

However,  in  1974  with  changing  patterns  of  fuel 
availability  and  prices,  the  relative  values  of  the  natural 
work  and  fossil  fuel-based  work  began  changing,  and  we  need 
an  understanding  of  the  trends  ahead.   As  trends  reverse, 
strong  marine  resources  may  be  needed  as  an  economic  cushion 

This  study  considered  the  general  problem  of  co- 
existence of  economic  development  and  marine  resources 
using  Franklin  County,  Florida  which  includes  the  mouth  of 
the  Apalachicola  River  and  a  very   valuable  resource  of  oy- 
ster and  shrimp  landings.   Major  external  developments 
affecting  Franklin  County  include  changing  prices  of  pur- 
chased goods  and  fuels,  river  flow  and  coastal  modifications 
development  of  coastal  islands  for  tourism,  and  changing 
populations  pressures  for  tourist  and  retirement  develop- 
ments.  To  gain  understanding  of  the  energetic  base  of  this 
area,  system  models  were  constructed  as  follows:   (1)  an 
oxygen  balance  model  of  production,  respiration,  exchange 
and  diffusion  in  Apalachicola  Bay,  with  special  emphasis 
on  the  role  of  the  river  as  a  source  of  nutrients  and  or- 
ganic matter;  and  (2)  a  regional  model  showing  trends 
caused  by  external  driving  functions  on  human  activity  on 


land  interacting  with  estuarine  and  fishery  acti vi ties  that 
together  constitute  the  economy  of  Franklin  County,  Florida 
Using  numerically  evaluated  models,  energy  analyses  were 
made  of  the  primitive  and  present  patterns  of  land-use  in 
the  county,  the  oyster  industry,  and  a  proposed  tourist 
development  on  the  barrier  island.   Land-use  maps  were  con- 
structed from  aerial  photographs  to  show  spatial  relation- 
ships of  the  primitive  condition  and  current  use  by  systems 
of  nature  and  man.   The  energy  quality  of  the  potential 
energy  of  river  head  was  also  calculated  and  used  in  the 
energy  analyses.   Lastly,  calculations  were  made  using  an 
energetic  basis  to  estimate  the  appeal  of  the  coastal  zone 
as  a  place  for  fossil  fuel-based  development. 

Theoretical  Coastal  Zone  Issues 


In  evaluating  the  energy  basis  of  coastal  Franklin 
County,  three  theoretical  issues  were  studied  as  follows: 
(1)  A  seasonally  pulsing  energy  flow  may  accelerate  the 
net  yield  of  a  system.   The  pulsing  discharge  of  the 
Apalachicola  River  may  be  an  example  affecting  net 
yields  of  the  bay.   What  is  the  energy  value  of  a 
pulsing  regime  with  net  yield?  How  does  it  compare 
to  other  river  programs  which  have  regular  flow  and 
less  net  yield?  Does  the  system  receiving  the  pulse 
benefit  from  it?  What  are  the  feedbacks  or  character- 
istics exported  to  insure  continued  pulsing  programs? 


Is  there  a  direct  or  indirect  rate  of  the  yield  re- 
lated back  to  the  pulsing  inflow? 

(2)  A  second  question  is  the  relationship  of  the  image  of 
the  coastal  zone  to  energy  flow.   Is  it  possible  that 
this  attractiveness  or  image  is  based  on  the  richness 
both  in  types  and  magnitudes  of  natural  coastal  zone 
energy  flows  which  find  viable  couplings  with  urban 
devel opment? 

(3)  The  third  question  concerns  the  forward  and  feedback 
exchange  of  value  between  estuarine  productivity  and 
its  coupled  fishery  harvest.   In  agriculture  there 

is  a  well-known  give  and  take  between  farmer  and  field 
with  yields  related  to  feedback  studies  of  soil  prepara- 
tion, fertilization,  breeding,  weed  control  and  the 
like.   What  are  the  feedbacks  from  urban  areas  re- 
quired to  maintain  yielding  estuarine  processes?  Can 
additional  feedbacks  be  established  to  promote  net 
yields  and  prevent  deterioration  of  fisheries? 

Description  of  Study  Area 


Franklin  County  is  located  on  the  northwest  Gulf 
Coast  of  Florida  (Fig.  1).   It  is  bounded  on  the  north  by 
Liberty  and  Wakulla  Counties,  on  the  west  by  Gulf  County 
and  on  the  east  and  south  by  the  Gulf  of  Mexico.   The  county 
is  rural  with  most  of  the  7,000  residents  living  along  the 
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coast.   The  gross  area  of  the  county  consists  of  361,000 
acres,  of  which  18,400  acres  are  freshwater  rivers  and 
lakes.   The  county  consists  mostly  of  flatlands,  with 
swampy  forests  covering  over  90%  of  the  land  area.   Apa- 
lachicola  Bay,  a  medium  salinity  estuary  separated  from 
the  Gulf  of  Mexico  by  several  offshore  islands,  extends 
nearly  the  entire  length  of  the  county. 

The  bay  is  approximately  36  miles  long,  1  to  14 
miles  wide  and  covers  an  area  of  104,320  acres  (163  square 
miles)  including  all  of  St,  Vincent  Sound  on  the  west.  East 
Bay  on  the  north,  Apalachicola  Bay  proper,  and  St.  George 
Sound  as  far  east  as  the  western  tip  of  Dog  Island.   Of 
the  104,320  acres,  9,400  acres  contain  submerged  vegetation 
most  of  which  occurs  in  St.  George  Sound.   Species  include 
Thalassia  testidinum,  Syringodi  urn  f i 1 i forme.  Pi  pi anthera 
wrighti  i ,  and  Halophi la  engelmanni  i .   Major  emergent  spe- 
cies cover  approximately  21,000  acres  of  marsh  and  include 
Spartina  al  ternif  lora  .  S^.  patens  ,  Juncus  roemarinus  , 
D i s t i  c h  1  i s  spicata  ,  and  S a  1  i c o r n i a  perennias  (McNulty  et_ 
al.  ,  1972). 

Gorsline  (1963)  estimated  that  river  water  remains 
in  the  bay  system  for  a  few  days  in  the  winter  and  up  to 
a  month  in  the  summer.   However,  the  oyster  producing  area 
of  Apalachicola  Bay  is  smaller  than  the  bay  area  defined 
by  Gorsline  and  the  retention  time  for  this  area  may  be 
less.   The  mean  depth  of  the  bay  system  at  low  water  is 


2.3  meters  (7.5  feet).   During  the  winter  season  the  bay 
appears  to  be  well  mixed  but  has  salinity  induced  strati- 
fication around  the  inlets.   During  the  warmer  seasons, 
stratification  is  generally  observed  with  salt  wedge  pene- 
tration extending  into  East  Bay  at  times  (Estabrook,  1973). 
This  condition  may  be  a  recent  feature  as  both  previous 
hydrologic  studies  (Dawson,  1955a;  Gorsline,  1963)  reported 
no  stratification  except  in  dredged  channels.   High  turbidi- 
ty was  associated  with  river  discharge,  being  highest  in 
the  winter  (Estabrook,  1973).   Surface  circulation  in  the 
bay  is  partly  wind-driven  and  partly  dri  ven  by  the  river  dis- 
charge flowing  to  the  westward  most  often.   Water  flows 
east  when  west  winds  are  strong.   Bottom  salinities  had  a 
different  distribution  from  the  surface  except  when  the 
bay  was  well  mixed  (Estabrook,  1973). 

The  Apalachicola  River  has  an  average  yearly  dis- 
charge of  26,713  cfs  (cubic  feet  per  second)(6.53  x  10 

q 

m  /day)  draining  a  watershed  of  18,000  square  miles  in- 
cluding parts  of  Georgia  and  Alabama  (Hawkins,  1973)(Fig. 
2).   River  flow  had  a  well-defined  seasonal  cycle  with 
highest  and  lowest  flows  occurring  in  March  and  September, 
respectively.   High  and  low  flow  monthly  averages,  over  a 
43-year  recording  period,  were  52,000  cfs  and  17,000  cfs, 
respectively.   The  river  is  the  major  source  of  both  coarse 
quartz  sands  and  clays.   Most  of  the  finer  sediments  are 
located  in  the  basin  of  Apalachicola  Bay  proper  (Kofoed 
and  Gorsl ine  ,  1 963) . 


Fig.  2.  Map  showing  location  of  Franklin  County  and  Apala- 
chicola  Bay  on  the  northwest  Florida  coast  and  the 
extent  of  the  Apal achi col  a  River  drainage  basin. 
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About  8%   of  the  land  area  is  owned  by  commercial 
pulp  and  paper  companies  or  is  part  of  the  Apalachicola 
National  Forest.   Most  of  the  land  is  unavailable  for  in- 
dustrial or  residential  development  in  the  foreseeable 
future.   Paper  company  activities  account  for  a  small  part 
of  the  economic  base  of  the  county  (Colberg  and  Windham, 
1965;  Rockwood,  1973).   The  local  economy  is  predominantly 
dependent  on  water  resource-based  activities  centered  on 
Apalachicola  Bay.   The  Bay  produces  90%  of  the  commercial 
oyster  harvest  in  the  state  and  supports  shrimp,  crab,  and 
finfish   industries  with  an  annual  dockside  value  of  several 
million  dollars  (Rockwood,  1973).   The  largely  undeveloped 
forest  lands,  inland  and  coastal  waters,  and  barrier  islands 
in  Franklin  County  support  a  tourist  business  and  provide 
recreational  opportunities  for  the  residents  (Florida  Tour- 
ist Study,  1970).   Agriculture  and  manufacturing  activities 
are  small  at  present.   A  large  tourist-retirement  community 
is  planned  for  one  of  the  offshore  islands.   Most  basic 
goods,  fuels,  and  services  are  purchased  from  outside  the 
county  with  money  derived  from  oyster  sales,  tourism,  and 
land  development  income.   Per  capita  income  is  low  compared 
to  other  areas  of  Florida  and  the  nation.   Population  levels 
have  changed  little  due,  in  part,  to  emigration  of  residents 
searching  for  more  productive,  stable  jobs  (Colberg  and 
Windham,  1965).   In  the  1930's  Apalachicola  was  a  boom  town 
of  the  Gulf  Coast,  second  only  to  Mobile  and  New  Orleans  in 
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size.   A  major  part  of  this  growth  was  supported  by  cotton 
trade  with  the  interior  which  reached  a  peak  during  the 
years  1835-1850.   Some  75,000  bales  of  cotton  were  shipped 
from  Apalachicola  in  1839,  and  125,000  bales  in  1843.   The 
population  of  Apalachicola  was  2,060  in  1837.   Timber,  rice 
and  sugar  cane  production  also  played  a  role  in  the  economy 
during  the  territorial  period.   No  mention  was  made  of  com- 
mercial fisheries  during  this  period  (Martin,  1944).   Fol- 
lowing the  Civil  War,  lumbering  (primarily  cypress)  was  a 
major  local  activity  with  oystering,  fin  fishing,  and  boat 
building  also  mentioned  (Sawyer,  1974).   Severe  hurricanes 
hit  the  county  in  1873  and  1898,  causing  considerable  damage 
A  fire  burned  most  of  the  business  district  in  1900.   The 
population  of  the  county  at  that  time  was  3,077  (Sawyer, 
1974).. 

Following  the  turn  of  the  century,  slow  growth,  un- 
employment and  underemployment  were  characteristic  of  the 
county.   Shipping  declined  on  the  Apalachicola  River  due 
to  a  shallow  channel  to  Columbus,  Georgia  and  development 
of  more  dependable  rail  transportation.   Oysters  and  fish 
were  mentioned  after  1900  as  important  products  of  the  area 
(Franklin  County  Overall  Economic  Development  Program, 
1965).   The  current  economy  retains  i  ts  dependence  on  fish- 
ery products  (oysters,  shrimp,  and  fin  fish)  with  additional 
income  generated  from  tourism,  government  expenditures,  and 
some  outside  investment  in  retirement  and  tourist  facilities. 


A  summary  of  county  trends  including  fuel  usage,  fishery 
yields  and  value,  population,  tourism,  and  income  are  given 
in  Figs.  30  through  34  in  the  Results  section.   A  simpli- 
fied model  suggesting  relationships  within  the  county  is 
given  in  Fig.  3.   This  diagram  shows  major  inputs,  aggre- 
gated storages  of  urban,  plankton  and  oyster  reef  systems 
and  depreciation.   Figure  28a  shows  the  model  with  the  com- 
plexity used  in  simulations. 

Previous  Studies  in  Franklin  County 


Early  studies  conducted  in  Franklin  County  were 
primarily  concerned  with  the  biology  of  the  local  oyster 
fishery.   Several  economic  and  planning  studies  have  been 
recently  completed.   Early  studies  by  Ingersoll  (1881), 
Swift  (1897),  and  Danglade  (1917)  discussed  oyster  reef 
position,  density  of  oysters  and  associated  oyster  reef 
fauna.   Ingle  and  Dawson  (1953)  conducted  a  survey  of  oy- 
ster reefs  including  information  on  spawning,  setting, 
growth,  and  condition  of  oysters  (Ingle,  1951;  Ingle  and 
Dawson,  1950,  1952).   Menzel  et  al.  (1957,  1966)  published 
data  on  oyster  abundance  in  relation  to  oyster  predation 
as  influenced  by  salinity.   Menzel  and  Cake  (1969)  re- 
ported data  on  temperature,  salinity,  diversity,  and  bio- 
mass  of  invertebrates  and  fish,  and  commercial  catch  for 
important  species.   Livingston  and  Thompson  (1974)  re- 
ported pesticide  concentrations  and  other  characteristics 
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Of  the  bay  (Livingston  e_t  al.  ,  1974).   Information  on 
nutrient  stocks,  phytopl ankton  biomass  and  productivity 
was   recently  developed  by  Estabrook  (1973).   Cox  (1970) 
and  Cox  and  Autfi  (1971)  compiled  information  on  water  qual- 
ity and  fish  populations  in  the  Apalachicola  River.   Bark- 
uloo  (1961)  has  reported  on  the  striped  bass  population  in 
the  river. 

Hydrographic  studies  were  conducted  by  Dawson  (1955a) 
and  Gorsline  (1963).   A  single  layer  spatial  simulation  model 
of  salinity  and  other  conservative  materials  was  produced 
by  Swallows  (1973).   Estabrook  (1973)  also  discussed  hydro- 
graphic  features  in  the  bay. 

Jordan  (1951)  described  the  geology  of  the  conti- 
nental slope  off  the  Apalachicola  River  and  reef  formation 
in  the  adjacent  Gulf  of  Mexico  (Jordan,  1952).   Kofoed  and 
Gorsline  (1963)  published  maps  of  the  sedimentary  geology 
of  Apalachicola  Bay.   Summaries  of  water  quantity  informa- 
tion, fishery  yields,  river  flow,  cl imatol ogical  data,  and 
coastal  land  uses  were  given  in  McNulty  e_t  aj_.  (1972)  and 
Jones  let  aj^.  (1973).   A  statistical  summary  and  maps  of 
Florida  coastal  counties  were  produced  by  the  Florida  Coastal 
Coordinating  Council  (1970). 

Economic  studies  emphasized  the  importance  of  the 
oyster  industry  to  the  local  economy  and  included  those  by 
Cclberg  and  Windham  (1965),  Colberg,  Dietrich,  and  Windham 
(1968),  and  Rockwood  (1973).   The  Rockwood  report  also 
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included  suggestions  for  oyster  industry  management.   Whit- 
field (1973)  described  the  success  of  construction  and  re- 
habilitation of  oyster  reefs  in  the  county  and  elsewhere. 
He  also  presented  estimates  of  the  value  of  submerged  lands 
based  on  present  and  projected  dollar  sales  of  marine  products 

Regional  studies  include  those  which  considered  eco- 
nomic and  industrial  development  in  Franklin  County  (Frank- 
lin County  Planning  and  Development  Committee,  1965),  com- 
mercial tourism  and  land  absorption  (Northwest  Florida  De- 
velopment Council,  1972),  and  an  economic  development  pro- 
file of  Apalachicola ,  Florida  (Florida  Dept.  of  Commerce, 
1973).   Recently,  the  RMBR  Planning/Design  Group  (1974), 
under  a  federal  grant,  produced  a  county  development  plan 
including  information  on  social,  economic,  and  environmental 
conditions.   Included  in  the  report  was  a  proposed  set  of 
land  development  regulations.   Hawkins  (1973)  and  Boynton 
ejt  aj^.  (1975)  developed  mathematical  simulation  models  show- 
ing relationships  between  the  estuarine  oyster  industry  and 
coastal  development.   The  Florida  Statistical  Abstract 
(1973)  contains  detailed  social,  cl imatol ogi cal  and  economic 
data  by  both  state  and  county  areas.   The  Florida  Division 
of  State  Planning  has  developed  a  bibliography  for  Apala- 
chicola Bay  and  the  Apalachicola  River  basin. 

Issues  for  Decision  in  Franklin  County 


From  discussions  with  local  leaders  in  Franklin 
County,  synthesis  of  data  collected  in  previous  and  ongoing 


studies  in  the  county,  and  from  general  principles  of  sys- 
tems ecology  (Odum,  1971a,  1973)  the  following  list  of 
questions  needing  understanding  and  decision  was  developed. 

(1)  What  is  the  energy  value  of  a  generally  higher  more 
stable  salinity  regime  on  an  oyster  producing  estuary 
adapted  to  seasonal  pulses  in  river  flow,  sunlight, 
and  temperature? 

(2)  What  is  the  value  of  proposed  coastal  island  develop- 
ment as  a  part  of  a  coastal  fishing  economy? 

(3)  What  are  the  prospects  for  future  growth  in  Franklin 
County  and  its  impact  on  the  oyster  industry? 

(4)  What  are  estuarine  responses  to  nutrient  and  detritus 
loading  from  river  and  local  sources? 

(5)  What  are  the  major  energy  sources  characterizing  Apala- 
chicola  Bay  estuarine  ecosystems  and  oyster  fishery? 


METHODS 

The  evaluation  of  the  energy  basis  of  Franklin 
County  included  system  models,  direct  field  measurements, 
energy  calculations,  and  aerial  maps.   Energy  flows  and 
storages  were  identified  and  their  interactions  charac- 
terized with  energy  circuit  diagrams.   The  model  diagrams 
indicated  critical  data  to  be  collected  in  the  field  and 
were  the  basis  for  simulations  and  energy  value  calcula- 
tions.  Field  measurements  were  made  of  total  metabolism, 
seasonal  salinity  patterns,  and  organic  matter  inputs  and 
stocks.   Nutrient  flows  and  stocks,  fishery  stocks  and 
yields,  and  county-wide  money  flows  were  available  from 
other  studies.   Maps  were  used  to  quantify  areas  of  urban 
and  natural  land-uses  and  to  show  changing  land-use  patterns 
Computer  simulations  of  simplified  energy  models  were  used 
to  test  our  concepts  of  interactions  of  processes  as  com- 
pared to  field  data.   After  sensitivity  checks  and  valida- 
tion, simulation  results  suggested  future  county  trends. 
Energy  calculations  were  made  to  determine  which  of  several 
development  alternatives  now  facing  the  county  has  the 
greatest  survival  contribution. 
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Model i  ng  Techni ques 


From  literature  reviews,  assembled  field  data,  and 
discussions  with  individuals  familiar  with  local  processes 
in  the  study  area,  an  initial  list  of  model  components  was 
prepared.   The  components  included  all  storages  of  materials 
in  the  system,  external  forcing  functions  and  pathways  of 
interaction  identified  as  important  to  the  study  area.   In 
translating  the  component  list  to  diagram  form,  storages 
or  state  variables  were  shown  by  the  tank  symbol,  and  forc- 
ing functions  or  external  energy  sources  which  affect  amounts 
in  storage  by  the  circle  symbol.   Transfer  relationships 
between  storages  and  between  storages  and  forcing  functions 
were  shown  with  solid  lines.   Money  flows  were  shown  as 
dashed  lines  and  flowed  counter-current  to  material  and 
energy  flows.   Interactions  were  shown  with  the  workgate 
symbol  and  were  some  function  of  the  interacting  pathways. 
Intersection  functions  may  be  additive,  multiplicative, 
logarithmic,  exponential,  or  switching.   An  explanation  of 
energy  circuit  language  symbols  is  given  in  Appendix  A. 
After  further  discussions  and  literature  reviews  an  im- 
proved diagram  was  made  summarizing  knowledge  about  form 
and  function  within  the  study  area. 

In  modeling  efforts  selection  of  system  boundaries 
determines  which  components  are  state  variables,  and  which 
are  forcing  functions  and  independent  of  model  behavior. 
If  modeling  boundaries  are  chosen  one  scale  larger  than 
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the  size  scale  containing  the  specific  component  of  inter- 
est, there  are  fewer  external  driving  functions.   For  ex- 
ample, if  understanding  oyster  production  is  a  modeling 
goal,  modeling  would  be  done  at  the  scale  of  the  estuary. 
In  this  way,  more  components  subject  to  change  or  inter- 
action are  contained  within  the  model.   The  procedure  avoids 
the  mistake  of  showing  a  storage,  subject  to  change,  as  a 
constant  source  independent  of  changes  within  the  model. 

When  the  full  model  was  diagramed,  numerical  data 
were  obtained  on  observed  chronological  records.   The  model 
was  evaluated  and  simulated  with  initial  values  for  an  ini- 
tial time  period.   In  cases  where  local  information  was  not 
available  for  evaluating  a  portion  of  the  model,  data  were 
taken  from  the  same  process  in  a  similar  area.   In  some  por- 
tions of  the  model  where  no  data  were  available,  some  rates 
and  storage  values  were  calculated  by  difference  as  if  the 
system  was  in  a  temporary  steady  state. 

When  all  components  of  the  system  have  been  evalu- 
ated numerically,  the  relative  sizes  of  storage  and  flows 
become  evident.   Indications  of  the  importance  of  flows 
and  the  stability  of  state  variables  can  be  obtained  by 
calculating  the  residence  time  of  materials  in  each  storage 
(calculated  by  dividing  the  amount  in  storage  by  the  total 
inflow).   Flows  that  are  of  a  similar  nature  are  compared. 
Those  that  are  several  orders  of  magnitude  larger  than  others 
are  carefully  considered  because  they  may  have  controlling 
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actions  in  the  model.   Flows  with  effects  several  orders 
of  magnitude  smaller  than  identical  flows  can  be  omitted. 

Simplification  of  Simulation  Models 

Fully  evaluated  models  were  simplified  prior  to 
simulation  for  several  reasons.   First,  in  the  evaluated 
model  many  pathways  were  small  and  thus  not  considered  to 
be  essential.   Rather  than  model  these  flows  by  themselves, 
they  were  either  lumped  with  similar,  larger  flows  or  ne- 
glected.  The  same  procedure  was  applied  to  state  variables. 
Thus  some  complexity  was  eliminated.   Second,  some  factors 
included  in  the  evaluated  model  may  not  change  or  limit 
any  process.   These  factors  can  be  eliminated  from  the 
models  while  their  effect  remains  implicit  in  the  calcula- 
tion of  transfer  coefficients.   Comparison  of  simulation 
results  from  condensed  models  with  field  and  literature 
data  provide  indications  of  model  realism. 

Writing  and  Scaling  Equations 

The  set  of  differential  equations  associated  with 
each  model  is  taken  directly  from  the  energy  language  dia- 
grams.  Scaled  terms  for  each  equation  are  given  in  Appen- 
dix Tables  C-3  and  D-3.   Odum  (1972a)  gives  theoretical 
derivations  of  mathematical  equivalents  of  the  energy  symbols 

In  analog  simulations  real  world  values  contained 
in  equations  must  be  magnitude-scaled  so  as  not  to  exceed 


the  voltages  available  in  the  analog  computer.   To  do  this, 
transfer  coefficients  must  have  scaled  values  between  0.001 
and  0.999,  and  state  variables  and  forcing  functions  must 
have  scaled  values  as  a  percent  of  full  scale  of  the  com- 
puter.  Each  state  variable  and  forcing  function  is  assigned 
a  maximum  value.   This  number  is  at  least  as  large  as  the 
maximum  value  expected  for  that  variable  but  preferably 
two  to  ten  times  larger  so  as  to  avoid  overloads  in  the 
early  modeling  stages.   The  calculated  transfer  coefficient 
associated  with  each  term  in  each  equation  is  multiplied  by 
the  maximum  value  (or  product  of  the  maximum  values  if  the 
term  is  multiplicative)  associated  with  that  term.   This 
product  is  then  divided  by  the  maximum  value  of  the  state 
variable  which  the  full  equation  describes.   In  this  fashion, 
proportionality  of  all  terms  is  maintained  while  each  term 
is  adjusted  to  stay  within  the  capacity  of  the  machine.   If 
after  this  procedure  is  completed,  scaled  transfer  coeffi- 
cients are  less  than  0.001  or  exceed  0.999,  time  scaling 
may  be  required.   This  procedure  increases  or  decreases  the 
time  period  used  to  calculate  flows.   For  instance,  if  a 
scaled  model  had  flows  calculated  in  units  of  per  year,  di- 
vision of  transfer  coefficients  by  a  factor  of  10  would 
reduce  the  transfer  coefficient  value  by  a  factor  of  10 
and  flows  to  units  of  per  0.1  years.   This  discussion  was 
intended  to  provide  an  overview  of  the  scaling  process  re- 
quired in  analog  computation.   Details  of  magnitude  and 
time  scaling  are  given  in  Peterson  (1967). 
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Simulation  Techniques 

The  regional  model  was  simulated  on  two  slaved 
Electronic  Associates,  Incorporated  Model  580  computers. 
The  oxygen  balance  model  was  simulated  on  two  slaved  Elec- 
tronic Associates,  Incorporated  MiniAc  analog  computers. 
For  each  model  an  analog  patching  diagram  was  developed 
directly  from  scaled  differential  equations  and  patched 
onto  appropriate  analog  boards.   Analog  patching  diagrams 
are  given  in  Appendix  Figs.  C-1,  D-1,  and  D-2. 

Following  checks  to  be  sure  all  transl ati onal  pro- 
cedures to  this  point  were  correct,  each  model  was  tested 
for  responses  to  expected  changes  in  the  study  areas.   Each 
computer  run  was  similar  to  a  controlled  experiment  in 
which  one  or  more  factors  were  changed.   Expected  changes 
were  first  calculated  and  then  translated  into  new  computer 
settings  and  run.   Model  output  was  either  photographed 
directly  from  the  oscilloscope  attached  to  the  computer  or 
obtained  from  a  plotter. 

Validation  and  Basis  for  Model  Configurations 


The  following  procedures  were  used  in  validating 
models.   For  some  variables  and  some  rates,  graphs  of  his- 
torical data  were  assembled  and  compared  to  model  output 
with  the  model  being  started  in  the  past.   Good  agreement 
between  model  output  and  graphs  indicates  an  absence  of 


25 


contradictions.   This  procedure  was  especially  useful  wfien 
systems  had  undergone  appreciable  changes  in  the  past. 
Simulations  of  natural  system  models  were  validated  with 
time  graphs  obtained  from  field  studies  in  Apalachicola  Bay. 
microcosm  studies  and  field  studies  in  similar  areas.   In  a 
second  validation  technique,  model  sensitivities  were  com- 
pared with  single  and  multiple  flow  rate  test  changes. 

This  section  also  provides  a  rationale  for  the  use 
of  the  energy  symbols  in  the  models  presented  in  this  dis- 
sertation.  Given  here  are  the  criteria  concerning  the 
selection  of  state  variables,  forcing  functions,  pathways, 
and  interactions.   Storages  (state  variables)  were  used  to 
define  functional  groupings  within  the  system  which  were 
expected  to  vary  with  time  and  which  had  interaction  with 
other  components  or  pathways  in  the  model.   Everything  in- 
cluded in  a  tank  had  characteristics  which  were  similar  for 
the  purposes  of  the  simulation  and  could  be  considered  as 
a  whole.   For  instance,  the  natural  land  tank  in  the  re- 
gional model  (Fig.  28a)  contained  a  variety  of  ecosystem 
types  (see  Appendix  B).   When  used  for  purposes  other  than 
simulation  this  storage  was  broken  down  into  8  specific 
ecosystems  (see  Figs.  6  and  7).   Each  of  those  could  be 
further  subdivided  if  a  finer  level  of  detail  was  required. 
However,  each  of  the  8  ecosystems  contained  in  the  natural 
land  storage  had  the  properties  of  occupying  space  in  the 
county  and  self-maintenance  without  urban  inputs. 
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Obviously,  not  all  storages  of  energy  within  the 
boundaries  of  each  model  have  been  included.   The  same  is 
true  for  forcing  functions  and  pathv/ays.   Those  which  were 
expected  to  change  were  identified  and  retained  in  the  model, 
Those  which  were  very    large  or  not  expected  to  change  and 
those  which  were  small  or  had  few  or  no  major  interactions 
in  the  county  were  deleted.   A  similar  process  was  followed 
in  arranging  pathways  and  forcing  functions. 

Three  types  of  interactions  were  used  in  the  models 
and  included  multiplicative  and  drag  action  workgates  and 
workgates  operated  by  a  force  from  a  sensor.   Multipliers 
were  used  extensively  because  they  had  the  property  of  re- 
quiring all  inputs  to  be  present  for  a  flow  to  occur.   Of- 
ten it  was  possible  to  identify  factors  required  for  a  flow 
to  occur  and  thus  multipliers  were  used  to  show  this  require- 
ment.  Secondly,  the  flow  produced  from  this  type  of  inter- 
action is  proportional  to  the  levels  in  the  input  storages. 
If  one  or  both  increased,  then  the  flow  increased  until 
one  became  limiting.   This  property  was  useful  because  it 
was  often  possible  to  identify  such  relationships.   In  some 
cases,  as  in  primary  production,  the  response  of  a  multi- 
plier is  of  the  limiting  factor  type  which  is  a  well-estab- 
lished photosynthetic  response  (Odum,  1971a;  Odum,  1972a; 
Rabinowitch  and  Govindjee,  1969). 

Drag  action  workgates  ware  used  in  cases  where  one 
flow  was  inhibited  by  increases  in  one  state  variable,  but 
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increased  with  increases  in  another.   These  were  useful  in 
approximating  situations  having  both  stimulating  and  in- 
hibiting factors  interacting  and  producing  a  flow.   For 
instance,  in  the  regional  model  shown  in  Fig.  28a,  juvenile 
oysters  feed  in  proportion  to  the  food  available  but  are 
inhibited  when  salinity  is  yery    low.   If  salinity  reached 
zero,  there  would  be  no  food  intake  by  the  juvenile  oysters. 
While  the  exact  nature  of  the  feeding-salinity  relationship 
was  not  known,  this  interaction  captured  the  basic  features. 
Lastly,  workgates  operated  by  a  force  from  a  sensor 
were  used  in  all  cases  where  energy  with  or  without  mater- 
ials flowed  in  proportion  to  some  other  flow.   All  river  in- 
puts were  done  this  way  as  well  as  inflows  of  capital  asso- 
ciated with  incoming  tourists  and  new  residents  as  shown  in 
the  regional  model  (Fig.  28a). 

Field  Study  Techniques 

Field  study  techniques  used  here  included  measure- 
ments of  community  metabolism  in  the  bay,  particulate  carbon 
concentrations  in  the  bay  and  river,  and  two  synoptic  salin- 
ity surveys  in  the  bay. 

Metabolic  Measurements 


Community  metabolism  was  measured  using  the  diurnal 
oxygen  sampling  method  of  Odum  and  Hoskins  (1958)  and  an 
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abbreviated  dawn-dusk-dawn  method  used  by  McConnell  (1952). 
Oxygen  was  measured  using  the  azide  modification  of  the 
Winkler  technique  (Strickland  and  Parsons,  1972)  with  a 
further  modification  of  thiosulfate  normality  which  allowed 
convenient  use  of  small  sampling  bottles.   Details  of  the 
procedure  used  here  are  given  by  McKellar  (1975). 

Between  June  22  and  September  30,  1974,  36  diurnal 
stations  were  occupied  with  stations  AA,  A,  C,  D,  F,  H,  L 
each  occupied  twice  (Fig.  4).   The  sampling  scheme  was  aimed 
at  generally  characterizing  the  magnitude  of  metabolism  and 
the  photosynthesis-respiration  (P-R)  balance  in  the  dominant 
bay  environments.   Stations  were  chosen  to  represent  the 
three  dominant  ecosystem  types  existing  in  the  bay.   Sta- 
tions J,  K,  L,  M,  Q,  Z,  Y  and  the  boat  basin  station  at  the 
Florida  State  University  Lab  were  established  in  areas 
where  respiration  predominated.   In  the  medium  salinity 
portion  of  the  bay  system  stations  AA,  A,  B,  C,  D,  E,  F,  H 
were  established  where  phytopl ankton  activity  predominated 
and  most  oyster  bars  were  located.   Along  the  shorelines  and 
clearer  areas  where  benthic  grasses  predominated  stations 
were  established  at  the  Florida  State  Lab  and  station  N. 
Two  stations  were  also  taken  where  wastes  were  suspected 
of  entering  the  bay  (P  and  I). 

Over  the  course  of  any  one  diurnal,  stations  were 
sampled  from  3  to  7  times  depending  on  travel  time  between 
stations  and  weather  conditions.   The  minimum  number  of 
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samples  required  for  a  complete  diurnal  was  a  dawn-dusk- 
dawn  sample  sequence,   "^ach  sampling  included  measurements 
of  oxygen,  salinity,  temperature,  and  depth. 

Diurnal  oxygen  data  were  plotted  and  calculated  as 
shown  by  Odum  and  Hoskins  (1958),  except  for  the  diffusion 
adjustment  to  the  final  rate  of  change  curve.   One  measure- 
ment of  diffusion  between  the  water  and  the  air  was  measured 
using  a  small  nitrogen-filled  plastic  dome  which  floated  on 
the  surface  (Hall,  1972,  based  on  original  work  by  Copeland 
and  Duffer,  1964).   The  diffusion  rate  as  g  02/ni^/hr  at  100% 
deficit  was  calculated  from  the  area  and  volume  of  the  dome, 
ambient  oxygen  concentration  in  the  water,  and  the  rate  of 
oxygen  increase  in  the  dome.   Oxygen  changes  were  measured 
using  a  membrane  oxygen  probe.   The  calculation  of  a  diffu- 
sion coefficient  followed  the  method  used  by  McKellar  (1975), 
In  calculating  rate  of  change  curves,  diffusion  coefficients 
were  estimated  from  a  graph  that  included  the  one  measure- 
ment made  in  Apalachicola  Bay  and  measurements  reported  by 
McKellar  (1975)  for  Crystal  River,  Florida. 


Salinity  Measurements 


In  September  of  1973  and  April  of  1974  synoptic 
salinity  surveys  were  made  in  Apalachicola  Bay.   Stations 
are  shown  as  dots  in  Figs.  8  and  9.   Each  survey  was  started 
about  two  hours  before  high  tide  and  was  completed  about  2 
hours  after  high  tide.   In  each  survey  top  and  bottom 
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salinity  and  temperature  measurements  were  taken  at  each 
station  using  a  Beckman  electrodeless  induction  salinometer. 

Dissolved  and  Particulate  Carbon  Analyses 

In  the  summer  of  1973  one  effort  was  made  to  esti- 
mate the  standing  .stock  of  particulate  and  dissolved  organic 
carbon  in  the  water  column  and  in  the  river  for  use  in  the 
simulation  models.   Sixteen  stations  were  sampled  including 
two  stations  in  the  Apalachicola  River  located  approximately 
1  and  2  miles  upstream  from  the  mouth.   Water  samples  were 
taken  approximately  1  meter  from  the  surface  for  all  bay 
stations  and  approximately  1  meter  from  the  top  and  bottom 
of  the  river.   Samples  were  stored  in  acid-washed  polyethe- 
lene  bottles  and  preserved  immediately  with  2-3  drops  of 
saturated  mercuric  chloride  solution  per  100  ml  of  sample. 
Analyses  for  total  organic  carbon  were  made  from  unfiltered 
water.   Particulate  carbon  measurements  were  made  with 
material  retained  on  precombusted  0.45  ym  glass  fiber  fil- 
ters.  Dissolved  organic  carbon  values  were  obtained  by 
difference.   Subsequent  techniques  in  the  analysis  followed 
those  specified  for  use  with  the  Oceanography  International 
Total  Carbon  System  (1972). 


Energy  Value  Calculations 

Energy  value  calculations  were  made  to  quantify 
total  work  contributions  from  various  components  of  man  and 
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nature.   Value  was  calculated  as  the  work  done  for  the 
total  economy  of  man  and  nature  when  all  work  processes 
were  expressed  as  work  equivalents  of  one  kind  of  energy. 
The  procedure  takes  as  its  theoretical  base  Lotka's  maximum 
power  principle  (Lotka,  1922)  as  developed  by  Odum  (1971a). 
This  principle  states  that  the  system  that  tends  to  prevail 
over  alternative  systems  is  the  one  that  maximizes  the  use 
of  all  energy  flows  available  to  it  and  is  then  one  which 
develops  useful  feedback  roles  for  all  participants  for 
the  purpose  of  assuring  continued  energy  flows  and  for  cap- 
turing any  additional  flows  that  come  available.   In  a  sys- 
tem containing  both  man  and  nature  (towns,  estuaries,  fish- 
eries, river,  industry,  and  forests),  the  principle  requires 
using  the  variety  of  components  in  partnerships  rather  than 
in  competition  to  obtain  maximum  work  per  time  (power) 
through  the  full  system.   Two  additional  concepts  were  used 
in  energy  calculations  and  are  presented  next,  followed  by 
the  general  procedure. 

Energy  Quality  Concept 


The  value  of  a  system  process  was  defined  as  the 
contribution  of  the  process  to  the  useful  work  of  the  sys- 
tem.  In  evaluating  energy  models,  all  pathways  were  first 
shown  in  units  of  ki 1  oca! ori es/area/t ime ,   Raw  energy  flows, 
as  measured  in  ki 1 ocal ories  of  heat,  do  not  represent  the 
ability  to  do  work  but  rather  show  only  the  heat  content  of 
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that  particular  flow.   Whereas  any  energy  flow  can  be  de- 
graded to  heat  with  100%  efficiency  the  ability  of  an  energy 
flow  to  do  useful  work  depends  on  the  packaging  or  concen- 
tration of  the  energy  flow.   For  instance,  the  kilocalories 
associated  with  wood  production  in  photosynthesis  represent 
the  concentration  to  wood  kilocalories  of  the  dilute  kilo- 
calories  of  unprocessed  sunlight  and  energy  of  collecting 
wood  to  one  location.   Approximately  99%  of  the  sunlight 
kilocalories  are  degraded  to  heat  to  make  one  kilocalorie 
of  wood  quality.   While  it  is  difficult  for  human  systems 
to  use  the  dilute  energy  contained  in  sunlight  directly, 
the  energy  contained  in  wood  has  been  concentrated  or  up- 
graded at  the  expense  of  sunlight  energy  and  thus  represents 
a  higher  quality  energy.   In  the  same  way  electrical  energy 
is  at  higher  concentration  than  the  energy  contained  in 
coal;  its  generation  requiris   approximately  four  kilocal- 
ories of  coal  type  energy  to  obtain  one  unit  of  electrical 
type  energy.   Three  kilocalories  of  coal  energy  are  contrib- 
uted from  the  coal  to  operate  steam  engines  and  one  kilo- 
calorie  is  expended  to  perform  the  work  of  constructing 
arid  maintaining  the  power  plant  structure.   Several  other 
examples  of  conversion  calculations  are  given  by  Odum  ejt 
a_l_.  (1974b)  for  relating  kilocalories  of  wind»  wood,  and 
electricity.   Ene^-gy  quality  relationships  relating  pro- 
ducers and  many  consumers  in  a  shallow  marine  ecosystem 
were  given  by  McKellar  (1975).   One  calculation  of  the 
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energy  quality  of  the  potential  energy  of  river  head  being 
converted  to  electricity  was  given  by  Young  ejt  aj_.  (1974). 
An  additional  calculation  of  river  head-electrical  relations 
is  contained  in  this  dissertation  for  the  dam  on  the  Apalachicola 
River.   Thus,  the  above  considerations  suggest  ways  to  com- 
pare varying  types  of  energy  flows  in  macroscopic  systems 
of  man  and  nature.   Before  comparisons  are  made,  each  flow 
must  be  converted  to  a  common  baseline  energy  quality.   In 
this  dissertation,  all  energy  flows  have  been  converted  to 
the  fossil  fuel  quality  level  (expressed  as  Kcalnpc).  A 
complete  list  of  conversion  factors  used  in  this  disserta- 
tion is  given  in  Table  1. 

In  this  study,  the  following  method  is  used  to  de- 
termine the  energy  quality  of  a  work  flow  relative  to  other 
types  of  flows.   The  method  involves  calculating  the  energy 
required  to  generate  the  high-quality  work  flow.   As  shown 
in  Fig.  5,  the  quality  factor  is  obtained  by  comparing  the 
output  of  the  conversion  process  with  the  input  energy  af- 
ter all  required  feedbacks  and  subsidies  have  been  converted 
to  equivalents  of  the  same  quality  and  subtracted.   All 
output  feedback  and  subsidy  energies  must  be  of  the  same 
quality  before  being  manipulated.   A  complete  calculation 
of  this  type  is  given  in  Table  13,  where  potential  energy 
of  river  head  is  converted  to  fossil  fuel  and  electrical 
work  equivalents.   This  method  assumes  that  maximum  energy 
conversions  done  through  one  pathway  equal   those  done 
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Table  1 
Energy  quality  factors  relating  different  work  processes. 


Energy  Conversion  Process 


Energy  Quality  Factor* 


Sunlight  to  gross  photosynthetic 
production 

Gross  photosynthetic  production 
to  wood 

Wood  to  fossil  fuel 

Wood  to  electricity 

Fossil  fuels  to  electricity 

Gross  photosynthetic  production 
to  fossil  fuel 

x-Sunlight  to  fossil  fuel 

''   Wind  to  fossil  fuel 

^  Tidal  energy  to  fossil  fuel 

Hydrostatic  head  to  fossil  fuel 

Fresh/salt  water  concentration 

Heat  gradient  to  fossil  fuels 

y   Waves  to  fossil  fuels 

Fossil  fuel  flow  to  dollar  flow 
in  U.S. 


/ 


(divide  by) 
100 

10 
2 
8 
3.6 

20 
2000 
7.7 
2.5 
0.63^ 
0.3 
"  10,000 
5 

25,000  Kcal/dollar 


Energy  quality  factor  is  the  ratio  of  energy  input 
of  one  type  required  to  generate  another  type  of  higher  qual- 
ity.  (See  Methods  section)   Energy  quality  factors  used  here 
are  preliminary  and  may  be  subject  to  revisions.   For  calcu- 
lations leadina  to  ratios  given  above  see  Odum  e_t  aj_.  (1974b), 
and  Young  e^  al.  (1974). 

Calculations  for  conversion  of  hydrostatic  head  to 
electricity  and  fossil  fuel  are  given  in  the  Results  section 
of  this  dissertation. 


Fig.  5.   Diagram  showing  example  of  the  upgrade  method  of 
calculating  energy  quality  ratios.   Note  that  a 
portion  of  the  gross  output  is  fed  back  to  the 
process,  via  the  general  economy,  to  cover  full 
cost  of  operating  the  process.   This  is  required 
to  get  a  true  net  relationship  between  X  and  Y 
forms.   The  eneray  quality  ratio  relating  X  to 
Y  is  20  in  this  example  (100/5). 
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SUBSIDIES 
REQUIRED 

FOR 
PROCESS 


GOODS    AND   SERVICES 
2  units/doy; 


X-QUALITY 
SOURCE    VOO""'ts/day 
INPUT 


GENERAL 
ECONOMY 


FUELS 
3  units/day 


OTHER 
EXCHANGES 


/  5  units/day 


10  units/doy 


Y- QUALITY 
5  units/day 


GROSS    OUTPUT 


NET  OUTPUT 


95  units/day 
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through  another  if  the  systems  are  under  competition  for 
survival  and  selected  to  operate  at  maximum  power. 

Thus,  the  application  of  energy  quality  ratios  to 
all  energy  flows  allows  for  equal  evaluation  of  all  types 
of  work  whether  they  be  associated  with  man  and  money  flows 
or  supplied  without  money  exchanges. 


Energy  Investment  Ratio  Theory 


In  some  energy  evaluations,  total  energy  flow  fol- 
lowing a  proposed  alteration  is  greater  than  before  the 
alteration.   By  criteria  of  maximum  power,  the  new  pattern 
would  seem  better  than  the  old  and  should  be  adopted.   How- 
ever, outside  goods,  services,  and  fuels   have  to  be  pur- 
chased, and  the  income  and  investments  for  these  purchases 
can  only  be  assured  if  the  resident  natural  energies  are 
enough  to  attract  the  required  capital.   The  energy  invest- 
ment ratio  is  the  ratio  of  high  quality  external  energies 
that  are  attracted  from  outside  the  system  and  paid  for  by 
some  kind  of  exchange  to  the  natural  energies  operating  in 
the  study  area.   In  1973  this  ratio  was  estimated  to  be  2.5 
to  1  for  the  United  States  (Odum  £t  al.  ,  1974c).   Systems 
that  have  lower  investment  ratios  can  match  high  quality 
external  energies  with  more  lower  quality  natural  energies 
and  thus  compete  well  with  those  things  offered  for  ex- 
change.  The  theory  suggests  that  as  the  local  investment 
ratio  exceeds  the  ratio  of  surrounding  systems,  the  local 
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system  generates  less  value  per  unit  of  high-quality  energy 
used.   This  decline  in  value  per  unit  could  be  reflected 
in  higher  prices  required  for  exports  and  thus  a  disadvan- 
tage in  competing  with  other,  less  developed  systems  for 
high-quality  energies.   The  ultimate  contribution  of  energy 
flows  depends  both  on  high-quality  purchased  flows  of  fossil 
fuels  and  resident  natural  energies  with  which  the  high-qual 
ity  flows  interact. 

Energy  Evaluation  of  a  Process  or  Alternative 


The  first  step  in  the  energy  evaluation  of  a  proc- 
ess or  an  alternative  is  to  construct  a  diagram  using  the 
energy  symbol  language  including  all  major  work  processes 
occurring  in  the  study  area.   These  diagrams  include  the 
work  done  by  nature  in  maintaining  soils,  water  flows, 
forests,  estuarine  and  river  systems,  beaches,  etc.,  as 
well  as  the  work  associated  with  man's  activities.   By 
putting  the  work  processes  of  both  man  and  nature  on  the 
same  diagram,  we  summarize  all  of  the  work  processes  that 
contribute  to  the  welfare  of  a  region.   The  summary  dia- 
gram avoids  the  common  mistake  of  counting  as  valuable 
only  those  processes  associated  with  money  (Odum  197Zb, 
1973). 

The  next  step  is  to  evaluate  all  pathways  in  the 
diagram  in  units  of  work  ( Kcal /area/yr) .   In  this  disser- 
tation, natural  system  work  was  evaluated  by  using  gross 
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community  metabolism  as  an  estimate  of  total  work.   Work 
done  by  physical  systems  such  as  tides,  waves,  wind,  etc., 
was  calculated  by  standard  formulas.   Work  done  by  human 
activity  was  often  available  as  dollar  expenditures  and 
converted  to  work  units  in  the  following  way.   In  the  U.  S. 
economy  as  a  whole,  there  is  an  average  exchange  ratio  be- 
tween work  done  and  money  flow  (Kcal /dol lar) .   After  sev- 
eral previous  calculations  (Odum  ejt  aj^.  ,  1972),  this  ratio 
in  the  United  States  was  estimated  to  be  approximately 
25,000  Kcalppg/dollar  in  1973  (Odum  and  Brown,  1975).   This 
figure  was  obtained  by  summing  total  fuel  usage  in  the  U.S. 
per  year  and  total  work  of  the  natural  systems  expressed 
as  Kcalppc  required  to  do  the  same  work.   This  figure  was 
then  divided  by  the  gross  national  product  for  that  year 
giving  an  average  ratio  of  work  that  accompanies  a  dollar 
expenditure  expressed  in  ki localories/dol lar.   Thus,  when 
data  on  human  acti vi ty  are  gi ven  in  dollars,  the  dollar 
figure  can  be  converted  to  kilocalories  using  the  ratio 
of  25,000  Kcal/dollars  if  the  nature  of  that  work  has  energy 
uses  similar  to  the  national  average. 

When  the  energy  diagram  is  fully  evaluated  a  table 
is  constructed  with  each  pathway  in  the  diagram  becoming 
an  entry  in  the  table.   Next,  each  work  flow  is  expressed 
in  units  of  equivalent  work  by  one  type  of  energy.   Thus, 
all  energy  flows  are  converted  to  a  common  work  quality 
level.   Once  all  entries  are  converted,  they  may  be  summed 
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to  a  total  showing  total  value  generated  in  the  region  per 
year  under  the  present  pattern  of  man  and  nature.   The  in- 
vestment ratio  can  also  be  calculated.   Following  this, 
separate  tables  can  be  constructed  for  the  totally  natural 
condition  as  it  was  thought  to  once  exist  and  for  antici- 
pated outcomes  or  patterns  of  man  and  nature  under  each 
alternative.   Investment  ratios  for  each  alternative  can 
be  calculated  and  compared  to  national  and  local  averages. 

Land-Use  Maps  from  Aerial  Photographs 


Two  sets  of  maps  were  constructed  to  show  spatial 
features,  one  of  the  primitive  condition  and  one  of  the 
present  pattern  of  man  and  nature  in  the  Franklin  County 
region. 

The  primitive  map  was  constructed  using  a  soils  map 
of  the  region.   The  assumption  was  made  that  soil  types  re- 
flect original  vegetation  types.   The  scale  on  the  map  was 
1  inch  equals  4  miles. 

The  map  showing  present  land-use  was  developed  from 
aerial  photographs  taken  in  1973  on  a  scale  of  1  inch  equals 
2000  feet.   The  map  was  verified  with  checks  on  the  ground 
and  with  the  Franklin  County  Agricultural  Extension  Agent. 
Each  system  type  shown  on  the  maps  is  defined  in  Appendix  B. 


RESULTS 

Maps  of  Principal  Ecosystems  in  Franklin  County 

Shown  in  Figs.  6  and  7  are  maps  of  the  ecosystems 
in  Franklin  County.   A  description  of  each  type  of  system 
is  given  in  Appendix  B.   Figure  6  is  the  map  showing  the 
primitive  pattern  prior  to  the  first  commercial  cutting 
of  the  forest  (circa  1838).   Areas  of  each  ecosystem  in 
the  primitive  pattern  are  given  in  Table  2. 

The  present  pattern  is  shown  in  Fig.  7.   Areas  of 
present  land-use  are  given  in  Table  3.   In  the  present 
pattern  urban  land  accounts  for  only  0.5%  of  the  total 
regional  area  (about  1.2%  of  the  land  area).   Approximately 
15%  of  the  original  bay  swamp  and  most  of  the  original  pine 
palmetto  flatwoods,  high  pinelands,  and  pine  scrub  have 
been  converted  to  pine  forests.   Cleared/planted  lands 
account  for  approximately  26%  of  the  land  area,  most  of 
which  was  originally  bay  swamps  and  pine  palmetto  flat- 
woods.   The  river  swamp  area  is  approximately  15%  smaller 
than  in  the  primitive  pattern.   In  the  estuary,  new  systems 
account  for  only  1.4%  of  the  bay  area  and  include  boat 
basins  and  harbors  (200  acres),  navigation  channels  (560 
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Table  2 


Area  and  percentage  of  total  area  of  each  ecosystem 
type  in  primitive  pattern  in  Franklin  County. 


Ecosystem   Type 

Area,   Acres 

Percentage   of  Total 
Areaa 

Marine 

Oligohaline  system 

14,842 

1.65 

Grass  flats 

19,415 

2.16 

Salt  marsh 

7.701 

0.86 

Oyster  reef 

7,539 

0.84 

Sand  bar 

1,534 

0.17 

Medium  salinity 
plankton  system 

82,971 

9.25 

Coastal  plankton  system 

406,400 

45.30 

High  velocity  channels 

718 

0.08 

High  energy  beach 

(11,596') 

(1.30;) 

Coastal  dunes 

4,426 

0.49 

Freshwater 

Shallow  woodland  ponds 

1,388 

0.15 

Turbid  piedmont  river 

3,644 

0.41 

Blackwater  river 

2,396 

0.27 

Terrestrial 

Pine  scrub 

43,118 

4.81 

High  pineland 

11,586 

1.29 

Pine-palmetto  flatwoods 

111,209 

12.40 

Bay  swamp 

118,044 

13.16 

River  swamp 

35,575 

4.00 

Freshwater  marsh 

13,052 

1.45 

TOTAL 

897,154 

100.00 

^Total  county  area  includes  all  land  and  water  eco- 
systems in  the  county  and  adjacent  Gulf  of  Mexico  to  the 
10-fathom  line. 
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Table  3 


Area  and  percentage  of  total  area  of  each  ecosystem 
type  in  present  pattern  in  Franklin  County. 


Ecosystem  Type 


Area,  Acres 


Percentage  of  Total 
Area^ 


Urban 


Scraped  land 

1,134 

Residential  —  1 ight 

1  ,789 

Residential  — med. 

dense 

478 

Commercial  —  indus- 

trial 

412 

Airports 

459 

0.13 
0.20 

0.05 

0.05 
0.05 


Marine 

Oligohaline  system 

Grass  flats 

Salt  marsh 

Oyster  reefs 

Sand  bars 

Medium  salinity 
plankton  system 

Coastal  plankton 
system 

High  velocity 
channel s 

High  energy 
beaches 


14 

,841 

19 

,415 

7 

,072 

8 

,870 

1 

,534 

80 

,217 

406 

400 

767 

n 

595 

1.64 
2.14 
0.78 
1  .00 
0.17 

8.85 

44.84 

0.08 

1  .28 


Total  county  area  includes  all  land  and  water  eco- 
systems in  the  county  and  adjacent  Gulf  of  Mexico  to  the  10 
fathom  line. 
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Table  3  -  continued 


Ecosystem  Type 


Area,  Acres 


Percentage  of  Total 
Area^ 


Marine  (cont.  ) 


Coastal  dunes 

4 

,782 

Boat  basins  and 
harbors 

200 

Dredged  navigation 
channel s 

560 

Freshwater 

Shallow  woodland 
ponds 

1 

,388 

Turbid  piedmont 
River 

3 

,644 

Blackwater  river 

2 

.396 

Terrestrial 

Pine  forests  103,303 

Cleared/planted 

lands  93,015 

River  swamp  30,539 

Bay  swamp  99,494 

Freshwater 

marsh  12.100 

TOTAL  906,401 


0.53 
0.02 
0.06 

0.15 

0.40 
0.26 

11  .40 

10.26 

3.37 

10.98 

1  .33 
100.00 
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acres).  Bob  Sikes  Pass  (49  acres)  and  planted  oyster  bars 
(1331  acres). 

Characteristics  of  Apalachicola  Bay 

This  section  brings  together  data  on  Apalachicola 
Bay  collected  in  this  project  and  assembled  from  other  in- 
vestigators.  Using  these  data,  an  aggregated  model  of  the 
diurnal  properties  of  the  bay  ecosystem  was  evaluated, 
simulated,  and  compared  with  field  data. 

Salinity  Measurements 


Figures   8  and  9  are  maps  showing  the  sur- 
face and  bottom  salinity  distribution  on  September  29, 
1973  and  April  20,  1974.   Both  sets  of  measurements  were 
taken  on  a  rising  tide.   The  September  measurements  were 
taken  when  river  flow  was  low  (15,060  cfs)  while  the  April 
measurements  were  taken  during  high  flow  (49,390  cfs). 
Salinity  stratification  was  evident  over  most  portions  of 
the  bay  in  the  fall  (Fig.  8).   This  observation  is  supported 
by  data  collected  by  Estabrook  (1973).   Area  weighted  aver- 
age surface  and  bottom  total  salt  content  for  the  bay  was 

9               9 
3.58  X  10   kg  and  6.0  x  10   kg,  respectively.  These  corre- 
spond to  average  surface  and  bottom  salinities  of  15  and 
26  ppt,  respectively. 

In  the  spring  sample,  both  surface  and  bottom 

salinities  were  lower,  averaging  4.9  and  11.5,  respectively 


Fig.  8.   Salinity  pattern  in  Apalachicola  Bay  on  September  29, 
1973,  (a)  surface;  and  (b)  bottom.   River  flow  was 
15,060  cfs  and  wind  was  from  the  southeast.   Salin- 
ity stations  are  indicated  by  dots. 
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SURFACE      SALINITIES 
Sept.    29,     1973 


-WEST     PASS 


BOB     SIKES     PASS 


BOTTOM      SALINITIES 
Sopt.     29,    1973 


BOB     SIKES      PASS 


' WEST      PASS 


Fig.  9.   Salinity  pattern  in  Apalachicola  Bay  on  April  20, 
1974,  (a)  surface;  and  (b)  bottom.   River  flow  was 
49,390  cfs  and  the  wind  was  from  the  east  at  10- 
15  knots.   Salinity  stations  are  indicated  by  dots 
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SURFACE      SALINITIES 
April     20,     1974 


-WEST    PASS 


BOB    SIKES     PASS 


BOTTOM      SALINITIES 
April     20,     1974 


BOB     SIKES      PASS 


-WEST     PASS 
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(Fig.  9).   However,  the  average  relative  difference  between 
surface  and  bottom  salinity  was  greater  in  the  spring  sample 
than  in  the  fall  sample.   High  salinity  bottom  water  (20- 
30  ppt)  was  evident  in  a  large  portion  of  the  bay  along  the 
barrier  island  on  both  sides  of  Bob  Sikes  Pass.   Very  little 
stratification  was  evident  nearer  the  river  mouth  and  in 
East  Bay  in  the  spring  sample.   There  was  no  evidence  of 
thermal  stratification.   Surface  and  bottom  temperatures 
seldom  differed  by  more  than  several  degrees. 

Diurnal  Metabolism  Measurements 


Thirty-six  diurnal  metabolism  measurements  were 
taken  in  the  bay  between  June  22  and  September  30,  1973. 
Station  locations  were  given  in  Fig.  4.   Data  from  two 
24-hour  diurnal  metabolism  surveys  (Station  I,  July  24-25, 
1973,  and  Station  C,  September  29,  30,  1973)  are  shown  in 
Figs.  10a  and  10b.   In  Fig.  10a  respiration  exceeded  day- 
time net  production  (see  rate  of  change  curve)  by  a  factor 
of  four  indicating  that  a  large  portion  of  the  respiratory 
activity  was  supported  by  organic  matter  supplied  from 

other  areas.   Photosynthesis  was  at  a  maximum  at  1100  hours 

2 
(1.05  g  Op/m  /hr)  and  declined  rapidly  in  the  early  after- 

2 
noon.   Respiration  averaged  about  0.5  g  O^/m  /hr  during 

the  night.   At  this  station  the  water  was  generally  under- 
saturated  with  oxygen.   Hence,  when  the  rate  of  change 
curve  was  corrected  for  diffusion,  hourly  rates  of 
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respiration  became  larger  as  shown  on  the  corrected  rate  of 
change  curve.  Most  stations  followed  a  pattern  similar  to 
the  one  at  Station  I,  except  that  the  relative  proportion 
of  photosynthetic  and  respiratory  activity  (P/R  ratio) 
varied  from  station  to  station.   Figure  10b  shows  results 
from  a  station  (C)  where  photosynthesis  was  predominant  and 
the  P/R  ratio  was  greater  than  1.0. 

A  summary  of  diurnal  metabolism  measurements  taken 
in  the  bay  in  the  summer  of  1973  is  given  in  Table  4.   Re- 
sults are  grouped  in  the  table  according  to  major  features 
characterizing  the  portion  of  the  bay  from  which  the  samples 
were  taken  (Odum  e_t  aj^.  ,  1974a).   During  the  sampling  pe- 
riod, metabolism  (daytime  net  photosynthesis  plus  night 
respiration)  ranged  from  3.1  g  O^/m  /24  hrs  at  Station  L 
to  21.6  g  Q^/m^/2^    hrs  at  the  FSUML  boat  basin.   The  latter 
station  was  the  deepest  station  where  a  measurement  was 
taken  and  this  in  part  accounts  for  the  large  metabolism, 
most  of  which  was  respiration.   There  was  some  variability 
at  stations  where  metabolism  was  measured  twice. 

Figure  11  shows  total  metabolism  plotted  against 

time.   No  clear  seasonal  trend  was  evident  over  the  four- 

2 

month  sampling  period.   Average  metabolism  was  8.5  g  02/m  / 

24  hrs.   Figure  12  shows  daytime  net  photosynthesis  plotted 
against  nighttime  respiration  for  each  station.   The  diago- 
nal line  represents  a  P/R  ratio  of  1.0.   When  plotted  in 
this  fashion  those  stations  in  East  Bay  (K,  L,  M,  ),  near 
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Fig.  n.   Graph  of  bay  productivity  for  June-September,  1973. 
Each  point  represents  sum  of  daytime  net  photosyn- 
thesis and  nighttime  respiration  from  a  24-hour  reC' 
ord  of  oxygen. 


65 


20  _ 


15    _ 


CVJ 

o 


10 


5    _ 


J  '  A 

MONTH 


Fig.  12.   Graph  of  daytime  net  photosynthesis  as  a  function 
of  nighttime  respiration  for  systems  in  Apalachi- 
cola  Bay  (grass  flats  x  ;  boat  basins  A  ;  medium 
salinity  plankton  system  •  and  oligohaline  system 

o  ).   The  diagonal  line  represents  a  P/R  ratio 
of  one. 
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the  river  mouth  (Z,  J)  and  stations  suspected  of  receiving 
man-related  wastes  (I,  r',  FSUML  boat  basin)  all  had  P/R 
ratios  less  than  1.0,   The  medium  salinity  plankton  system 
stations  (A,  AA,  B,  C,  D,  E,  F,  H,  R,  S,  T,  and  U)  had  P/R 

ratios  of  about  1.0  or  somewhat  greater.   Average  metabolism 

2 
in  this  area  was  6.8  g  O^/m  /24  hrs.   The  grass  flat  stations 

(N,  and  FSUML  stations)  all  had  P/R  ratios  greater  than  1.0. 

Average  metabolism  at  the  grass  flat  station  was  9.6  g  0^/ 

m^/24  hrs. 

Organic  Carbon  Measurements 

Fourteen  stations  were  sampled  for  particulate  and 
dissolved  organic  carbon  between  September  20  and  September 
30,  1973.   Organic  carbon  data  are  summarized  in  Table  5. 

The  average  concentration  of  total  and  particulate  organic 

3 
carbon  in  the  bay  area  was  5.9  and  0.8  g  C/m  ,  respectively. 

Surface  and  bottom  concentrations  of  total  and  particulate 

3 
organic  carbon  in  the  river  were  4.5  and  0.6  g  C/m  and 

3 
3.3  and  0.62  g  C/m  ,  respectively.   The  dissolved  to  particu- 
late organic  carbon  ratio  for  an  average  of  all  bay  stations 
was  6.4  This  ratio  was  5.3  for  the  river.   In  general, 
organic  matter  concentrations  in  the  bay  were  higher  than 
thoseintheriver. 

Data  Assembled  from  Other  Sources 


Given  in  Fig.  13  is  a  plot  of  monthly  average  river 
flow  from  October  1928  through  October,  1972  (Mann,  1975). 
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Table  5 


Summary  of  organic  carbon  measurements  taken  in  Apala- 
chicola  Bay  in  September,  1973.   Station  locations  are 
given  in  Fig.  4.   Numbers  in  parentheses  indicate  the 
number  of  replicates  averaged. 


Station 


Date    Total  organic   Particulate 
carbon,        organic 
mg  C/1     carbon,  mg  C/1 


1973 


FSUML' 


(boat  basin) 


Sept.  20     8.5  (4) 


1.3  (2) 


FSUML 
(1/4  mi .  offshore)   Sept.  20 


Station  I 


John  Gorrie  Bridge 
(western  end) 


River  Mouth 
Station  C 
Station  AA 
Station  D 


Sept.  20 


John  Gorrie  Bridge    Sept.  20 
(eastern  end) 


Sept.  20 

Sept.  30 

Sept.  30 

Sept.  30 


8.6  (3) 
7.8  (1) 
4.8  (4) 


Sept.  20     5.2  (4) 


6.1  (2) 
4.7  (3) 
5.4  (3) 
5.4  (2) 


1.2  (2) 
0.7  (3) 

0.6  (3) 

0.4  (2) 
0.7  (2) 
0.9  (1) 
0.6  (1) 


Florida  State  University  Marine  Lab  at  Turkey  Point, 
Franklin  County,  Florida. 

The  John  Gorrie  Bridge  crosses  Apalachicola  Bay  con- 
necting the  towns  of  East  Point  and  Apalachicola. 
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Table  5  -  continued 


Station 


p.  .    Total  organic   Particulate 
,Q7^      carbon,       organic 
'   "^      mg  C/1     carbon,  mg  C/1 


Station  F 

Station  H 

Station  L 

River 
(bottom) 


River 
(surface) 


Sept.  30  4.2  (1) 

Sept.  30  5.0  (1) 

Sept.  30  5.5  (3) 

Sept.  30  4.5  (4) 

Sept.  30  3.3  (3) 


0.6  (1 

0.7  (1 

0.7  (1 

0.6  (1 

0.6  (1 


The  river  station  was  taken  about  a  half  mile 
downstream  of  the  railroad  bridge  crossing  the  Apalachi- 
c  0 1  a  River. 


Fig.  13.   Monthly  average  flow  (cfs)  of  the  Apalachicola 
River  at  Chattahoochee,  Florida  from  October, 
1928  to  October,  1972. 
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Yearly  average  flow  for  this  period  is  shown  along  the 
dotted  line.   During  this  period,  average  flow  was  21,690 
cubic  feet  per  second  (cfs).   The  maximum  and  minimum  flows 
were  of  293,000  and  4,950  cfs  recorded  in  March,  1929  and 
October,  1954,  respectively.   The  flow  rates  presented  in 
Fig.  13  were  recorded  at  Chattahoochee,  Florida.   An  esti- 
mate of  river  flow  into  Apalachicola  Bay  was  made  by  adding 
the  flow  of  the  Chipola  River  to  that  recorded  at  Chatta- 
hoochee and  increasing  the  total  by  10  percent  (Dawson, 
1955a).  With  this  adjustment,  the  44-year  average  flow  be- 
comes 25,179  cfs. 

Figure  14a   shows  average  monthly  river  flow  over 
the  period  1929  through  1972.   Flow  was  adjusted  to  approxi- 
mate the  volume  which  enters  the  bay  as  described  above. 
Peak  flow  occurred  in  March  (52,107  cfs)  and  minimum  flow 
in  September  (12,801  cfs).   Peak  flow  occurred  as  a  pulse 
from  January  through  April  with  the  remaining  monthly  aver- 
age flows  near  or  below  the  yearly  average.   Figure  14b 
shows  monthly  average  river  flow  for  15-year  periods  before 
and  after  construction  of  Jim  Woodruff  Dam  at  Chattahoochee, 
Florida.   The  combined  average  flow  for  the  full  43-year 
recording  period  is  also  shown.   All  three  curves  are  sim- 
ilar.  Apparently,  the  relatively  small  storage  capacity  of 
Lake  Seminole  behind  Jim  Woodruff  Dam  did  not  affect  the 
timing  of  discharge. 


Fig.  14.   Monthly  river  flow,  (a)  averaged  for  the  period 
1929  through  1972;  and  (b)  averaged  for  the  15- 
year  period  before  and  after  construction  of  Jim 
Woodruff  Dam  at  Chattahoochee,  Florida. 
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Given  in  Figs.  15-17  are  estuarine  data  developed 
by  Estabrook  (1973)  in  a  study  of  phytopl ankton  ecology 
and  hydrography  of  Apalachicola  Bay.   They  were  redrafted 
and  presented  here  for  use  in  both  the  diurnal  activity  and 
the  regional  models. 

Figure  15a  shows  monthly  average  insolation  for 
1972-1973.   The  monthly  averages  show  a  seasonal  cycle 
typical  of  Florida  with  maximum  insolation  in  May  and 
lower  values  later  in  the  summer  due  to  rainy  season  cloud 
cover.   Minimum  insolation  occurred  in  December.   Figure  15b 
shows  average  surface  and  bottom  salinities  for  the  bay. 
The  average  salinity  of  the  bay  was  lowest  in  February  dur- 
ing high  river  flow  and  highest  in  September  during  low 
river  flow.   The  average  difference  between  surface  and 
bottom  salinities  was  8.5  ppt  with  a  seasonal  range  of  1 
to  18  ppt.   Detailed  salinity  maps  are  given  in   Figs.  8 
and  9  and  in  Estabrook  (1973). 

Average  water  temperature  is  shown  in  Fig.  16a. 
The  maximum  average  temperature  was  about  60°C  in  September 
with  a  minimum  of  about  10°C  in  February.   Figure  16b  shows 
average  values  for  phytopl ankton  production  as  measured 
by  Estabrook  (1973)  using  the  C-14  technique.   The  dashed 
lines  are  95%  confidence  limits.   Mean  daily  phytopl ankton 

production  ranged  between  1.69  g  C/m  /day  in  April  to  0.06 

2 
g  C/m  /day  in  February.   However,  there  was  a  wide  range  in 

productivity  in  the  bay  at  any  one  sampling  time  as  indi- 
cated by  the  large  confidence  intervals. 


Fig,  15.   Monthly  average  insolation  (a)  and  average  surface 
and  bottom  salinity  for  Apalachicola  Bay  (b).   Data 
are  from  Estabrook  (1973). 
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Fig.  16.   Average  seasonal  water  temperature  (a)  and  average 
daily  phytoplankton  production  in  Apalachicola  Bay 
(b).   Dashed  lines  in  (b)  are  95%  confidence  limits. 
Data  are  from  Estabrook  (1973). 
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Shown  in  Figs.  17a  and  17b  are  seasonal  concentra- 
tions of  nitrate  and  ammonia  in  the  surface  and  bottom 
waters  of  the  bay  and  from  the  surface  of  the  river  (Esta- 
brook,  1973).   Nitrate  concentrations  in  the  river  were 
rather  constant  ranging  from  0.18-0.20  mg  N/1.   Average 
nitrate  concentrations  in  the  bay  showed  a  seasonal  cycle 
with  high  concentrations  in  February  of  about  0.18  mg  N/1 
when  river  flow  was  high.   Low  concentrations  of  about  0.05 
mg  N/1  occurred  in  the  fall  when  river  flow  was  low.   Bottom 
nitrate  concentrations  were  generally  less  than  surface 
concentrations  except  when  river  flow  was  low,  at  which 
time  there  was  little  difference  between  surface  and  bottom. 

Ammonia  concentrations  are  shown  in  Fig.  17b.   Con- 
centrations in  the  bay  ranged  from  7-40  yg  N/1.   No  clear 
seasonal  trend  was  evident.   Bottom  water  concentrations 
of  ammonia  were  generally  higher  than  surface  water  levels. 
Ammonia  concentrations  in  the  river  varied  considerably 
ranging  from  about  7-35  yg  N/1  over  the  year.   At  times 
ammonia  concentrations  in  the  bay  were  higher  than  those 
in  the  river. 

Shown  in  Fig.  18  are  monthly  estimates  of  fish 
and  benthic  invertebrate  biomass  from  15  stations  in  Apa- 
lachicola  Bay  (Menzel  and  Cake,  1969).   Oyster  reef  bio- 
mass was  not  included  in  these  estimates.   Biomass  was  in 

2 
g/m  wet  weight.   While  there  was  considerable  variation 

in  biomass  between  samples  from  any  one  sampling  data. 


Fig.  17.   Seasonal  record  of  nitrogen  concentrations  in  the 

Apalachicola  River  and  in  surface  and  bottom  waters 
of  Apalachicola  Bay.   Data  are  from  Estabrook  (1973). 
(a)  Nitrate,  (b)  Ammonia. 
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Fig.  18.  Seasonal  estimates  of  combined  fish  and  benthic 
invertebrate  biomass  in  Apalachicola  Bay.  Data 
were  adapted  from  Menzel  and  Cake  (1969). 
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biomass  was  larger  in  the  warmer  periods  of  the  year.  Aver- 

2  -    ? 

age  biomass  ranged  from  0.12  g/m   in  February  to  b  g/m 

in  August. 

Since  no  direct  estimates  of  oyster  biomass  for 
Apalachicola  Bay  were  available,  several  indirect  methods 
of  calculating  biomass  were  used.   These  are  given  in  Table 
6.   In  the  first  calculation  biomass  was  estimated  from 
1970  oyster  harvest  data  with  corrections  for  predation  mor- 
tality and  harvest  losses.   Total  weight  in  pounds  was  then 
converted  to  grams  dry  weight.   Corrections  for  predation 
and  harvest  losses  were  estimated  by  Menzel  (1973). 

In  the  second  method  harvest  per  acre  was  multi- 
plied by  total  area  of  oyster  reef  and  then  corrected  for 
predation  and  harvest  losses.   This  calculation  was  re- 
peated for  both  marginal  and  very   productive  oyster  reefs. 
In  the  third  method,  areas  of  three  grades  of  oyster  bars 
were  planimetered  from  an  oyster  reef  survey  map  (U.S.  Army 
Corps  of  Engineers,  1968).   From  data  developed  by  Menzel 
et^  aj[.  (1966)  concerning  oyster  size  and  numbers  occurring 
on  reefs,  total  weights  of  both  adult  (defined  as  being  of 
harvestable  size)  and  juvenile  (all  oysters  less  than  har- 
vestable  size)  oysters  were  calculated.   Total  weights 

Q 

varied  from  3.1  x  10   g  in  the  first  calculation  to  29.0 

o 

X  10   g  in  (b)  of  calculation  II.   Calculation  III  yielded 

o 

an  intermediate  estimate  of  14.2  x  10  g.   The  estimate 
from  calculation  III  was  used  in  simulation  models  both 
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because  the  result  was  intermediate  in  the  range  of  values 
obtained  and  because  it  gave  an  estimate  of  both  adult  and 
juvenile  biomass.   The  estimate  obtained  in  11(b)  (715  bu/ 
acres)  is  a  ^ery   high  yield  and  probably  not  representative 
of  average  conditions.   Calculations  I  and  Ila  and  lib  in- 
clude only  adult  biomass  and  thus  may  underestimate  total 
oyster  biomass.   Adult  biomass  estimates  Ila  and  III  are 
similar  (9.0  x  10^  gms  and  8.7  x  10^  gms )  . 

Model  of  Diurnal  Activities 


The  diurnal  activity  model  shown  in  Fig.  19  has 
component  storages  of  phytopl ankton  biomass  (Q,),  consumers 
(Q2),  dissolved  oxygen  (Q3),  detritus  (Q^),  and  nitrogen 
(Qg).   Tracing  flows  from  left  to  right  across  the  diagram 
inputs  from  external  sources  of  sunlight  (J  ),  nitrogen 
(N),  detritus  (D),  and  advected  oxygen  (A)  enter  storage 
modules.   At  the  top  of  the  diagram,  oxygen  in  the  atmo- 
sphere (As)  exchanges  with  oxygen  dissolved  in  the  water 
depending  on  the  saturation  gradient  (As-Q^)  and  the  amount 
of  turbulence  (T).   On  the  left,  sunlight  and  nutrients 
interact  mul ti pi i cati vely  to  produce  oxygen  (Jp)  and  phyto- 
plankton  biomass  (J-i).   Losses  from  the  phytoplankton 
storage  (Q^)  include  respiration  (J4),  death  (J^y).  and 
exports  (^13)-   Inputs  to  the  nitrogen  storage  (Q^)  come 
from  an  external  source  (N),  and  from  recycle  along  path- 
ways ^12'  '^]S'    ^^^   '^16'       Losses  include  flushing  (v]^^) 
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and  utilization  in  photosynthesis  (J3).   Inputs  to  the 
detritus  storage  (Q,)  come  from  phytopl ankton  death  (Jny) 
and  an  external  source  (D).   Losses  from  this  storage  in- 
clude microbial  and  consumer  utilization  (J,  and  Jo)  and 
flushing  (J-jg).   Microbe  biomass  is  part  of  the  detritus. 
The  oxygen  content  in  the  water  (Q  )  is  determined  by  the 
net  balance  of  respiratory  utilization  (Jg,  J,,  Jg),  phyto- 
plankton  production  of  oxygen  (J2),  advection  of  oxygen 
into  (Ja)  and  out  of  the  system  {<^ia)>    and  by  reaeration 
(J12).   All  inputs  to  the  model  are  constant  flow  sources 
except  sunlight  which  enters  as  a  chopped  si  ne  wave  or  square 
wave  to  simulate  daily  sunlight  conditions.   Temperature 
and  salinity  are  held  constant  over  24  hours   because  only 
small  changes  were  observed  in  the  field.   Justifications 
for  the  model  configuration  are  given  in  Fig.  20.   Calcu- 
lations and  assumptions  used  in  evaluating  the  model  are 
given  in  Table  C-1.  Scaled  equations  are  given  in  Table 
C-2.   The  analog  patching  diagram  is  given  in  Fig.  C-1. 

Simulation  Results  and  Comparisons  with  Field  Data 


Simulation  results  obtained  from  the  diurnal  acti 
vity  model  are  given  in  Figs.  21  and  23-27.   Comparisons 
of  "control  condition"  model  output  with  field  data  are 
given  in  Fig.  22.   Table  7  is  a  tabular  summary  of  model 
output  for  al 1  runs  . 


Fig.  20.   Justification  of  the  configuration  used  in  the 
diurnal  activity  model  shown  in  Fig.  19a.   (a) 
phytopl ankton;  (b)  consumers;  (c)  oxygen;  (d) 
detritus;  (e)  nitrogen. 


Fig.  20a.   Phytopl ankton  were  chosen  as  a  state  variable  be- 
cause they  were  the  major  biological  producers  of 
oxygen  and  organic  matter  in  the  bay,  and  because 
the  stock  of  phytopl ankton  were  subject  to  frequent 
changes  (Estabrook,  1973).   Photosynthesis  was  mod- 
eled with  a  multiplicative  workgate  because  the  re- 
sponse is  of  the  limiting  factor  type  observed  for 
photosynthesis.   In  this  model  either  nitrogen  or 
sunlight  could  be  the  limiting  factor.   The  special 
equation  for  Or  is  given  in  Appendix  C-2.   Phyto- 
plankton  respiration  was  assumed  to  be  a  multipli- 
cative function  of  biomass  and  oxygen.   Nutrient 
recycle  was  assumed  to  be  a  constant  proportion  of 
the  respired  material.   Phytoplankton  death  and 
excretions  were  accounted  for  by  linear  drain  J, y. 
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Fig,  20c.   Oxygen  was  chosen  as  a  state  variable  because 
its  diurnal  behavior  was  a  main  interest  in 
this  model.   Photosyntheti c  production  of  oxy- 
gen (J2)  and  multiplicative  losses  (J5,  Js.  and 
Jg)  have  already  been  explained.   The  reaera- 
tion  pathway  (0-13)  is  the  same  as  the  standard 
turbulent  diffusion  equation.   Oxygen  advection 
(J/\)  was  assumed  to  be  constant  over  a  diurnal 
cycle.   Linear  pathway  Ji4  accounts  for  a  diur- 
nal cycle.   Linear  pathway  J-]4  accounts  for  oxy- 
gen exports  via  flushing.   Pathways  J2.  J5>  Js> 
and  Jg  are  nonlinear.   See  Fig.  19  for  details. 


99 


OXYGEN     IN     ATMOSPHERE 


DISSOLVED 
OXYGEN 


MIXING    ENERGY 


^EXPORTS 


FROM    PHYTOPLANKTON 

TO    PHYTOPLANKTON 


TO   CONSUMERS 


TO   DETRITUS 


Fig.  20d.   Detritus  was  chosen  as  a  state  variable  because 
it  was  the  major  storage  of  organic  matter  in 
the  water  column,  a  major  source  of  bound  nitro- 
gen, and  had  several  interactions  with  oxygen. 
The  hexagon  group  symbol  was  used  to  recognize 
that  detritus  in  this  model  included  some  living 
components  (microbes).   The  group  symbol  does  not 
imply  additional  pathways  beyond  those  shown. 
Detritus  input  (Jq)  from  an  external  source  (D) 
was  assumed  to  be  constant.   Input  due  to  phy- 
toplankton  death  and  excretion  was  assumed  to  be 
in  proportion  only  to  the  phytopl ankton  stock. 
Detritus  export  was  assumed  to  be  in  proportion 
to  the  detritus  stock.   The  respiration  pathway 
was  dependent  on  both  oxygen  and  detritus  which 
are  necessary  for  any  aerobic  respiratory  proc- 
ess . 
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Fig.  20e.   Nitrogen  was  chosen  as  a  state  variable  for  sev- 
eral reasons.   First,  nitrogen  is  required  in 
photosynthesis  and  in  turbid  coastal  plain  es- 
tuaries of  the  United  States  may  be  the  primary 
limiting  nutrient  (Pomeroy  ejt  aj_.  ,  1972).   Sec- 
ondly, the  nitrogen  stock  could  be  expected  to 
change  due  to  changes  in  recycle  rates  or  due  to 
changes  in  the  input  rate.   Nitrogen  input  from 
an  external  source  was  assumed  to  be  constant 
during  each  simulation  run.   Estabrook  (1973) 
found  fairly  constant  concentrations  of  total 
inorganic  nitrogen  in  the  river  (Fig.  17). 
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Computer  simulations  for  the  control  condition  are 

shown  in  Fig.  21.   Sunl-'ght  is  shown  on  the  top  of  Fig.  21. 

2 
Maximum  daily  insolation  was  estimated  to  be  700  Kcal/m  /hr 

at  1200  hrs  (Estabrook,  1973).   Day  and  night  were  divided 

into  equal  12-hour  segments  which  closely  approximated  the 

day  length  for  late  September.   Consumer  and  detritus  stocks 

remained  at  initial  condition  levels  through  the  4.5  days 

which  the  model  was  run.   However,  phytoplankton  biomass 

2 
showed  a  diurnal  change  ranging  from  0.8  g/m  at  dawn  to 

2 
2.9  g/m  in  the  late  afternoon.   Indications  of  diurnal 

changes  in  phytoplankton  stocks  have  been  reported  (Piatt 
and  Conover,  1971).   Dissolved  oxygen  ranged  between  15.4 
and  19.9  g/m   (water  column  of  2.1  meters)  with  low  and 
high  values  at  0700  and  1600  hours,  respectively.   Phyto- 
plankton production  was  the  major  source  producing  the 
daily  pulse.   Advection  in  this  run  was  small.   The  sharp 
increase  in  oxygen  in  the  morning  was  due  to  high  phyto- 
plankton activity  under  conditions  when  nitrogen  was  high 
and  sunlight  increasing  while  respiratory  activity  was  low 
due  to  depressed  phytoplankton  stocks.   Diffusion  was  also 
small  during  the  morning.   The  afternoon  and  evening  decline 
in  oxygen  was  not  as  sharp  as  the  morning  increase.   Nitro- 
gen also  exhibited  a  diurnal  change  ranging  between  a  high 

2 
of  0.20  at  dawn  to  0.04  g  N/m  at  sunset.   The  daytime  de- 
cline in  nitrogen  was  sharper  than  the  nighttime  recovery. 
Gross  photosynthesis  reached  a  maximum  rate  of  0.95  g 


Fig.  21.   Graphs  of  diurnal  activity  simulated  with  the  model 
Fig.  19a,  using  transfer  coefficients  calculated 
from  data  in  Fig.  19b. 
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02/ni  /hr  at  1100  hours.   The  afternoon  decline  in  gross 
photosynthesis  was  more  gradual  than  the  early  morning 
increase.   This  was  due  both  to  declining  sunlight  and  a 
sharply  declining  nitrogen  stock. 

Figure  22  compares  control  condition  model  output 
to  field  data  for  which  diurnal  measurements  were  avail- 
able.  Neither  detritus  nor  consumers  changed  appreciably 
over  the  four-day  simulation  run  retaining  values  close  to 
those  measured  or  estimated  (see  Tables  5  and  6  and  Fig.  18). 
In  general,  seasonal  changes  were  larger  than  diurnal  changes 
for  these  storages.   Nitrogen  stocks  (NH^  +  NO^)  changed  by 

a  factor  of  5  over  a  24-hr  period  ranging  from  0.04  to  0.20 

2 
g  N/m  .   No  field  data  are  available  to  estimate  the  actual 

diurnal  change  of  available,  labile  nitrogen.   In  a  model 

developed  prior  to  the  one  reported  here,  the  nitrogen 

storage  included  labile  and  less  labile  fractions  as  total 

nitrogen  (NO^  +  NH,  +  particulate  and  dissolved  organic 

nitrogen)  and  diurnal  changes  were  much  smaller.   In  that 

model,  conversion  of  organic  nitrogen  compounds  to  NH„  was 

made  a  function  of  both  sunlight  and  the  total  nitrogen 

stock. 

Figures  23-27  show  simulation  output  of  model 

variables  under  a  range  of  conditions  which  may  occur  in 

the  study  area.   These  included  reduction  in  sunlight  due 

to  increased  turbidity  (Fig.  23),  increased  nitrogen  (Fig. 

24),  and  detritus  (Fig.  25)  input  rates,  changes  in  microbe 


Fig.  22.   Comparisons  of  field  data  with  computer  output  from 
the  model  given  in  Fig.  19a.   Computer  output  was 
obtained  using  transfer  coefficents  calculated  from 
data  used  to  evaluate  the  model.   (a)  oxygen,  (b) 
percent  saturation,  (c)  gross  photosynthesis. 
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respiration  under  nutrient  and  detritus  loading  (Fig.  26), 
and  finally  a  mixture  of  these  changes  under  reduced  rates 
of  oxygen  exchange  between  the  atmosphere  and  water  (Fig. 
27).   The  results  are  summarized  in  Table  7. 

In  Fig.  23  sunlight  was  changed  to  a  square  wave 

2 
with  an  average  input  of  700  Kcal/m  /hr.   The  change  caused 

only  small  changes  in  variable  responses  most  notably  in 

gross  photosynthesis  where  a  sharp  morning  peak  (1.33  g 

2 
O^/m   /hr)  developed.   This  peak  was  followed  by  a  sharp 

2 
decline  to  lower  values  (0.46  g  02/m  /hr)  for  the  after- 
noon.  The  decline  was  caused  by  the  decreasing  nitrogen 
stock.   Other  variables  were  comparable  to  control  condition 
levels.   The  slight  increase  in  most  variables  was  a  result 
of  more  sunlight  entering  the  system  due  to  the  difference 
in  area  beneath  a  square  wave  and  sine  wave  of  the  same  am- 
plitude.  When  sunlight  was  decreased  to  0.5  of  the  initial 
value,  gross  primary  production  did  not  show  a  sharp  morning 
peak  and  the  nitrogen  limitation  was  not  as  sharp.   Nitrogen 
stocks  increased  slightly  but  other  variables  remained  about 
the  same.   At  0.25  of  the  initial  sunlight  input  rate  the 
diurnal  range  in  phytopl ankton  and  oxygen  values  both  de- 
creased but  minimum  values  were  similar  to  the  control 

situation.   Gross  primary  production  was  constant  at  0.33 

2 
g  Op/m  /hr  and  did  not  show  effects  of  nitrogen  limitation. 

The  diurnal  peak  in  the  nitrogen  stock  was  80%  greater  than 

in  the  control  condition.   Major  changes  in  model  output 


Fig.  23.   Simulation  results  from  the  model  given  in  Fig.  19a 
with  input  of  sunlight  as  a  square  wave.   Families 
of  curves  are  shown  for  each  variable  with  desig- 
nations of  1 ,  2,  and  3  corresponding  to  sunlight  in- 
puts of  normal,  1/2  normal,  and  1/4  normal. 
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under  conditions  of  decreasing  available  sunlight  were 
mainly  decreased  gross  photosynthesis,  and  increased  stocks 
of  nitrogen. 

Diurnal  responses  to  increased  nitrogen  input  rates 
are  given  in  Fig.  24.   Under  these  conditions  consumers  and 
detritus  increased  only  slightly  above  control  values. 
Gross  primary  production  increased  by  a  factor  of  1.6  be- 
yond the  control  value  for  both  treatment  cases.   The  nitro- 
gen stock  also  increased  but  the  magnitude  of  the  diurnal 
change  decreased  from  a  factor  of  5  in  the  control  run  to 
a  factor  of  3  in  the  run  with  a  5x  input  rate.   Oxygen  in- 
creased both  in  amount  and  in  the  amount  of  change  over  a 
daily  cycle.   In  the  cases  with  decreased  sunlight  (Fig.  23) 
and  increased  nitrogen  input  (Fig.  24),  the  amount  of  oxygen 
in  the  water  did  not  fall  below  control  values  at  any  time. 

Fig.  25  shows  model  responses  when  the  detritus 
input  rate  was  increased  by  factors  of  2,  5,  and  10  times 
the  normal  input  rate.   The  detritus  stock  increased  by 
factors  of  about  1.6,  3,  and  6  over  the  control  condition. 
All  three  input  rates  had  the  effect  of  lowering  oxygen 
levels  below  control  condition  levels  and  increasing  the 
range  of  diurnal  oxygen  changes.   In  the  control  condition, 
oxygen  ranged  from  16.4-14.2  (AO2  =  2.8  ppm)  whereas  in 
the  lOx  input  rate  case  oxygen  ranged  from  11.9-16.5 
(AOp  =  4.6  ppm).   Gross  photosynthesis  increased  by  factors 
of  1.1,  1.5,  and  1.6  above  the  peak  control  rate  of  0.90g 


Fig.  24.   Simulation  results  from  the  model  given  in  Fig. 
19a  with  varying  rates  of  nitrogen  input.   Fami- 
lies of  curves  are  shown  for  each  variable  with 
designations  of  1,  2,  and  3  corresponding  to  in- 
puts of  normal,  twice  normal,  and  five  times  the 
normal  rate. 
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Fig.  25.   Simulation  results  from  the  model  given  in  Fig. 

19a  with  varying  rates  of  detritus  input.   (1)  con- 
trol condtions;  (2)  2x  normal  input  rate;  (3) 
5x  normal  input  rate;  (4)  lOx  normal  input  rate, 
(a)  detritus,  oxygen,  percent  saturation,  photo- 
synthesis and  phytopl ankton  ;   (b)  nitrogen  and 
consumer  stock. 


117 


liib- 

DETRITUS 

ru— -i 

'^g  100- 

"^ 

^^^^t-""^^"'^ 

en 

^^.^^■■'^^^ 

75- 

^ - 

t   50- 

■ — ■ — 

h- 

Ul 

r^ 

Q    25- 

— 

V| 

0- 

1 

1—           r          -I — ' 

30- 

OXYGEN 

^  20- 

*""^S5»s.,_  ^»{^^^!^=^ 

/--=w^ 

^ 

e^ 

0XY6EN, 

O 

1 

^**^^^!^        ^^^ 

C= 

r~ 

0^ 

1 

i                                      1 

1 — 

rl50-j  PERCENT     SATURATION 


o 

H- 

Z  W 
<    G 


Oi 


6- 

PHYTOPLANKTON 

4- 

^s. 

/5C^4     / 

S^  ^ 

2- 

"X,...J^^W^ 

r*^^^-^ 

"<^J^ 

0- 

1                 

"  1                  1 

^^^ 

TIME,    DAYS 
(a) 


118 


NITROGEN 


119 


02/m  /hr.   The  increase  in  photosynthesis  observed  under 
increased  detritus  input  rates  was  not  so  large  as  the  in- 
crease observed  for  similar  increases  in  nitrogen  input 
rates.   Consumer  biomass  in  the  model  increased  by  factors 
of  1.1,  1.4,  and  1.7  beyond  the  control  case. 

Figure  26  shows  model  output  under  conditions  where 
microbial  respiration  of  detritus  (Jy)  was  increased  by  a 
factor  of  2.0.   Microbial  respiration  was  not  known  and 
was  increased  to  check  the  sensitivity  of  the  model  to  this 
type  of  change.   Microbial  use  of  oxygen  (J^)  and  nitrogen 
recycle  (J-ic)  were  also  increased  by  factors  of  2.0.   In 
Fig.  26a  nitrogen  input  was  at  2x  the  control  rate  and  in 
Fig.  26b  detritus  input  was  at  2x  the  control  rate.   In  gen- 
eral, only  small  changes  were  observed  when  these  curves 
were  compared  to  curves  in  Fig.  24  and  25  corresponding 
to  2x  input  rate  increases  in  nitrogen  and  detritus,  re- 
spectively.  In  both  cases,  minimum  oxygen  levels  were 
lower  than  observed  under  the  same  increased  input  condi- 
tions with  no  increased  microbial  activity.   Little  change 
was  observed  in  maximum  values. 

Figure  27  shows  simulation  graphs  for  oxygen.   In 
all  three  graphs  the  rate  of  oxygen  exchange  between  the 
atmosphere  and  water  was  decreased  to  0.1  of  the  control 
condition  value.   In  all  cases  oxygen  levels  dropped  below 
control  conditions.   In  case  (3)  maximum  and  minimum  levels 
decreased  by  10.8  ppm  and  9.6  ppm  respectively,  below 


Fig.  26.   Results  of  simulating  the  model  given  in  Fig.  19a 
with  microbe  respiration  increased  by  a  factor  of 
two  above  control  values.   In  (a)  nitrogen  input 
was  twice  the  control  rate.   In  (b)  detritus  input 
was  twice  the  control  rate. 
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Fig.  27.   Graphs  of  oxygen  resulting  from  simulating  the  model 
in  Fig.  19a  with  the  reaeration  rate  reduced  to  10 
percent  of  control  conditions,  varying  inputs  of 
detritus,  nitrogen,  and  sunlight,  and  with  sunlight 
as  a  square  wave. 

1.  control  conditions 

2.  detritus  input  4x  normal  rate 

3.  detritus  input  lOx  normal  rate 

4.  detritus  and  nitrogen  inputs  4x  normal  rate 

5.  control  conditions  but  with  sunlight  at  1/2 
normal  input  rate 

6.  sunlight  input  1/2  normal  rate;  detritus  and 
nitrogen  input  rates  at  4x  normal  rate. 
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control  levels.   The  oxygen  range  also  increased  in  cases 
(2)  and  (3)  over  the  control  case  (1).   Oxygen  in  case  (4) 
had  a  diurnal  range  of  7.7  ppm.   In  (c)  both  plots  show 
oxygen  levels  below  control  values.   In  case  (5)  the  diurnal 
range  of  oxygen  was  reduced  to  4.3  ppm.   In  case  (6),  both 
maximum  and  minimum  oxygen  values  were  lower  than  in  any 
other  case.   The  diurnal  range  of  oxygen  values  increased 
over  case  (5)  and  was  similar  to  the  control  condition  range, 

Characteristics  of  Franklin  County  Region 

Results  of  characterizing  the  regional  system  of 
Franklin  County  were  organized  with  the  regional  model  in 
Fig.  28.   Data  were  assembled  for  evaluating  the  main  fea- 
tures of  the  model  in  Table  8  and  Figs.  30  to  34.   Included 
are  graphs  for  fuel  use,  tourism,  population,  income, 
sales,  and  oyster  and  other  fishery  yields.   Included 
with  the  model  is  a  description,  justification  of  the 
model  configuration  (Fig.  29),  simulation  graphs  (Figs. 
35  to  49),  and  comparisons  of  simulation  output  with  field 
data . 

Regional  Model 


Shown  in  Fig.  28  is  the  detailed  model  of  Franklin 
County  and  Apalachicola  Bay,  Florida.  This  is  an  elabora- 
tion of  the  energy  circuit  language  diagram  given  in  Fig. 
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(b) 


Q.= 


^46^5^11%   "    ^44^5''3    "    ^37^4    "    ^9^6    "    ^9617^5^2 


-    S3'^7'Ql^4''l 


=    ^30^8^3Q5    '    hl%    '    S7^]2Q3'^5    "    ^33^6    "    ^4^6    "    ^35^6 
=    h2h^3   '    ^13^6   "    ^4^9    '    "^IsQlO    "    ^7^7    "    ^S^Io    " 


^02^7^9 
■^2^    ■    '^11^! 


•10 


'n 


^^4^1    ^    h'^^o    '    h%    -    ^01^9^7    -    ^6^9    (^-^8^ 
k25Q9{l-Q8)    -    k^gQio    -    ^6^10    -    1^7^10^12 
^9^10    -    ^8^11^12    -    40^11^5 


^12    =    ^7^10^12    '    ^8'^11^12    -    ^^12    -    ^0^12^8   '    Kho 


'^3 


"^21^1^2    "    ^22^6    "    "^23^13 
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3  to  summarize  structure,  function,  inputs,  exchange,  and 
interactions  in  the  county.   The  model  is  given  in  equation 
form  in  Fig.  28b.   Each  transfer  coefficient  (k)  has  a  sub- 
script that  corresponds  directly  to  the  subscripted  flows 
(J's)  given  on  the  energy  circuit  diagram.   The  pathways 
in  Fig.  28  can  thus  be  traced  to  the  differential  equations, 
Values  used  in  evaluating  the  model  are  given  in  Table  8 
with  all  calculations,  assumptions,  and  references  given 
in  Appendix  Table  D-1.   The  scaled  differential  equations 
which  describe  model  behavior  are  given  in  Appendix  D-2. 
The  analog  patching  diagram  is  given  in  Appendix  Figs.  D-1 
and  D-2.   Parts  of  the  model  are  isolated  in  Fig.  29  with 
more  comments  on  the  reasons  for  the  configuration  used. 

The  regional  model  is  composed  of  three  major  group- 
ings:  the  urban  system,  water  and  plankton  system,  and 
the  oyster  reef  system.   Land-use  distri buti on  i s  shown  by 
two  storages  located  to  the  upper  left  of  the  urban  sector 
of  the  model  (Q^  and  Q^) .      Urban  sources  include  tourists 
(1^2)  ^"d  tourist  money  (1^3).  new  residents  (Ig)  and  new 
resident  money  (Ig),  investment  (I^)  and  government  money 
(Ig)  and  goods,  fuels  and  services  (I-,).   Urban  storages 
include  image  (Q3),  capital  (Q^),  structure  (Q^)  and  popu- 
lation ( Q  g ) . 

In  the  urban  sector,  image  had  inputs  from  all 
major  productive  processes  of  the  region  including  one 
from  the  input  to  structure  i^^-^'),    development  of  new 
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Table  8 

Numbers  used  to  evaluate  regional  model  of  Franklin  County 
and  Apalachicola  Bay  given  in  Fig.  28a. 


Designation  on 
model  (Fig.  28a) 
and  in  Appendix 
Table  D-1 


Storages 


•10 


'11 


'12 


'13 


Name 


Natural  lands 

Developed  lands 

Image 

Capital 

Structure 

Population 

Nitrogen 

Freshwater  and 
turbidity 

Detri  tus 

Juvenile  oysters 

Adult  oysters 

Predators 

Col i  forms 


Forcing  Functions 


I, 


Sun! i  ght 

River  flow 

Col i  form 

concent  rat  ion 


Numerical  value 


9.26  X  lO^m'^ 

5.84  X  lO^m^ 

2.03  X  lO^^Kcal 

7.77  X  10^  dollars 

1.41  X  lO^^Kcal 

7423  people 

7.5  X  10^  g  N 

5.2  X  lO^m^ 

1.06  X  10^°g 
5.53  X  lO^g 
8.68  X  lO^g 
3.11  X  lO^g 

5.51  X  10    coliforms 


4.5  X  lO^^Kcal/20  days 
1 .04  X  lO^m^/20  days 


2.5  X  10   coliforms/m 
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Table  8  -  continued 


Designation  on 
model  (Fig.  28a) 
and  in  Appendix 
Table  D-1 


Name 


Numerical  value 


Forcing  Functions 
(cont. ) 


M2 
^13 

^14 
Flows 

J) 


Nutrient 

concentration 

Detritus 

concentration 

Capital  investment 

Government  spending 

Goods,  fuels  and 
services 

Immigration 

Money  associated 
with  new 
residents 

Tourist  inflow 

Tourist  money  in- 
flow 

Predator  inflow 


Sunlight 

Detritus  inflow 

River  flow 

Phytopl ankton 
production 


0.25  g  N/m"" 

8.0  g/m"^ 
$3.29  X  10^/20  days 
8.2  X  10^/20  days 

1  .96  X  10^°Kcal/20  days 
20.6  people/20  days 

$6.12  X  10^/person 
576  people/20  days 

$7  X  10^/20  days 
1 .64  X  10^  g/20  days 


4.5  X  lO^^Kcal/20  days 
8.32  X  10^  g/20  days 
1 .04  X  lO^m^/20  days 

2.72  X  10^°g/20  days 
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Table  8  -  continued 


Designation  on 
model  (Fig.  28a) 
and  in  Appendix 
Table  D-1 


Name 


Numerical  value 


Flows  (cont.) 
J. 


10 


11 
12 


13 


14 


15 


16 


17 


Community 
respi  ration 

Detritus  export 

Food  uptake  by 
juvenile  oysters 

Predation  of 

juvenile  oysters 

Predation  of 
adult  oysters 


2.30  X  10^^  g/20  days 
8.50  X  10^  g/20  days 


1 .46  X  10^  g/20  days 


1  .21  X  lO'  g/20  days 

1  .90  X  10^  g/20  days 
7 


Predator  respiration  2.34  x  10  g/20  days 


Predator  stress 
respi  ration 

Bay  flushing 

River  nitrogen  input 

Nitrogen  from  towns 

Nutrient  regenera- 
tion 

Oyster  nutrient 
regeneration 

Oyster  respiration 

ii'  u  t  r  i  e  n  t  export 

Nutrient  uptake  in 
photosynthesi  s 


7.7  X  10°  g/20  days 
1  .04  X  lO^m^/20  days 
2.6  X  10^  g  N/20  days 
2.77  X  10^  g  N/20  days 


1  .28  X  lO''  g  N/20  days 


7.2  X  lO'  g  N/20  days 

1.3  X  10^  g/20  days 

1 .34  X  10^  g  N/20  days 


1 .36  X  10^  g  N/20  days 
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Table  8  -  continued 


Designation  on 
model  (Fig.  28a) 
and  in  Appendix 
Table  D-1 


Name 


Numerical  value 


Flows  (cont.) 
^21 


22 


23 


25 


101 
102 


'27 


29 
'30 
'31 


32 


'33 


34 


35 


Adult  oyster 
recrui  tment 

Oyster  harvest 

Col i form  input 
from  river 

Col i form  input 
from  towns 

Coliform  losses 


Oyster  uptake 
of  detritus 


2.91  X  10'  g/20  days 

1 .06  X  10^  g/20  days 

2.6  X  10^^ 

coliforms/20  days 

1 .89  X  10^^ 

coliforms/20  days 

2.79  X  10^^ 

coliforms/20  days 

1  .46  X  10^  g/20  days 


Microbe  respiration   2.56  x  10  g/20  days 


Microbe  nutrient 
utilization 

Estuarine  input 
to  image 

Image  depreciation 

Immi  gration 

Birth  rate 

Tourist  inflow 

Death  rate 

Tourist  outflow 

Enii  gration 


1.28  X  10°  g  N/20  days 

6.2  X  10^  Kcal/20  days 
2.95  X  10^°Kcal/20  days 
20.6  people/20  days 
6.6  people/20  days 
576  people/20  days 
4.6  people/20  days 
576  people/20  days 
18.2  people/20  days 
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Table  8  -  continued 


Designation  on 
model  (Fig.  28a) 
and  in  Appendix 
Table  D-1 


Name 


Numerical  value 


Flows  (cont.) 

^37 
•^38 


Return  to  natural 
land 


5.6  X  lO'^m^/20  days 


39 


40 


'80 


'41 


'41 


'42 


'43 


'45 


46 


'92 


'92 


Investment  return   $1.64  x  10/20  days 


Government 
spending 

Money  outflow  with 
emi  gration 

Tourist  money 
inflow 

Money  inflow  with 
immi  grati  on 

Maintenance  of 
structure 

Urban  input  to 
image 

Linear  depreciation 
of  structure 

Square  drain  depre- 
ciation of  struc- 
ture 

Income  earned  out- 
side the  county 

Income  due  to  oy- 
ster sales 

Land  development 

Loss  of  natural 
1  ands 


.2  X  10^/20  days 


$1.13  X  10^/20  days 


$7.0  X  10^/20  days 


$1.38  X  10^/20  days 
1  .89  X  lO^^Kcal/20  days 
1.89  X  lO^^Kcal/20  days 
1.5  X  10^°Kcal/20  days 


2.86  X  lO^Kcal/20  days 


$1 .02  X  10^/20  days 


$4.76  X  10^/20  days 
5.6  X  10"^m^/20  days 


5.6  X  lO^m^/20  days 
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Table  8  -  continued 


Designation  on 
model  (Fig.  28a) 
and  in  Appendix 
Table  D-1 


Name 


Numerical  value 


Flows  (cont.) 

^93' 
^93 


95 


'96 


'97 


Goods,  fuels,  and 
services  used  to 
develop  land 

Capital  outflow  for 
development  of 
new  lands 

New  land  develop- 
ment input  to 
image 

Investment  rate 

Capital  outflow  for 
goods,  fuels,  and 
servi  ces 

Oyster  harvest 
effort 


7.4  X  10°Kcal/20  days 


$2.96  X  10^/20  days 


7.4  X  10°Kcal/20  days 
$3.29  X  10^/20  days 


$7.6  X  10^/20  days 


1 .74  X  lO^Kcal/20  days 


Predator  inflow 


1.64  X  10°  g/20  days 


Fig.  29.   Justification  of  the  configuration  used  in  the 
regional  model  of  Franklin  County  and  Apalachi- 
cola  Bay  shown  in  Fig.  28a. 

Fig.  29a.   Natural  land  (Qi)  was  chosen  as  a  state  variable 
because  Franklin  County  leaders  felt  land  avail- 
able for  development  played  a  role  in  attracting 
'investment  and  people  to  the  county  and  could  in 
the  future  be  a  factor  limiting  areal  urban  growth, 
A  developed  land  storage  (Q2)  was  included  to  ac- 
count for  all  land  in  the  county.   The  double  work- 
gate  was  chosen  because  three  factors  were  sug- 
gested as  being  the  minimum  number  required  for 
urban  growth  (J92>  J93)  to  occur.   These  were 
capital  to  pay  for  fuels,  space,  and  a  source  of 
fuels  (I7).   Other  factors  which  are  required 
(labor,  information,  etc.)  were  assumed  to  be 
nonlimiting.   Pathway  J35  was  included  to  show 
the  continual  effect  of  depreciation  and  was  made 
linear  because  no  important  interactive  affects 
seemed  to  be  present.   The  price  (P])  paid  for 
goods,  fuels,  and  services  was  controlled  by  out- 
side factors  because  Franklin  County  does  little, 
even  in  the  oyster  industry,  to  affect  average 
prices  (Colberg  and  Windham,  1965).   Both  natural 
and  developed  land  storages  are  enclosed  by  the 
hexagon  group  symbol  in  recognition  of  the  self- 
maintaining  nature  of  natural  ecosystems.   The 
use  of  the  group  symbol  here  does  not  add  additional 
pathways  or  interactions  to  the  model. 
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PRICE 


■TO    IMAGE 


GOODS, 
FUELS,       ^ 
AND  -J-r* 

SERVICES 


CAPITAL 


e>TO  STRUCTURE 


TO    IMAGE 


Fig.  29b.   Image  (Q3)  is  the  long-term  ability  of  the  area 
to  attract  outside  goods,  fuels,  and  services, 
investment,  tourists,  and  new  residents.   Image 
was  chosen  to  be  a  state  variable  to  integrate 
the  diverse  factors  that  may  attract  people  and 
capital  to  the  county.   As  a  state  variable  image 
had  a  decay  time  and  could  influence  other  path- 
ways for  a  period  of  time  after  inputs  stopped. 
It  was  assumed  that  image  had  inputs  from  all 
major  work  activities  in  the  area  (J4i  .  ^26'  "^93" 
The  linear  drain  (J29)  was  included  to  show  de- 
preciation effects.   It  was  assumed  that  losses 
from  image  were  negligible  in  attracting  outside 
materials  to  the  county  because  promotion  and  ad- 
vertising were  not  large  county  activities.   No 
special  feedback  interactions  were  used  between 
image  and  image  inputs  for  this  same  reason. 
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Fig.  29j.   Juvenile  and  adult  oysters  and  predators  were  rep- 
resented by  hexagon  symbols  to  show  the  self-main- 
taining nature  of  these  components.   In  each  case 
the  full  hexagon  symbol  was  simplified  to  show  on- 
ly processes  considered  essential  to  the  simula- 
tion.  Food  input  to  juvenile  oysters  was  a  multi- 
plicative function  of  food  availability  (Qg)  and 
the  effect  of  low  salinity  in  retarding  oyster 
feeding  (l-Qs)-   While  the  exact  relationship  be- 
tween feeding  and  salinity  was  not  known,  the 
configuration  used  captured  at  least  the  extremes 
of  the  relationship  (Gunter,  1953).   An  auto- 
catalytic  feedback  was  not  shown  on  juvenile  oys- 
ters because,  even  in  poor  harvest  years,  set  of 
juveniles  was  farily  good  and  therfore  was  assumed 
to  be  nonlimiting  (Menzel,  1973).   Juvenile  oyster 
respiration  was  approximated  by  a  linear  drain 
(Jlc).   Adult  oyster  recruitment  (J19)  was  assumed 
to  Bepend  only  on  the  size  of  the  juvenile  stock. 
Predation  loss  was  assumed  to  depend  on  both  the 
biomass  of  prey  (Qiq  and  Qn)  and  predators  (Q12). 
In  the  case  of  predators,  the  autocatalytic  feed- 
back was  required  because  external  input  pulses 
of  predators  could  have  significant  impact  on  oys- 
ter stocks.   If  predator  biomass  increased  for  this 
reason,  predation  would  increase  because  of  the 
autocatalytic  arrangement  used.   Pathway  Jp  was  a 
constant  source  of  predator  biomass  coming  from 
the  Gulf  of  Mexico.   Without  the  seeding  provided 
by  Jp,  predators  would  be  eliminated  permanently 
if  they  ever  went  to  zero.   Predator  losses  in- 
cluded both  a  linear  respiratory  drain  (J4)  and 
a  multiplicative  loss  which  was  a  function  of  preda- 
tor biomass  and  freshwater  in  the  bay  (Jio)- 
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lands  (Jg3"),  natural  land  metabolism  (Joe)  (shown  to  be 
dependent  on  the  area  of  land  available),  and  primary  pro- 
duction in  the  estuary  (J27).   Image  had  one  density  de- 
pendent drain  to  account  for  losses  due  to  depreciation. 
Flows  extend  out  from  image  to  workgates  controlling  in- 
vestment money,  tourist,  and  new  resident  inputs.   Flows 
such  as  these  which  come  from  the  sensor  symbol  attached 
to  a  storage  indicate  that  the  flow  was  not  an  appreciable 
drain  from  that  tank.   A  sensor  action  causes  a  flow  along 
the  pathway  with  which  it  interacts  by  drawing  energy  from 
it. 

The  structure  storage  has  two  inputs.   The  first 
(J^^)  is  dependent  on  the  availability  of  capital,  the 
level  of  structure  present  and  the  availability  of  goods, 
fuels  and  services.   This  flow  represents  the  input  required 
to  maintain  the  level  of  structure  and  to  add  new  elements 
to  the  present  structure.   The  second  input  (Jgo)  which 
is  dependent  on  natural  lands  and  capital  represents  new 
structure  which  is  added  to  the  county  when  new  lands  are 
developed.   Both  of  these  work  flows  have  inputs  to  image. 
Additionally,  as  new  structure  is  added  along  pathway  Jq- 
natural  lands  are  used  (J92'^  ^^^   converted  to  developed 
lands  (Jg2)-   Developed  lands  revert  to  their  natural  state 
along  pathway  J^g.   Losses  from  the  structure  storage  in- 
clude a  density  dependent  drain  (J42)'  ^  square  drain  (J43) 
which  becomes  more  important  as  structure  becomes  larger 
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and  more  complicated,  and  a  loss  which  is  proportional  to 
the  effort  involved  in  the  oyster  harvest.   External  labor 
is  indicated  by  a  sensor  pathway  (J^^). 

Inputs  to  the  capital  storage  include  a  direct  in- 
put from  government  sources  (J30)  (net  input  beyond  taxes 
paid),  and  from  external  labor  markets  where  input  was  de- 
pendent on  price  (P^)  paid  for  labor  (J45)-   Money  also 
entered  with  new  residents  (Jog)  and  tourists  (J4Q)  to  the 
county  and  by  investment  spending  in  the  county  (Jgc)' 
Lastly,  money  entered  in  exchange  for  oysters  harvested 
(J,g)  with  an  externally-fixed  oyster  price  (P^).   Money 
left  the  county  with  emigrating  residents  (Jog),  interest 
on  investments  and  repayment  of  principal  on  loans  (J07). 
purchases  of  goods,  fuels,  and  services  used  for  growth 
and  maintenance  (Jqc)*  and  for  the  development  of  new  land 
(Jgo').   Prices  (P,  and  Pp)  were  also  controlled  externally 

Inputs  to  the  population  storage  include  new  resi- 
dents (J3g),  tourists  (Joo)'  and  births  (Joi)-   Losses 
include  death  (J33).  tourists  leaving  the  county  (J34)  .  and 
residents  emigrating  from  the  county  (^35)- 

Inputs  to  the  water  and  plankton  system  include 
coliforms  (Jo-i).  nutrients  (nitrogen)  (J-jp).  and  detritus 
(J,)  which  enter  in  proportion  to  the  input  of  freshwater 
from  the  Apalachicola  River  (Jp).   River  flow  varies  sea- 
sonally (Fig.  14).   Coliforms,  nutrients,  and  detritus 
concentrations  were  assumed  constant  in  each  run  but  could 
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be  changed  between  simulation  runs.   Sunlight  (J^)  also 
varied  seasonally  with  high  and  low  values  in  May  and  De- 
cember, respectively  (Fig.  15).   Predators  were  also  intro- 
duced from  the  Gulf  of  Mexico  (J  ). 

P 

The  amount  of  freshwater  in  the  bay  is  indicated  by 
storage  Qg.   When  this  tank  is  full,  the  bay  has  a  salinity 
of  0.0  o/oo.   This  tank  is  also  used  as  an  index  of  turbidi- 
ty, the  turbidity  being  high  when  the  salinity  is  low.   This 
relationship  was  based  on  data  from  Estabrook  (1973).   The 
storage  has  one  input  from  the  river  (J^)  and  one  density- 
dependent  drain  (J-ii). 

The  nitrogen  storage  (Q^)  has  inputs  from  the  river 
(J^^)  urban  population  (J^^)'^"^  recycle  from  respiration 
(J^4  and  J^g).   Losses  from  this  tank  include  uptake  in 
photosynthesis  (J^g),  bacterial  utilization  (Jnn?^'  ^"^ 
flushing  from  the  bay  (J-it). 

The  coliform  storage  (Q^3)  has  inputs  from  the  river 
(J2^)  and  urban  population  (J22)  a"^  losses  due  to  flush- 
ing and  death  (J23)-   The  number  of  coliforms  in  the  bay 
was  used  to  control  a  logic  switch  (L)  regulating  oyster 
harvest.   When  coliform  numbers  exceeded  the  Public  Health 
Service  standard  (70  MPN)  the  switch  closed  and  adult  oy- 
sters were  not  harvested.   Harvest  was  possible  again  when 
coliform  numbers  fell  below  the  legal  limit.   The  logic 
switches  of  the  model  provide  for  harvest  when  coliforms 
are  below  the  threshold  and  oyster  harvesting  is  in 
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season  (S).   Oyster  harvesting  was  permitted  from  September 
through  May. 

The  detritus  storage  (Qg)  has  inputs  from  the 
river  (J^)  and  phytopl ankton  production  (J3).   Losses  in- 
clude microbial  respiration  (J^oi^'  ouster  feeding  (Jg), 
use  by  other  herbivores  (J^),  and  flushing  to  the  Gulf  of 
Mexico  (J5).   Phytoplankton  production  was  proportional  to 
the  product  of  sunlight,  nutrients,  and  turbidity  (J3). 

The  oyster  reef  system  includes  juvenile  oysters 
(Q^q),  adult  oysters  (Q^^),  and  oyster  predators  (Q,^). 
Input  to  juvenile  oysters  (all  oysters  less  than  3")  came 
from  the  detritus  storage  (Ji5)-   Losses  include  respira- 
tion (J^g)  and  predation  (J^).   The  adult  oyster  storage 
receives  input  from  the  juvenile  oysters  via  recruitment 
and  has  losses  to  predators  (Jg)  and  harvest  (Jon^-   Oyster 
predators  receive  input  from  both  juvenile  and  adult  oysters 
and  have  losses  due  to  respiration  (Jg)  and  freshwater 
stress  (J^o^'   Pi^edator  biomass  is  seeded  from  the  Gulf  of 

Mexico  (J  ). 
P 

Assembled  Data  on  Franklin  County 


Data  used  in  evaluating  the  model  in  Fig.  28a  were 
given  in  Table  8  with  calculations  and  assumptions  in  Ap- 
pendix Table  D-1 .   Data  with  time  sequences  are  given  in 
Figs.  30-34.   Given  in  Fig.  30  are  graphs  of  the  annual 
consumption  of  electric  power,  gasoline,  and  kerosene. 


Fig.  30.   Yearly  consumption  of  fuels  in  Franklin  County, 
(a)  Electricity,  (b)  Gasoline,  (c)  Kerosene. 
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Electric  power  use  has  increased  steadily  over  the  period 
1958-1972  with  some  indication  of  an  increasing  rate  of 
use  in  1972.   Consumption  of  electric  power  almost  tripled 
between  1950-1970.   Gasoline  use  doubled  in  the  five-year 
period  prior  to  1972.   Kerosene  use  was  small  compared  to 
gasoline  and  electricity  use  and  annual  use  has  declined 
in  recent  years.   Time-series  data  were  not  available  for 
bottled  gas,  but  bottled  gas  use  was  estimated  to  be  approx- 
imately 7.5  X  10^  gal/yr  in  1972  (Stapleton,  1973). 

Figure  31  shows  trends  in  the  yearly  number  of 
visitors  to  Florida  and  Franklin  County  (Florida  Tourist 
Study,  1970).   Unfortunately,  tourism  data  for  Franklin 
County  have  not  been  recorded  for  a  long  period  of  time  and 
after  1970  were  no  longer  recorded  on  a  county  basis  for 
any  county  in  Florida.   The  1970  estimate  of  about  11,000 
visitors  to  the  county  probably  underestimates  total  visi- 
tors because  of  special  events  such  as  the  Apalachicola 
Seafood  Festival,  which  attracts  thousands  of  people  for 
a  short  period  of  time.   Estimated  tourist  expenditures  in 
the  county  amounted  to  $1.7  x  10  /yr  in  1970. 

Figure  32  shows  annual  trends  in  population,  total 
and  per  capita  personal  income,  and  gross  sales  in  Franklin 
County  (Florida  Statistical  Abstract,  1973).   County  popu- 
lation has  not  changed  dramatically  in  the  last  twenty 
years.   Both  small  increases  and  decreases  were  recorded 
during  this  period.   From  1900  to  1960  there  was  a  35% 


Fig.  31.  Annual  number  of  visitors  to  (a)  Florida,  (b) 
Franklin  County,  and  (c)  a  summary  of  average 
length  of  visit  and  expenditures. 
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SUMMARY      OF     AVERAGE      LENGTH       OF     VISIT      AND        EXPENDITURES 


AVERAGE      LENGTH      OF     STAY     PER      PERSON  12.4    DAYS 

AVERAGE       EXPENDITURE     PER      PERSON      PER     DAY  $12.63 

AVERAGE      EXPENDITURE     PER      PERSON      PER     STAY  $  156.61 


(c) 


Fig.  32.   Urban  trends  in  Franklin  County, 
(b)  Total  and  per  capita  income. 


(a)  Population, 
(c)  Gross  sales. 
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increase  in  population  (0.6%/yr)  from  4800  to  6500  people. 
Growth  in  population  in  the  county  has  been  very   slow  com- 
pared to  other  areas  of  Florida.   Both  total  personal  in- 
come and  per  capita  personal  income  have  increased  steadily 
since  1955.   Between  1965  and  1970  per  capita  income  increased 
at  a  rate  of  about  $100/yr  (8%/yr).   Per  capita  income  in 
Franklin  County  in  1970  ($1,928)  was  low  compared  to  most 
counties  in  Florida  and  was  well  below  both  the  Florida 
average  ($3,661)  and  the  United  States  average  ($3,933). 
Gross  sales  in  the  county  increased  slowly  from  1952  ($3.42 
X  10^/yr)  through  1969  ($9,24  x  10^/yr).   The  average  annual 
increase  was  about  $340,000/yr.   From  1971  through  1974 
gross  sales  almost  doubled,  going  from  $9.32  x  10^/yr  to 
$17.1  X  10^/yr  (Stalvey,  1975). 

Figure  33  gives  annual  oyster  harvest  in  terms  of 
quantity,  dockside  value  and  price  per  pound  in  1970  dol- 
lars (Rockwood,  1973).   Annual  harvest  has  varied  consider- 
ably from  year  to  year.   In  general,  annual  harvests  were 
larger  in  the  past  10  years  than  in  previous  periods.   Dock- 
side  value  of  the  oyster  harvest  has  increased  considerably 
since  1930.   In  general,  value  tracked  the  harvest.   Price 
per  pound  for  oysters  (in  1970  dollars)  increased  by  a  fac- 
tor of  about  3  between  1930  and  1971.   Prices  were  virtually 
constant  from  1960  to  1968  at  about  $0.35/lb  (meat  weight 
at  dockside).   Prices  increased  slowly  from  1968  to  1971 
reaching  about  $0.43/lb  in  1971.   Total  harvest  declined 
during  this  same  period. 


Fig.  33.   Annual  oyster  harvest  data  for  Franklin  County  in- 
cluding quantity,  dockside  value,  and  price  per 
pound  in  constant  1970  dollars. 
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Figure  34  shows  annual  landings  and  dockside  value 
of  finfish,  shrimp,  and  blue  crabs  landed  in  Franklin  Coun- 
ty.  Finfish  harvest  includes  both  food  fish  and  trash  fish 
(Division  of  Marine  Resources,  1970).   The  sharp  peaks  in 
landings  in  1952  and  1958  were  the  result  of  large  menhaden 
landings.   In  recent  years  the  dockside  value  of  finfish 
landings  has  been  less  than  $100,000/yr.   Total  shrimp 
landings  have  varied  considerably  from  year  to  year,  often 
by  as  much  as  a  million  pounds  between  successive  years. 
Since  1950  shrimp  landings  have  averaged  about  2  x  10^  lbs/ 
yr.      Dockside  value  prior  to  1970  averaged  about  $500,000/ 
yr.      Since  then  dockside  value  has  increased  to  about 
$1  ,500,000/yr  in  1972.   Blue  crab  landings  have  also  varied 
considerably  from  year  to  year.   From  1962  through  1972 
landings  averaged  about  1  x  10^  Ibs/yr  and  had  a  dockside 
value  of  about  $60,000/yr. 

A  great  deal  of  additional  information  was  required 
to  evaluate  the  regional  model  above  and  beyond  that  given 
in  this  section.   Data  for  which  no  time  series  information 
was  available  are  given  in  Appendix  Table  D-1. 

Simulation  Results  and  Comparisons  with  Assembled  Data 

Shown  in  Fig.  35a  and  35b  are  graphs  from  simula- 
tions of  the  regional  model  of  Franklin  County  and  Apala- 
chicola  Bay  with  starting  conditions  calculated  for  1970. 
Seasonal  trends  from  field  data  are  plotted  as  dots  on  the 


Fig.  34.   Annual  landings  and  dockside  value  of  major  com- 
mercial species  in  Franklin  County.   (a)  Finfish, 
(b)  Shrimp,  (c)  Blue  crabs. 
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Fig.  35.   Simulation  results  from  the  model  in  Fig.  28a  us- 
ing transfer  coefficients  calculated  from  1970 
data.   Data  points  are  from  Figs.  14  and  15  and 
Appendix  Table  D-1.   (a)  Inputs  of  river  and  sun- 
light and  responses  for  estuarine  storages.   (b) 
Responses  within  urban  sectors.   All  storages  re- 
fer to  amounts  in  areas  which  are  on   nontimber 
company  lands. 
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diagram.   Forcing  functions  of  sunlight  and  river  flow  were 
set  3  months  out  of  phar^e  with  peaks  in  June  and  March,  re- 
spectively.  Peak  sunlight  input  was  5000  Kcal/m^/day. 
Peak  river  flow  corresponds  to  50,000  cfs.   Minimum  points 
for  sunlight  and  river  flow  occurred  in  December  and  Septem- 
ber, respectively,  and  were  2400  Kcal/m^/day  and  12,000  cfs. 
Notice  that  most  estuarine  variables  showed  large  seasonal 
variations.   Urban  variables,  with  the  exception  of  capital, 
showed  no  noticeable  seasonal  variations.   Seasonal  changes 
in  capital  were  caused  by  open  and  closed  seasons  for  oyster 
harvest.   The  steady  state  pattern  shown  for  urban  variables 
is  in  agreement  with  qualitative  observations  of  trends  in 
the  county. 

Model  responses  to  both  increased  and  decreased 
concentrations  of  detritus  entering  the  bay  via  river  flow 
are  shown  in  Fig.  36.   This  modification  was  examined  to 
anticipate  possible  effects  of  changing  detritus  concentra- 
tions due  to  land-use  changes,  increased  urban  inputs,  or 
damming  of  the  Apalachicola  River.   Detritus,  nitrogen, 
juvenile  oysters  and  predators  all  showed  higher  seasonal 
maxima  and  minima  under  conditions  of  increased  detritus 
input.   Decreased  minima  and  maxima  were  observed  for  simu- 
lation runs  in  which  detritus  concentrations  were  reduced 
below  the  control  level.   Notice  that  under  conditions  of 
decreased  detritus  input,  adult  oyster  biomass  increased 
slightly.   However,  when  detritus  input  was  increased  to 


Fig.  36.   Simulation  results  from  the  model  in  Fig.  28a 
with  six  different  rates  of  detritus  input  as- 
sociated with  river  flow.   Normal  input  rates 
are  given  in  Appendix  Table  D-1.   (1)  normal 
rate,  (2)  twice  normal  rate,  (3)  three  times 
normal  rate,  (4)  one-half  normal  rate,  (5)  one- 
fifth  normal  rate,  (6)  no  detritus  input. 
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2  and  3  times  the  normal  input  rate,  adult  oyster  biomass 
developed  fairly  large  seasonal  changes  and  both  peak  and 
minimum  adult  oyster  biomass  was  depressed  below  control 
levels.   Under  conditions  of  detritus  loading  the  larger 
than  usual  pulse  in  young  oyster  biomass  stimulated  a  rapid 
increase  in  predator  biomass  as  soon  as  river  flow  began  to 
decline.   When  adult  oyster  biomass  began  to  increase, 
predator  biomass  was  quite  large  and  predation  was  heavy 
and  occurred  when  the  recruitment  of  juvenile  oysters  into 
the  adult  stock  was  approaching  a  seasonal  minimum.   The 
seasonal  pattern  of  capital  was  depressed  slightly  under 
conditions  of  detritus  loading. 

Model  responses  to  increased  and  decreased  concen- 
trations of  nitrogen  entering  the  bay  via  river  flow  are 
given  in  Fig.  37.   This  modification  was  examined  in  anti- 
cipation of  possible  nitrogen  concentration  changes  due  to 
changes  in  agricultural  runoff  and  increased  urban  inputs 
of  nitrogen.   In  general,  storages  were  not  very   sensitive 
to  either  increased  or  decreased  nitrogen  input  rates  and, 
except  for  the  nitrogen  storage,  were  similar  to  responses 
observed  in  Fig.  35  (control  conditions).   Note  that  while 
the  river  was  the  major  source  of  imported  inorganic  nitro- 
gen, most  of  the  nitrogen  used  in  photosynthesis  was  sup- 
plied via  nitrogen  recycled  from  respired  detritus  (com- 
pare pathways  J^^  and  J^^  in  Appendix  Table  D-1).   Thus, 
if  the  composition  estimates  used  were  correct,  the  original 


Fig.  37.   Simulation  results  from  the  model  in  Fig.  28a  for 
the  response  of  estuarine  variables  and  capital  to 
three  different  rates  of  nitrogen  input  associated 
with  river  flow.   Normal  input  rate  is  given  in 
Appendix  Table  D-1.   (1)  normal  rate,  (2)  twice 
the  normal  rate,  (3)  no  input  of  nitrogen. 
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source  of  most  inorganic  nitrogen  used  in  this  system  en- 
tered as  nitrogen  bound  in  organic  matter. 

Shown  in  Fig.  38  are  time  responses  of  all  state 
variables  to  two-fold  increases  in  the  input  rates  of  coli- 
forms,  nitrogen,  and  detritus.   This  modification  was  ex- 
amined because  changes  in  water  quality  often  include  sev- 
eral factors.   Because  of  possible  interactive  effects,  the 
outcome  under  conditions  including  several  pollutants  may 
be  different  than  when  impact  is  examined  with  one  pollutant 
at  a  time.   In  this  case,  the  response  of  the  estuarine 
storages  (nitrogen,  detritus,  juvenile  and  adult  oysters, 
and  predators)  was  very   similar  to  those  seen  in  Fig.  36 
under  conditions  where  the  detritus  input  rate  was  twice 
the  control  rate.   However,  a  two-fold  increase  in  coliform 
input  was  sufficient  to  exceed  the  state's  legal  coliform 
limit.   In  this  simulation,  this  resulted  in  closure  of  the 
fishing  grounds  in  January.   Thus,  oysters  were  harvested 
and  sold  for  only  a  part  of  the  season.   Because  of  this, 
capital  reached  a  lower,  but  still  seasonally  varying 
steady  state.   The  low  point  for  capital  in  this  run  was 
about  20%  lower  than  in  the  control  run  (Fig.  35).   in  re- 
sponse to  the  decline  in  capital,  image,  structure,  and 
population  declined  to  lower  steady  states. 

Figure  39  shows  adult  oyster  and  predator  biomass 
under  conditions  of  no  oyster  harvest  and  no  predation. 
This  modification  was  examined  to  show  the  relative  impact 


Fig.  38.   Simulation  results  from  the  model  in  Fig.  28a  for 
the  response  of  all  state  variables  to  a  two-fold 
increase  in  the  input  rate  of  coliforms,  nitrogen, 
and  detritus.   River  flow  was  not  adjusted  from 
the  normal  condition. 
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Fig.  39.   Simulation  results  from  the  model  in  Fig.  28a  for 
the  response  of  adult  oysters  and  predators  to  var- 
ious conditions  of  harvest  and  predation.   (a)  nor- 
mal harvest  and  predation,  (b)  no  harvest,  normal 
predation,  (c)  normal  harvest,  no  predation. 
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of  harvest  and  predation  on  the  adult  oyster  stock.   When 
oysters  were  not  harvested,  oyster  biomass  increased  only 
slightly,  but  peak  predator  biomass  increased  by  50%.   How- 
ever, when  predation  was  removed,  oyster  biomass  increased 
by  approximately  a  factor  of  two,  this  indicating  the  domi- 
nant effect  of  predation  in  controlling  oyster  stocks. 

Figure  40  shows  the  response  of  state  variables  to 
changes  in  both  river  flow  volume  and  river  pulse  amplitude. 
There  were  three  types  of  river  flow  modifications  examined. 
These  included  changing  the  pulse  height  (cases  2-5),  elim- 
inating the  pulse  (cases  6-8),  and  eliminating  river  flow 
entirely  (case  9).   Case  1  was  the  normal  condition.   Plots 
of  river  flow  curves  are  shown  at  the  top  of  Fig.  40.   No- 
tice that  adult  oyster  biomass  increased  with  river  flow 
patterns  that  had  seasonal  pulses  smaller  than  the  control 
pulse  but  that  had  approximately  the  same  amount  of  total 
freshwater  input  (cases  5  and  6).   However,  river  flows 
much  higher  (7)  and  much  lower  (cases  9  and  8)  than  the 
normal  case  produced  smaller  adult  oyster  stocks.   Smaller 
than  normal  predator  stocks  were  observed  in  cases  5,  6, 
and  7  due  to  freshwater  stress.   With  no  river  flow,  pred- 
ator biomass  was  small  because  their  food  source  was  small. 
Whenever  the  juvenile  oyster  stock  and  predator  stock  were 
above  normal,  the  adult  stock  was  below  normal.   Model  re- 
sponses for  capital,  residents,  and  structure  followed  the 
response  of   adult  oysters.   When  adult  oyster  biomass 


Fig.  40.   Simulation  results  from  the  model  in  Fig.  28a  with 
varying  flow  and  amplitude  of  river  pulse.   Total 
yearly  river  flow  volume  and  pulse  height  are  giv- 
en below  as  percentages  of  average  river  flow  and 
pulse  height.   (a)  Estuarine  storages,  (b)  Urban 
storages. 

Yearly  volume.    Pulse  height. 
Flow  regime      %   of  average      %  of  normal 
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increased,  each  of  the  above  variables  increased.   In  the 
control  case  these  variables  were  in  steady  state.   Even 
slight  decreases  in  oyster  biomass  caused  declines  in  urban 
variables. 

The  response  of  adult  oyster  biomass,  capital  and 
structure  to  three  different  adult  oyster  recruitment  rates 
is  shown  in  Fig.  41.   This  modification  was  examined  both 
because  recruitment  rates  may  vary  and  because  earlier  sim- 
ulations indicated  that  all  urban  variables  were  sensitive 
to  changes  in  the  adult  oyster  stock.   In  this  simulation 
the  price  paid  for  oysters  and  for  goods,  fuels  and  services 
was  inflating  at  10%  per  year.   Note  that  increased  recruit- 
ment rates  increased  the  seasonal  variation  in  adult  oyster 
biomass  while  reduced  recruitment  rates  leveled  seasonal 
changes  and  depressed  the  standing  stock  of  oysters.   Capi- 
tal followed  a  similar  pattern.   Note  also  that  when  prices 
paid  for  oysters  (an  income  source  to  the  county)  and  for 
goods,  fuels,  and  services  (a  cost  to  the  county)  both  in- 
creased at  an  annual  rate  of  10%,  county  structure  went 
into  a  gradual  decline  and  was  about  30%  lower  after  10 
years.   The  decline  in  structure  was  more  rapid  when  oyster 
recruitment  was  at  a  tenth  of  the  normal  rate. 

Shown  in  Fig.  42  are  model  responses  to  four  dif- 
ferent rates  of  government  spending  in  the  county.   In- 
creased government  spending  (cases  4  and  5)  increased  the 
capital  in  the  county  and  caused  small  increases  in  structure 


Fig.  41.   Simulation  results  from  the  model  in  Fig.  28a  for 
the  response  of  (a)  adult  oysters,  (b)  capital, 
and  (c)  structure  to  three  different  rates  of 
adult  oyster  recruitment.   Prices  for  oyster  sales 
and  goods,  fuels,  and  services  purchases  were 
both  inflating  at  10%  per  year.   (1)  normal  re- 
cruitment rate,  (2)  5  times  normal  recruitment 
rate,  (3)  one-tenth  normal  recruitment  rate. 
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Fig.  42.   Simulation  results  from  the  model  in  Fig.  28a 

for  the  response  of  (a)  capital,  (b)  structure, 
(c)  population,  and  (d)  adult  oysters  to  four 
different  rates  of  government  spending  in  the 
county.   All  prices  were  constant.   (1)  normal 
rate  of  government  spending,  (2)  no  government 
spending,  (3)  5  times  the  normal  rate  of  govern- 
ment spending,  (4)  10  times  the  normal  rate  of 
government  spending. 
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and  population.   In  case  5  there  was  a  25%  increase  in 
structure  and  tripling  of  population  over  the  5  years  the 
model  was  run.   Increased  government  spending  resulted  in 
lower  adult  oyster  biomass  because  of  increased  harvest 
effort.   The  changes,  however,  were  small. 

Shown  in  Fig.  43  are  model  responses  to  four  dif- 
ferent rates  of  tourist  and  tourist  money  inputs.   Changes 
in  capital,  population  and  structure  were  small  for  an  in- 
crease to  5  times  the  normal  rate  and  likewise  declines  were 
small  when  tourist  entry  was  zero.   However,  exponential 
growth  occurred  when  tourist  and  tourist  money  input  was 
increased  by  a  factor  of  10  beyond  the  control  case.   Note 
that  in  year  5  capital  takes  on  a  saw-tooth  shape  as  a  re- 
sult of  partial  oyster  industry  closure  due  to  excessive 
coliform  levels.   Note  also  that  money  earned  from  oyster 
sales  was  only  a  small  part  of  the  capital  stock  when  the 
partial  closure  occurred.   This  simulation  assumed  that  all 
external  urban  sources  were  in  infinite  supply  and  that  all 
prices  were  constant. 

Figure  44  shows  simulation  results  for  capital, 
structure,  and  image  under  conditions  where  the  rate  of 
investment  money  entering  the  county  is  increased  but  only 
for  one  year.   The  investment  rate  was  readjusted  to  its 
control  level  after  one  year.   This  modification  was  ex- 
amined because  large,  but  short-term  investments  in  tourist 
facilities  may  occur  in  the  county.   Could  such  an  investment 


Fig.  43.   Simulation  results  from  the  model  in  Fig.  28a 

for  the  response  of  (a)  capital,  (b)  population, 
and  (c)  structure  to  four  different  rates  of 
tourist  and  tourist  money  inflow.   All  prices 
were  constant.   (1)  normal  rate,  (2)  no  tourist 
or  tourist  money  inflow,  (3)  5  times  the  normal 
rate  of  tourist  and  tourist  money  inflow,  (4) 
10  times  the  normal  rate  of  tourist  and  tourist 
money  inflow. 
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initiate  long-term  growth?   Note  that  all  variables  under 
both  cases  (2  and  3)  of  increased  investment  responded  with 
increased  levels.   The  response  of  capital  was  the  sharpest 
to  increased  investment  both  in  increasing  and  returning  to 
control  levels.   The  increases  observed  for  structure  and 
image  were  smaller,  took  longer  to  reach  a  peak  and  declined 
to  control  conditions  more  slowly  than  capital  stocks.   In 
this  case,  large,  but  short-term  injections  of  capital  did 
not  result  in  permanent  growth  or  permanently  higher  levels 
of  structure  capital  or  image.   If  the  investment  continued 
and  inflation  was  not  severe,  growth  did  occur  as  shown  in 
Fig.  42. 

Shown  in  Fig.  45a  are  model  responses  when  prices 
paid  for  oysters  (an  income  source  for  the  county)  and 
prices  paid  for  goods,  fuels,  and  services  (a  cost  to  the 
county)  were  both  increasing  at  10%  per  year.   Notice  that 
over  the  ten-year  period  inflation  caused  a  decline  in 
structure  of  about  50%.   Peak  levels  in  capital  declined 
only  slightly  but  the  seasonal  change  in  capital  increased 
reaching  a  minimum  lower  than  observed  in  the  control  case. 
Inflation  caused  more  rapid  turnover  of  capital  because 
both  income  and  expenses  increased.   Because  of  this,  the 
capital  storage  became  sensitive  to  interruptions  in  money 
flows  such  as  those  associated  with  the  seasonal  closure 
of  the  oyster  industry.   The  seasonal  change  in  adult  oy- 
ster biomass  increased  slightly  but  the  average  stock 


Fig.  45 


Simulation  results  from  the  model  in  Fig.  28a  for 
the  response  of  capital,  structure,  predators  and 
adult  oysters  to  an  inflation  rate  of  10%  per  year 
for  both  purchases  and  sales  and  changes  in  har- 
vest effort.   (a)  normal  harvest  effort,  (b)  twice 
normal  harvest  effort. 
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remained  about  the  same.   In  Fig.  45b  a  doubling  in  harvest 
effort  reduced  the  decline  in  structure  to  30%  over  the  ten- 
year  period  but  was  not  sufficient  to  maintain  original 
stocks.   The  peak  level  of  capital  was  the  same  as  in  the 
control  case  but  the  seasonal  range  in  capital  increased  be- 
yond that  observed  in  Fig.  45a.   Both  adult  oyster  and  preda- 
tor biomasses  decreased  by  factors  of  about  two  below  those 
recorded  for  the  control  case.   Adult  oyster  biomass  under 
conditions  of  twice  the  normal  harvest  effort  became  a  lim- 
iting factor.   This  simulation  suggests  that  attempting  to 
increase  work  efforts  in  periods  of  inflation  will  not  stimu- 
late growth  or  even  maintain  past  levels  if  energy  sources 
are  limiting  (oysters  in  this  case).   Such  a  plan  may  lead 
to  instability  as  indicated  by  the  capital  storage. 

In  Fig.  46  harvest  effort  was  increased  to  four 
times  the  control  effort.   Sales  and  purchases  were  both 
inflating  at  a  rate  of  10%  per  year.   Under  these  conditions 
structure  still  showed  some  decline  (20%  over  ten  years). 
However,  the  seasonal  range  in  capital  reached  "boom  and 
bust"  proportions  and  it  may  be  doubtful  that  an  economy 
could  operate  under  such  extreme  conditions.   Notice  that 
under  these  extreme  harvest  conditions  predators  are  vir- 
tually eliminated  much  the  same  as  pests  are  eliminated 
in  intensive  agriculture.   Oyster  stocks  developed  a  larger 
seasonal  range  and  slightly  higher  peak  levels  but  were 
still  the  major  limiting  factor.   In  Fig.  46b  harvest  effort 


Fig.  46.   Simulation  results  from  the  model  in  Fig.  28a 
for  the  response  of  structure,  capital,  preda- 
tors, and  adult  oysters  to  harvest  effort  four 
times  greater  than  in  the  control  condition.   In 
(a)  both  purchases  and  sales  had  inflation  rates 
of  10%  per  year.   In  (b)  purchases  of  goods,  fuels 
and  services  (J93',  Jge)  inflated  at  10%  per  year, 
while  oyster  sales  inflated  at  5%  per  year. 
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was  still  at  four  times  the  control  level  but  the  inflation 
rate  of  oyster  sales  was  reduced  to  5%  per  year  while  the 
inflation  rate  of  goods,  fuels  and  services  remained  at 
10%  per  year.   In  this  case  the  decline  in  structure  was 
greater  than  in  Fig.  44a  but  the  seasonal  range  in  capital 
decreased.   Adult  oyster  and  predator  biomasses  were  about 
the  same  as  in  Fig.  46a. 

Shown  in  Fig.  47a  are  simulation  results  for  struc- 
ture, capital,  adult  oysters  and  predators  under  conditions 
of  10%  per  year  price  inflation  for  both  purchases  and  sales 
but  with  increased  investment  rates.   This  modification  was 
examined  to  estimate  the  level  of  investment  required  to  off- 
set the  declines  caused  by  inflation.   Structure  and  capital 
have  patterns  similar  to  the  control  case  for  the  first  five 
years  of  the  simulation  and  then  both  show  slow  declines. 
The  seasonal  range  in  capital  was  only  slightly  larger  than 
in  the  control  case.   In  this  simulation  an  investment  rate 
ten  times  the  control  rate  offset  the  declines  observed  in 
urban  variables  in  simulation  runs  with  normal  investment 
rates  but  with  similar  inflation  rates.   Figure  47b  shows 
response  under  the  same  inflation  rate  conditions  as  in 
47a  but  with  an  investment  rate  20  times  the  control  rate. 
Under  these  conditions  some  growth  in  population  and  struc- 
ture was   observed  during  the  first  six  years  of  the  simula- 
tion.  However,  at  the  end  of  ten  years  the  seasonal  pattern 
in  capital  was  about  the  same  as  in  the  previous  case. 


Fig.  47.  Simulation  results  from  the  model  in  Fig.  28a  for 
the  response  of  variables  to  inflation  rates  of 
10%  per  year  for  purchases  and  sales  with  increased 
investment.  (a)  10  times  normal  investment  rate, 
(b)  20  times  normal  investment  rate,  (c)  20  times 
normal  investment  rate  with  temporary  3-year  clo- 
sure of  the  oyster  fishery. 
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Structure  and  population  were  higher  than  in  the  previous 
case  but  both  were  declining.   Figure  47c  shows  the  level- 
ing effect  that  the  loss  of  the  oyster  fishery  has  even 
when  the  investment  rate  was  20  times  the  control  rate. 
The  results  suggest  that  during  early  stages  of  growth, 
the  oyster  fishery  remains  a  key  feature  of  the  local  econ- 
omy.  Loss  of  the  fishery,  at  that  stage,  could  completely 
reverse  trends  as  shown. 

Figure  48  shows  the  response  of  capital,  structure, 
population,  and  image  to  an  inflation  rate  of  20%  per  year 
for  both  sales  and  purchases  and  investment  rates  of  10, 
30,  and  50  times  the  control  rate.   In  all  three  graphs, 
inflation  1 eveled  growth  and  caused  declines  to  levels  be- 
neath peak  values.   Declines  were  sharpest  in  Fig.  48c  be- 
cause of  temporary  oyster  fishery  closure  after  the  third 
year  of  the  simulation.   Closure  was  caused  by  coliform 
levels  that  exceeded  state  health  limits.   In  Figs.  48b 
and  48c  harvest  effort  increased  because  of  increased  struc- 
ture.  However,  oysters  soon  became  limiting  and  thus  oy- 
ster sales  did  not  offset  the  declines  caused  by  inflating 
prices  . 

Given  in  Fig.  49  are  long-term  plots  of  structure, 
capital,  image,  population,  adult  and  juvenile  oysters, 
and  predators.   This  simulation  was  done  to  show  that  county 
growth  could  be  supported  by  the  oyster  fishery  and  that 
growth  dependent  on  the  resource  was  stable  and  leveled 


Fig.  48.   Simulation  results  from  the  model  in  Fig.  28a  for 
the  response  of  variables  to  an  inflation  rate  of 
20%  per  year  for  purchases  and  sales  with  increased 
investment  inflows.   (a)  10  times  the  normal  in- 
vestment rate,  (b)  30  times  the  normal  investment 
rate,  (c)  50  times  the  normal  investment  rate. 
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Fig.  49.   Historical  plot  obtained  from  simulation  of  model 
in  Fig.  28a.   Initial  conditions  for  adult  and 
juvenile  oysters  and  predators  was  not  readjusted 
for  this  run.   Oyster  harvest  was  initiated  in 
1954.   Initial  conditions  for  structure,  capital, 
population  and  image  were  adjusted  downward  to 
levels  approximating  those  in  1920.   Historical 
data  are  given  below  and  plotted  as  dots  on  the 
diagram.   (a)  structure  and  capital,  (b)  image 
and  population,  (c)  adult  and  juvenile  oysters, 
(d)  predators 
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soon  after  the  fishery  was  established.   The  oyster  fishery, 
as  a  major  activity  in  the  county,  was  estimated  to  have  be- 
gun in  1954  (Fig.  33).   In  this  simulation  prices  were  held 
constant.   With  the  beginning  of  the  oyster  fishery,  all 
urban  variables  showed  increases  and  leveling  at  new  but 
higher  steady  states.   Historical  data  for  urban  variables, 
when  available,  have  been  plotted  on  the  diagram. 

After  the  beginning  of  the  oyster  fishery  the  adult 
oyster  stock  declined  by  about  20%  due  to  new  harvest  losses. 
There  was  no  noticeable  change  in  juvenile  oyster  stocks. 
Both  nutrient  and  detritus  levels  (not  shown  in  Fig.  49) 
increased  slightly  in  response  to  increased  population  which, 
in  this  model,  resulted  from  increased  inputs  of  both  coli- 
forms  and  nitrogen  in  the  bay.   Predator  stocks  decreased 
slightly.   This  was  due  to  the  new  competition  that  developed 
for  food  (oysters)  between  predators  and  human  harvestors. 
Unfortunately,  historical  data  for  ecological  variables  in 
this  area  are  very   scarce  and  it  was  not  possible  to  plot 
field  data  points  to  support  the  simulation  results  obtained 
for  these  variables. 

Energy  Evaluations  for  Franklin  County 


Presented  in  this  section  are    calculations  of  the 
energy  basis  for  the  primitive  and  present  pattern  of  ac- 
tivities in  Franklin  County,  the  oyster  fishery,  and  a 
proposed  tourist  development  on  St.  George  Island. 
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Calculation  of  Energy  Basis 


An  inventory  of  energy  flows  entering  the  county  was 
drawn  in  Figs.  50  and  51,  tabulated  for  both  the  primitive 
and  present  patterns  (Tables  9  and  10),  and  converted  to 
fossil  fuel  equivalent  flows  using  energy  quality  factors 
given  in  Table  1.   Depending  on  the  resident  energy  flows 
available  for  matching  with  fossil  fuels  (purchased  flows) 
an   energy  investment  ratio  was  calculated.   Additionally, 
an  energy-based  carrying  capacity  was  estimated.   This  was 
done  by  multiplying  the  natural  energy  flows  in  fossil  fuel 
equivalents  by  2.5  thus  giving  the  amount  of  purchased  ener- 
gies the  county  could  attract  if  it  were  to  develop  to  a 
level  equal  to  the  national  average. 

Primitive  pattern.   The  energy  basis  of  the  primi- 
tive pattern  prior  to  the  first  cutting  of  local  forests 
in  about  1838  was  modelled  in  Fig.  50  and  tabulated  in 
Table  9.   The  pathways  shown  suggest  where  input  energies 
are  used,  but  not  all  internal  pathways  are  shown.   Calcu- 
lations, assumptions  and  data  sources  are  given  in  notes 
to  Table  9.   In  the  primitive  pattern,  energy  was  contributed 
from  rain  (as  it  amplifies  photosynthesis),  wind,  waves, 
tide,  river  head,  energy  due  to  the  mixing  of  fresh  river 
water  and  ocean  water,  and  sunlight  as  it  developes  temper- 
ature gradients  between  the  ground  and  above  the  vegetation 
and  as  it  drives  photosynthesis  in  the  variety  of  ecosystems 
covering  Franklin  County  and  Apalachicola  Bav. 
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The  energy  contributions  of  rain  in  generating  a 
hydrostatic  head  and  mixing  energy  are  shown  on  the  rain 
input  pathway.   The  energy  generated  by  rain  interacting 
with  other  energy  sources  in  photosynthesis  is  shown  in 
the  production  symbols. 

In  the  primitive  condition,  the  fossil  fuel  equiva- 
lents of  all  inflowing  energy  was  3.86  x  10^^  Kcalpp^/yr. 
About  60%  of  the  fossil  fuel  equi  valent  kilocalorles  of  sunlight 
was  accounted  for  in  the  estimated  metabolism  of  natural 
ecosystems.   On  an  area  basis  the  largest  system  (coastal 
plankton)  had  the  smallest  metabolism  per  area.   The  smallest 
system  (freshwater   system)  had  the  highest  per  area  metabo- 
lism.  Most  of  the  potential  work  of  the  physical  or  chemical 
type  was  done  via  wind,  and  the  free  energy  of  mixing  fresh 
and  salt  water.   Tides,  hydrostatic  head,  waves,  and  heat 
gradient  declined  in  magnitude  in  the  order  presented. 
Combined,  they  represented  about  3%  of  the  total  inputs. 

Present  pattern.   The  diagram  with  the  present  pat- 
tern of  energy  inputs  is  given  in  Fig.  51  and  tabulated  in 
Table  10.   Total  energy  input  was  greater  in  the  present 
pattern  than  in  the  primitive  pattern,  totalling  4.58  10^^ 
Kcalpp^/yr.   Man-related  additions  of  fuels,  goods,  and 
services  were  0.721  x  10^^  Kcalpp^/yr.   Natural  metabolic 
work  flows  decreased  slightly  because  of  conversion  of 
swamp-type  forests  to   cleared/planted  lands  or   pine  planta- 
tions.  These  new  systems  may  have  higher  net  production 


OJ 

SZ 

r^ 

t/> 

(/) 

ja 

+J 

fO 

to 

c: 

r — 

3 

o 

•I— 

o 

1 — 

to 

, — 

cr 

-M  cC 

QJ 

:3 

Q. 

1 — 

C 

. 

CJ 

•f- 

c 

3 

o 

4- 

>5 

•1— 

O) 

en 

1 — 

t- 

0) 

•1— 

(U 

S- 

1/1 

c 

(/) 

<u 

>->  o 

+-> 

4- 

4- 

c 

O 

3<NJ 

O 

— 

c 

O 

O 

4->  t— 

•,— 

(O  ^ 

Cl  c 

QJ 

(O 

S- 

<u  s- 

fO 

>  u. 

•1— 

E 

+J  <u 

ra 

•r-  x: 

s- 

E  •<-> 

ai 

220 


ococ 

-J  55 

<o 

4 


^«^^d3Ald.^' 


O    C    (O 

T3    O    <U 

>> 

o  c 
o  <u 

O)  J=    Q. 

^  </> 

ra  (/)  +J 
r-  S  C 
•r-  O  <U 
Id  ■—  r— 
>  t-    (O 

<a        > 

(/)  I—  3 
+J  <  O" 
3  (U 

Q. 
c     •■— 

•r-     C    O) 
O    3 

cn  a> 

S-    (U  r— 

O)  i-  •!- 
<!;>,</) 

4->    o 

O    3 

OCM 
C  Or— 
S-  O 
<U  C  I— 
4J  T- 
4->  I—    C 

Ql  C 

^-^  s-  s- 

CO  U_    ro 

o^  <u  E 

I—  ^    (C 

en 
+->  c  to 
c  T-  •.- 
cu       -o 

(U    S-    (U 

s-  o  x: 

Q.    s  4-> 


222 


^^BAJ^^" 


223 


-•-' 

>, 

(U 

c 

>»■•-> 

S-  I— 

13 

Ol 

•r- 

O  JQ 

O 

s. 

r* 

■M    «0 

(U 

to 

O  1— 

c 

3 

lO 

Lu  cru-  a; 

c 

(U 

,_ 

CO 

e: 

i. 

(O 

>i 

i. 

r— 

u. 

(O 

<T3 

3 

(U 

C 

s- 

<u 

C 

#t 

la 

<  -ii 

+J 

s- 

r~ 

1 — 

o 

(O 

t— 

S- 

^    0) 

W»    >> 

S-     3 

+J  -^ 

O  Lu 

c    u 

3 

(U    u. 

r— 

^™  'r* 

<o  1— 

rtJ  to 

>    «0 

3    <0 

•r-    O 

C    O 

3  i^ 

C  U- 

orvj 

ct 

LUi— 

>> 


c 

<Dt^ 

<  a.  o 

■" 

M- 

9* 

o 

F 

OJ   w 

ro 

+->    <D 

CU 

tn    S- 

s- 

>^  U 

<C  00  <c 

^ 

o 

o 
o 


00 

o 


o 
o 
o 


o 


o 
o 


ro  ID 

s 


4-)  -O 

fC    (O 


i/l  ro 

O)  3 

XI  4-> 

•1-  (/I 


LlJ  CJ3 

<l 

O^ — ' 

1    o 

c 

O  -r- 

•r-     >, 

i-   4J  -U 

X    O) 

X)    ro    ro 

•1-    S- 

>>+->    (U 

s:  <u 

31    (/)  31 

224 


r— 

s- 

^ 

(U 

W  s^ 

i~ 

3 

■M-C: 

o 

U- 

C     Ll 

2 

OJ     U. 

r— 

r—      U. 

r— 

•r" 

<a  .— 

(0 

(/> 

>    (O 

3 

V) 

•r-    O 

C 

O 

3  iii 

LU 

OTVJ 
LUi— 

C 

o 

>>      ai 

>,-(->    S-r- 

cn-r-  o  ja 

S-  I—  ■!->  ro 
Qj  <o  o  I— 
c  3  m 

UJ  O-U.  <D 
(U 
t>0 


eC  ^  ^ 

S-  I— 
1—    O    n3 


00 
O 


J- 

1 

>> 

(tJ  ^ 

(U  cu 

Jbif    S-    ^ 

S-  ct    o 

<u 

O           fO 

c 

s-  -^ 

c 

ojr^ 

<i:  Q.  o 

*"■ 

«4- 

„ 

O 

F 

(U   (O 

(0 

+J  <u 

<u 

U1    i. 

s- 

>>o 

<  c/J  < 

5 

o 

>>o 

(O   LO 

3 


O  -r- 
>54->     OJ 


>—  o 

•t- 

Irt 

(O    t/l 

4-> 

CO 

'O  ■)->   E 

i_ 

F 

3  E 

t>1 

F 

+->  -i^    QJ 

CO 

a) 

^  <U 

0> 

0) 

(fl    C  -fJ 

3 

+-> 

{/)   4J 

s- 

■»-> 

03   ns   t/) 

+J 

lo 

tt)    t/5 

s- 

CO 

O  ^    >, 

to 

>>, 

s-    >^ 

<u 

>. 

O  Q.  C/5 

UJ 

1/5 

Ll_  00 

1— 

C/1 

225 


r-  $- 

S-  3  -M  -^ 
O  U-    C    LU 

3  <U    U. 

I—  1—    U_ 

1—  T-  ro  i— 
(O  10  >  rd 
3  to  T-  O 
C  O  3  i^ 
C  LU    OTVJ 

<:         LlIi— 


O 
O 


CO 
O 


+-> 

>,      <u 

c 

>,-)-J    S-r- 

3 

Cn-r-    O  JD 

o 
o 

S-  r—    +J    <d 

0)    IT3    O  1— 

C    3    (O 

loj  cyi-i.  <u 

c: 

<u 

T- 

CO 

■ — 

' 

.^ 

c: 

S- 

itj 

>> 

s- 

r^~           ^^ 

U- 

(O           (O 
3           0) 

O) 

C      •«  (O 

J= 

<  ^  --^ 

4-> 

S-  I— 
.—    O    (0 

c 

<a  3  o 

4J           i^ 

•r* 

O        CM 

o 

1—     .— 

c 

o 

s. 

' — 

(U 

0) 

I— 

+J 

o 
o 
o 


c 

0) 

— , 

<  Q.  O 

<4- 

O 

f= 

<u 

V) 

03 

+J 

<u 

(U 

10 

i. 

S- 

>, 

o 

< 

00  <a: 

^ 

o 

o 


r—  O 


•r—  to     ro 


O  -r- 

S-    4->  T3 


226 


1 — 

t. 

^ 

O) 

</l     >, 

S- 

3 

4->  <: 

o 

u. 

C     LU 

3 

(U     U. 

r— 

r—     U. 

r^ 

•r— 

(0  r— 

fO 

(/> 

>    fO 

3 

10 

•r-     O 

C 

O 

:3  v: 

C 

LL. 

Otvj 

< 

LU.— 

vo 
o 
o 


o 


00 

o 


>,4->  S-  I— 
Dl-r-  O  J3 
t.  r—  4J  (O 
Oi  (O  U  h- 
C  3  ro 
LU  CTLi.    O) 


00 
CVJ 

C) 


o 
o 


<u  <u 

^    5-  i. 

J-  <  o 

O  ro 


3        :^ 

«a:  Q-  o 


00 

CO 


>d-  >— 


M- 

•> 

O 

i^ 

a; 

m 

CO 

+-> 

OJ 

(U 

00 

S- 

S- 

>^ 

o 

ct  00  ct 

^ 

o 

o 

o 


■0(0(0 


(O  (O  l/l 
O  r-  >, 
C_>  Q.  OO 


(O    (/) 

s  e 

(/I  -M 
(1)    (/> 

i-     >5 

Ll_  I/) 


(O  uo 

2 


227 


r— 

i. 

^    <U 

</l  >, 

S_    3 

■M  -^ 

O  U. 

C     LU 

s 

(U     U. 

r^ 

1—     li_ 

1—  "r- 

IB  1— 

to    I/) 

>    (0 

3    (A 

••-    O 

C    O 

3  ^ 

C  U- 

OtVJ 

<: 

UJr— 

c 

o 

a^-i-  o  jQ 

S-  I—  -l-J  (O 
O)  (O  O  h- 
C    3    ra 

uj  cyu-  (u 
CO 


<  -i<£  ^ 

J-  1— 

1—0(0 


o 

<u  <u 

r— 

^    S-    i- 

s-  <t  u 

<u 

O           (0 

r— 

3  -M  -.^ 

^ 

•r-  r^ 

lO 

.—    E    ro 

H- 

<0  =3    O 

3          i^ 

c  i--^ 

c  (ur-» 

<  Q-  O 

' 

«»-     « 

O  E 

O)    (/) 

(T3  -M    <U 

<u  to  S- 

i.    >,  O 

<  C/)  <£ 

03 

lO 

CO 

^ 

^ 

<£) 

o 

r^ 

CM 

>— 

o 

o 

^a- 

t^ 

O 

o 

o 

o 

O 

00 


o 
o 


<0  LO 


+->   to 

00  ■!- 

o  ra 

■o  > 

CO  CJ 

O    S- 

o  <u 

o  00 

228 


Notes  to  Tables  9  and  10 


Estimated  metabolism  of  coastal  plankton  system. 
Average  metabolism  was  estimated  to  be  9  Kcal/m2/day  (McKel- 
lar,  1975). 

(9  Kcal/m^/day)(365  days/yr) (4047  m^/acre)  = 

1  .33  X  10^  Kcal/acre/yr 

Area  weighted  metabolism  estimate  for  all  estuarine 
systems.   Included  here  are    the  oligohaline,  grass  flat,  salt 
march,  oyster  reef,  sand  bar,  medium  salinity  plankton,  high 
energy  beach,  high  velocity  channels,  high  energy  beach,  and 
coastal  dune  systems.   For  metabolism  estimates  see  Young 
(1974),  Lehman  (1974),  McFarland  (1963),  and  Table  4. 

Estimated  metabolism  for  all  freshwater  systems. 
Included  here  are  the  turbid  peidmont  river,  shallow  woodland 
pond,  and  blackwater  river  systems.   Average  metabolism  was 
estimated  to  be  76  Kcal/m^/day  (Odum,  1956). 

(76  Kcal/m^/day)(365  days/yr ) (4047  m^/acres)  = 

11 .2  X  10^  Kcal/acre/yr 

Area  weighted  metabolism  estimate  for  all  terres- 
trial systems.   Included  here  are  the  pine  scrub,  high  pine- 
land,  pine-palmetto  flatweeds,  bay  swamp,  river  swamp,  and 
freshwater  marsh  systems.   For  metabolism  estimates  see 
Odum  ejt  al.  (1972),  Ovington  (1961),  Woodwell  and  Whittaker 
(1968),  Whittaker  and  Woodwell  (1969),  Westlake  (1963), 
Bayley  and  Odum  (1973),  Bayley  and  Burnes  (1974),  Cowles 
(1974)  ,  and  Young  (1974). 


Total  sunlight  was  estimated  by  multiplying  average 
yearly  sunlight  by  total  area.   Average  sunlight  input  was 
estimated  to  be  4  x  103  Kcal/m^/day  (Odum,  1971). 

(4.0  X  10^  Kcal/m^/day)(8.97  x  10^  acres) 
(4.047  X  10^  ni^/acre)(365  days/yr)  = 
5.30  x  10^^  Kcal/yr 
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Local  heat  gradient  work  estimated  by  multiplying 
average  sunlight  input  per  area  per  year  times  Carnot  ratio 
(AT  =  2°C)  times  the  total  study  area  (Odum,  1974b). 

(4000  Kcal/m^/day)(g3^_^)(365  days/yr) (4047  m^acre) 

X  (897,154  acres)  =  36.78  x  10^^  Kcal/yr 

^Potential  energy  of  tides  calculated  by  using  the 
following  formula.   Average  tidal  ampliture  was  from  Gorsline 
(1963). 

W  =  l/2pgh^A 

Where  p  =  density  (1.0  g/cm  )  2 

g  =  gravitational  acceleration  (980  cm/sec  ) 

h  =  tidal  amplitude  (46  cm)  p 

A  =  unit  surface  area  (1.0  m  ) 

P.E.  Tides= 

(1/2)(1.0  g/cm^)(980  cm/sec^)(46  cm^)(10^  cm^/m^)  x 

(2  tides/day)(2  cycl es/tide ) ( 2 .4  x  10"^^  Kcal/erg)  x 

(365  days/yr)(4047  m^/acre) ( 1 27  ,020  acres)  = 

0.190  X  10^^  Kcal/yr 


Wave  energy  absorbed  on  shoreline  was  calculated 
usinq  the  following  formula.   Average  wave  height  from  Walton 
(1974). 

Power/Length  of  Shoreline  =  l/8pg^/^H^^^ 

Where  p  =  density  (1-0  g/cm  )  2 

g  =  gravitational  acceleration  (980  cm/sec  ) 
h  =  average  wave  height  (21.3  cm) 

Wave  Energy/yr  = 

(1/8)(1.0  g/cm^)(980  cm/sec^)  "^/^  (21  .  3)^/^ 

(2.4  X  10"^^  Kcal/erg)(3.15  x  10^  sec/yr)(100  cm/m) 

(11.17  X  10^  meters  of  beach)  =  6.77  x  10^°  Kcal/yr 


230 


Chemical  free  energy  released  due  to  mixing  of 
solutions  with  different  concentration  was  calculated  using 
the  following  formula. 

C2 

AG  cal/g  solute  =  nRT  In  7?^ 

^1 

Where  n  =  number  of  moles  (35  g/mole) 

R  =  gas  constant  (1 .99) 

T  =  temperature  (300°K) 

C,  =  concentration  of  solute  in  river  water 
'    (35,000  ppm) 

Cp  =  concentration  of  solute  in  ocean  water 
^        (120  ppm) 

AG  =  (1/35)(1.99)(300)(2.3  log  ■^^^) 
=    -96. 7g  cal/g  solute 

Kcal/yr  =  (96. 7g  cal/g)(120  g/m^)(2.21  x  10^°m3/yr) 

X  (1  Kcal/lOOOg  cal)  =  25.6  x  10^°  Kcal/yr 

"^Potential  energy  of  river  head  was  calculated  using 
the  formula  given  below  (Healy,  1974). 

PE/yr  =  pghd 

Where  p  =  density  (1,0  g/cm  ) 

2 
g  =  gravitational  acceleration  (980  cm/sec  ) 

h  =  river  elevation  change  in  Franklin  County 
(45  cm) 

d  =  yearly  river  discharge  (2.21  x  10   cm  /yr) 

PE/yr  =  (1.0  g/cm^)(980  cm/sec^)(45  cm)  x 

(2.21  X  10^^  cm"^/yr)(2.4  x  lO"^^  Kcal/erg) 

=  2.34  X  10^°  Kcal/yr 
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Yearly  work  done  by  winds  was  based  on  the  kinetic 
energy  of  wind.   An  eddy  diffusion  coefficient  of  10,000 
cm^/sec  was  used.   Average  wind  velocity  was  8.7  mph(Florida 
Statistical  Abstract,  1P72).   This  wind  speed  was  assumed 
to  occur  10^  cm  above  ground. 

Work/yr  =  1,2  x  10""^  g/cm'^)(371  cm/sec)^(l  x  10^  cm^/sec)  x 
(2.39  X  10"^^  Kcal/erq)(3.15  x  10^  sec/yr)(3  x  10^^  cm^) 

(2)(1  X  10*  cm) 

=  1.87  X  10^^  Kcal/yr 

Annual  work  done  by  rain  was  divided  into  three 
categories.   The  work  done  in  photosynthesis  was  accounted 
for  in  the  terrestrial  and  freshwater  metabolism  measure- 
ments and  was  not  recounted  here.   The  potential  energy  of 
rainwater  relative  to  river  water  due  to  the  concentration 
difference  (mixing  energy)  was  calculated  using  runoff  water 
but  not  water  lost  due  to  evaporation  or  transpiration.   Po- 
tential energy  of  water  due  to  its  position  relative  to  sea 
level  was  calculated  using  the  volume  of  runoff  water. 

(a)   Mixing  energy 

^2 
AG  g  cal/g  solute  =  nRTln^r- 

^1 

where  n  =  number  of  moles  (35  g/mole) 

R  =  gas  constant  (1.99  g  cal/deg-mol e) 

T  =  temperature  (300°K) 

C-,  =  concentration  of  solute  in  river  water  (120  ppm) 

Cp  =  concentration  of  solute  in  rainwater  (1.2  ppm) 
AG  =  (^  g/mole)(1.99  g  cal/deg-mole)(300°K)(2.3  log  |||) 

AG  =  -78  g  cal/g  solute 
Runoff  flow  estimated  to  be  43  cm/yr  (Mann,  1975) 
Total  flow/yr  =  (43  cm/yr)(total  area) 

=  (43  cm/yr)(13.8  x  10^^  cm^/county) 

1 4    3 
=  5.93  X  10   cm  /county  area/yr 

8  3 
=  5.93  X  10  m  /county  area/yr 

Kcal/yr    =  (78  g  cal/g  solute)(1.2  g/m'^)(5.93  x  loV/yr) 

=  5.55  X  10^°  cal/yr 

=  5.58  X  10^  Kcal/yr 


232 

(b)   Hydrostatic  head 

PE/yr  =  1/2  pghd 

o 

where  p  =  density  (1.0  g/cm  ) 

g  =  gravitational  acceleration  (980  cm/sec  ) 
h  =  average  elevation  above  sea  1 evel  (610  cm) 
d  =  yearly  runoff  flow  (5.93  x  10^^  cm"^/yr) 
PE/yr  =  (1/2)(1.0  g/cm^)(980  cm/sec^)(610  cm)  x 

(5.93  X  10^^  cm^/yr)(2.4  x  lO"^^  Kcal/erg)  =  4.26  x  10^  Kcal/yr 

Estuarine  metabolism  was  calculated  for  the  pres- 
ent pattern  as  in  (b)  but  with  area  weighted  adjustments 
for  metabolism  of  boat  basins  and  harbors  and  navigation 
channels.   There  were  al so  some  area  adjustments  in  other 
estuarine  ecosystems  (Table  3).   Estimated  metabolism  in- 
creased only  slightly. 

Terrestrial  metabolism  was  calculated  for  the 
present  pattern  as  in  (d)  but  with  area  weighted  adjust- 
ments for  the  metabolism  of  planted  lands  and  pine  forests. 
Estimated  metabolism  decreased  slightly  compared  to  the 
original  pattern. 

Total  electricity  use  was  obtained  from  Florida 
Power  Corporation  (Crosthwai te ,  1974). 

(3.8  x  10'^  KWH/yr)(860.5  Kcal/KWH)  =  3.27  x  10^°  Kcal/yr 

Kilocalories  of  electrical  energy  were  converted  to  KcalppE 
by  dividing  by  0.28  (Table  1).  Cost  of  electric  power  was 
$9.26  X  105/yr. 

^Gasoline  use  was  obtained  from  the  Department  of 
Revenue,  Tallahassee,  Florida  (Thomas,  1973).   Cost  of  gaso- 
1 ine  was  $1  .82  x  lO^/yr. 

(5.2  X  10^  gal/yr)(127,654  BTUs/gal )(0.253  Kcal/BTU)  =  1.68  x  10^^  Kcal/yr. 

Kilocalories  of  gasoline,  bottled  gas,  and  kerosene  were 
assumed  to  be  Kcalpp£  and  needed  no  further  conversions. 

^Kerosene  use  was  obtained  from  The  Revenue  Office, 
Department  of  Agriculture,  Tallahassee,  Florida  (Rutledge, 
1973).   Cost  of  kerosene  was  $2.4  x  104/yr. 
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(1.61  X  10^  gal/yr)(141,785  BTUs/gal ) (0.253  Kcal/BTU) 

=  5.7800^  Kcal/yr 

"^Bottled  gas  use  was  estimated  by  Stapleton  (1973). 
Cost  of  bottled  gas  was  estimated  to  be  $2.33  x  lO^/yr. 

(7.5  X  10^  gal/yr)(4.24  lbs/gal ) (2. 1 7  x  10^  BTU/16)  x 

(0.253  Kcal/BTU)  =  1.75  x  10^°  Kcal/yr 

The  energy  flow  associated  with  goods  and  services 
was  estimated  by  multiplying  the  total  dollar  outflow  times 
25,000  Kcal/dollar  (Table  1). 

(Total  Dollar  Outf 1 ow) (25 ,000  Kcal/dollar) 

($16.44  X  10^/yr)(25,000  Kcal/dollar)  =  4.11  x  10^^  Kcalpp^/yr 

Energy  investment  ratio  for  present  pattern  in 
Franklin  County. 

I  =  Outside  energy  being  invested  (summation  of  all  purchased 
energy  inputs) 

R  =  Energy  flow  generated  from  resident  system  (summation  of 
all  natural  energy  inputs) 

1  2 
Investment  Ratio  =  0-721  x  10   Kcalpp^/yr  ^  Q^g 

3.86  X  10^2  Kcali:t:^/yr 
r  r  e 


Preliminary  carrying  capacity  for  present  pattern 
in  Franklin  County. 

Carrying  Capacity  =  (2 . 5)  ( natural  energy  input) 

=  (2.5)(3.86  X  10^^  Kcalpp^/yr) 

=  9.65  X  10^^  Kcalppp/yr 
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than  the  swamp  types,  but  their  gross  production  is  less 
(see  "m"  in  Notes  to  Table  9).  Other  natural  work  flows 
remained  as  in  the  primitive  condition. 

The  energy  investment  ratio  for  the  present  pattern 
was  calculated  to  be  0.19  for  the  Franklin  County  region 
(Table  10,  Note  "t").   High  quality  energy  (purchased  goods, 
fuels,  and  services)  that  could  be  attracted  to  Franklin 
County  (carrying  capacity)  to  bring  activity  to  the  national 
average  was  calculated  to  be  9.65  x  10^^  KcaUp  /yr  or  some 
13  times  the  present  input  (Table  10,  Note  "u"). 

Oyster  Fishery 


An  energy  basis  calculation  of  the  oyster  industry 
was  made  to  evaluate  the  combination  of  natural  (not  paid 
for)  and  purchased  energies  that  act  in  support  of  the  in- 
dustry.  The  calculation  also  suggested  the  amount  of  fossil 
fuel-based  activities  elsewhere  that  may  be  matched  by 
natural  energies  originating  in  Franklin  County.   Simpli- 
fied energy  circuit  diagrams  of  the  Franklin  County  oyster 
fishery  are  given  in  Fig.  52.   Figure  52a  was  evaluated  in 
energy  flows  measured  in  heat  equivalents.   Figure  52b  was 
evaluated  in  fossil  fuel  equivalents  to  show  the  ability  of 
individual  flows  to  do  work  and  amplify  low-quality  energy 
flows.   Notes  follow  the  figures  with  calculations  and  ref- 
erences . 

In  each  figure  energies  are  shown  entering  the  es- 
tuary where  they  interact  in  the  production  of  oysters  and 
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other  products.   Energies  associated  with  the  cultch  plant- 
ing program  and  sewage  (nitrogen)  are  shown  as  additional 
work  subsidizing  estuarine  production.   The  model  of  oyster 
harvest  is  shown  as  the  product  of  the  oyster  industry  and 
estuarine  oyster  stocks.   Money  is  received  for  the  oyster 
harvest  and  spent  for  goods,  fuels  and  services  for  the  in- 
dustry.  Once  the  oyster  harvest  leaves  the  county  (indi- 
cated by  the  larger  rectangle  surrounding  the  estuary  and 
oyster  fishery  components),  additional  value  is  added  to  the 
harvest  before  it  reaches  final  retail  sales  and  consumption 

In  Fig.  52b  inputs  to  the  estuary  are  dominated  by 
the  natural  energies  (1.2  x  10   Kcalpp^/yr).   The  fossil 
fuel  equivalents  provided  by  the  cultch  program  and  sewage 
effluent  are  quite  small  compared  to  the  natural  sources. 
In  the  oyster  fishery,  approximately  1.2  x  10   Kcalppr/yr  of 
natural  work  (h)  interact  with  1.68  x  10^^  Kcalpp^/yr  of 
high-quality  energies  (1,  m,  and  n).   In  the  process  about 

o 

6.13  X  10  heat  kilocalories  of  oysters  are  exchanged  for 
$6.6  X  10  /yr.   Purchased  energies  thus  represent  about 
14%  of  the  total  input  energies.   Assuming  that  this  money 
will  be  exchanged  for  a  mixture  of  goods,  fuels,  and  ser- 
vices which  have  an  energy  equivalent  similar  to  the  na- 
tional average,  then  the  oyster  income  represents  1.62  x 
10   Kcalpp^/yr.   The  energies  expended  in  harvest  returned 
almost  ten  times  the  equivalent  energies  represented  by  the 
money  paid  for  the  oysters.   Most  of  the  oysters  produced  in 


Fig.  52.   Simplified  energy  circuit  language  diagrams  show- 
ing purchased  and  natural  energy  flows  involved 
in  the  Franklin  County  oyster  industry.   (a)  En- 
ergy flows  as  heat  equivalents,  (b)  Energy  flows 
as  fossil  fuel  equivalents. 
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Fig.  52,  footnote  e 

Fig.  52,  footnote  o 

Fig.  52,  footnote  p 

Fig.  52,  footnote  q 

Fig.  52,  footnote  r 
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Notes  to  Fig.  52 


Total  natural  energy  input  acting  in  support  of  the 
oyster  reefs  in  the  bay  was  calculated  by  converting  annual 
oyster  uptake  of  organic  matter  back  to  sunlight  equival- 
ent  ki localories.   Organic  matter   kilocalories  were  con- 
verted to  gross  photosynthesis   kilocalories  by  multiply- 
ing by  10.   These   kilocalories  were  converted  to  sunlight 
kilocalories   by  multiplying  by  100  (Table  1).   Oyster  up- 
take of  organics  is  given  in  Appendix  Table  D-1 . 


Kcal/yr  =  (2. 4x1 0^ ^ Kcal /yr) ( 1 0) ( 1 00)  =  2.4x10^^ 

Energy  input  associated  with  planting  cultch  to 
catch  oyster  larvae  was  estimated  by  converting  yearly  cost 
of  this  operation  to  kilocalories  using  25,000  Kcal/$.   Data 
are  from  Whitfield  (1973). 

($10,000/yr)(25.000Kcal/$)=  2.5  x  1 O^Kcal pp^/yr 

The  energy  equivalent  of  nitrogen  entering  the  bay 
from  a  secondary  treatment  plant  was  calculated  by  convert- 
ing nitrogen  to  organic  matter  via  photosynthesis.  (Nitro- 
gen data  are  given  in  Table  D-1). 

(5.05  X  10^  g  N/yr)(20  g  organic  matter/g  N) 


(4.5  Kcal/g)  =  4.6  x  10^  Kcal/yr 

The  energy  equivalent  of  the  oyster  harvest  effort 
included  replacement  costs  of  boats,  motors,  tongs  and  other 
boat  gear,  shucking  and  packing  houses.   Data  were  from 
Rockwood  (1973). 

($3.04  X  10^/yr)(25,000  Kcal/$)  =  7.6  x  10^  Kcalpp^/yr 

q 
Expenditure  of  fuel  was  estimated  to  be  8.8  x  10  Kcali-cc/y^ 
(Rockwood,  1973).  ^^^ 

(2.75  X  10^  gal/yr)(1.27  x  10^  BTU/gal ) ( Kcal /BTU ) 


8.8  X  10^  Kcal/yr 


10 


Total  harvest  effort  =  1.64  x  10   Kcal^^^/yr 

FFE  -^ 


239 


Oyster  harvest  for  1970  was  converted  directly  to 
kilocalories  using  4.5  Kcal/g   (Table  D-1). 

(3  X  10^  lbs/yr)(4.5  Kcal/g)(454  g/lb) 

(0.1  g/g  dry  weight)  =  6.13  x  10^  Kcal/yr 

Money  earned  due  to  oyster  sales  was  calculated 
as  given  in  Appendix  Table  D-1. 

$6.5  x  10^/yr 

^Final  retail  value  of  Franklin  County  oysters  was 
estimated  from  data  given  by  Col  berg  and  Windham  (1965). 
The  money  earned  by  oyster  sales  in  the  county  (f)  rep- 
resented about  63  percent  of  final  oyster  value.   Thus, 
$10.3  X  10°  dollars  were  eventually  paid  for  Franklin 
County  oysters. 

Total  purchase  of  goods,  fuels,  and  services  by 
those  directly  and  indirectly  involved  in  the  oyster  in- 
dustry.  Only  a  portion  of  the  purchase  is  used  in  the 
harvest  effort  (compare  (h)  and  (d)).   It  was  assumed  that 
all  money  earned  was  spent  (savings  were  small  compared  to 
spending).   The  money  flow  to  the  outside  economy  was  con- 
verted to  energy  flow  using  the  ratio  of  25,000  Kcal/$. 

($5.5  X  10^/yr)(25,000  Kcal/$) 

=  16.25  X  10^°Kcalppg/yr 

^As  stated  in  the  text,  approximately  37%  of  the 
full  value  of  oysters  did  not  enter  the  county.   This  was 
estimated  to  be  $3.8  x  10°/yr  (Colberg  and  Windham,  1965). 
This  was  converted  to  energy  flow  using  the  ratio  of  25,000 
Kcal/$. 

($3.8  X  10^/yr)(25,000  Kcal/$) 

=  9.5  X  10^°  Kcalppp/yr 


^In  "heat  calories"  the  total  input  flow  was  large 
compared  to  yields  (e)  and  subsidies  (b,  c,  and  d).   Thus, 
the  heat  lost  from  the  system  was  approximated  by  the  inputs 
and  was  assigned  an  approximate  value  of  2.4  x  1 0^ '^Kcal/yr . 
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Natural  energies  (a)  were  converted  to  fossil  fuel 
equivalents  by  dividing  by  2000  (Table  1). 

2.4  X  10^^  Kcal/yr   ,  „   -_n  ^   ,    , 
o '-^^^  =    1.2  X  10   Kcalt-cr/yr 

2  X  10^  ^^^ 

The  value  given  in  (b)  for  the  cultch  planting 
program  was  assumed  to  be  in  fossil  fuel  equivalents. 

2.5  X  10^  Kcalpp^/yr 

The  energy  content  of  nitrogen  input  to  the  bay 
from  sewage  wastes  was  calculated  in  (c).   This  value  was 
at  the  gross  photosynthesis  level  of  quality  and  was  ad- 
jsuted  to  fossil  fuel  quality  by  dividing  by  20  (see  Table 


4.6  X  lO^Kcal/vr  _  ,  ^ 
20        "  ^'-^ 


X  10°Kcal 


FFE 


/yr 


The  value  given  in  (d)  for  harvest  effort  was  in 
fossil  fuel  equivalents. 

1.64  X  10^°KcaU^^/yr 
t"  r  h 


The  energy  quality  of  oysters  in  a  dispersed  state 
in  the  bay  was  calculated  by  assuming  that  10%  of  oyster 
intake  of  organics  was  assimilated  into  new  tissue. 


Energy  Quality  of  Oysters  =  assimilated  material 
(dispersed ) 


2.4  X 
2.4  X 


total  intake 
10, 


10'"Kcal/yr  _  f.    , 
11  -U.I 

10' 'Kcal/yr 


With  this  assumption  of  assimilation  efficiency,  oysters  in 
a  dispersed  state  have  an  energy  quality  about  10  times  that  of 
net  photosynthetic  production.   Net  photosynthetic  produc- 
tion is  a  factor  of  at  least  2  less  in  quality  than  fossil 
fuels.   Thus,  this  calculation  indicates  that  oysters  dis- 
persed in  nature  have  a  quality  about  five  times  that  of 
fossil  fuels. 
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^The  energy  quality  of  oysters  concentrated  to  dock- 
side  was  calculated  by  taking  the  ratio  of  organic  matter 
consumed  to  yields  (both  harvest  and  predation).   This  ratio 
was  used  because  it  reflects  the  energy  input  required  to 
produce  a  yield  which  can  be  passed  up  the  food  web  or  to 
urban  systems.   Simulation  results  in  Fig.  45  suggested  that 
the  present  harvest  rate  was  close  to  the  carrying  capacity 
of  oyster  stocks. 


Energy  Qual ity 


oysters  at 
dockside 


Yield 


Total  input 


where 

yield  =  harvest  yield  +  predation 

harvest  yield  =  6 . 1 3x1 O^Kcal/yr 
(as  in  (e)  of  this  table) 

predation  =  (5x68x1 O^g/yr ) (4. 5  Kcal/g) 

=  2.56xlO^Kcal/yr 

(Appendix  Table  D-1  for 
predation  losses). 

Total  input  =  2  .  4x1 0^ ^ Kcal /yr 

(as  in  (a)  of  this  table) 


Energy  Qual ity 

oysters 
harvested 


oysters        ^  3 . 1 7x1 O^Kcal /yr  _  ^   ^,- 

•11  -  U  .  U  1  o 

2.4x10'  'Kcal/yr 


Thus,  harvest  yield  is  related  to  input  by  a  factor  of  about 
75  /   1   \    Input  energy  (net  production)  is  lower  in  qual- 

^0.013^- 
ity  than  fossil  fuels  by  a  factor  of  at  least  2  (Table  1). 
Thus  oyster  yields  are  35  times  fossil  fuel  quality.   The 
heat  calories  of  oyster  yiel d  were  thus  multiplied  by  35  to 
adjust  them  to  fossil  fuel  equivalents. 

10, 


(6.13xlO°)(35)  =  2.15xl0"^KcaU^^/yr 

FFE  •' 


242 


^The  energy  content  of  oysters  at  the  point  of  pro- 
cessing when  they  leave  the  county  was  assumed  to  be  re- 
flected by  the  price  paid  for  the  product. 


($6.5  X  10°/yr)(25,000  Kcal/$) 
=  16.3  X 


10^°Kcalpp^/yr 


The  energy  content  of  the  oyster  after  final  retail 
sale  was  estimated  from  the  price  paid  for  the  product  (g, 
this  table)  using  the  ratio  25,000  Kcal/$. 


($10.3  X  10^/yr)(25,000  Kcal/$) 
=  25 


8  X  lO^^Kcalpp^/yr 


This  flow  represents  the  energy  flow  resulting 
from  final  processing  of  oysters  outside  of  Franklin  Coun- 
ty.  The  money  flow  was  taken  as  the  difference  between 
(q)  and  (r)  and  converted  to  energy  flow  using  25,000 
Kcal/$. 


($3.8  X  10°/yr)(25,000  Kcal/$) 
=  9.5  X 


10^°Kcalppp/yr 


This  flow  represents  the  purchase  of  goods,  fuels, 
and  services  by  the  oyster  industry.  Calculation  was  as  in 
(h)  of  this  table. 


16.3  X  10^°Kcal 


FFE 


/y'T 
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Franklin  County  are  consumed  outside  the  county.   Based 
on  value  added  estimates  made  by  Colberg  and  Windham  (1965), 
oyster  exports  ultimately  were  exchanged  for  $10.3  x  10^. 
Approximately  37%  of  this  money  did  not  enter  the  county. 
Thus,  the  eventual  use  of  oysters  supported  an  additional 
$3.8  X  10   (9.5  X  10^°  Kcalpp^)  worth  of  work  in  some  other 
area  of  the  country. 


Tourist  Development 


Energy  calculations  were  made  anticipating  the  de- 
velopment of  St.  George  Island  as  a  tourist  resort.   Energy 
calculations  were  made  to  estimate  the  impact  of  adding  this 
new  feature  to  the  full  area  of  St.  George  Island  and  to  the 
local  area  of  St.  George  Island  as  presently  proposed. 

Three  simplified  models  are  given  in  Fig.  53  show- 
ing natural  and  purchased  energy  inputs  to  Franklin  County 
under  present  conditions  and  under  conditions  that  a  tourist 
community  is  developed  on  St.  George  Island  as  proposed. 
Figs.  53b  and  53c  show  urban  inputs  associated  with  the 
tourist  development.   The  models  show  natural  energy  in- 
puts on  the  left  side  of  the  diagram,  and  purchased  sources 
entering  from  the  top  of  the  diagram.   All  flows  are  in 
fossil  fuel  equivalents. 

Natural  energy  flows  and  energy  flows  associated 
with  the  present  urban  system  are  shown  in  Fig.  53a  with 
documentation  given  in  Table  10.   Natural  energy  flows  were 


Fig.  53.   Simplified  energy  circuit  language  diagrams  showing 
aggregated  inputs  of  natural  and  purchased  energies 
and  effect  of  tourist  developments  on  St.  George  Is- 
land,  (a)  Present  pattern  of  the  county,  (b)  Pattern 
with  full  tourist  development  of  St.  George  Island, 
(c)  One  unit  of  tourist  development  with  surrounding 
ecosystems  as  defined  in  Fig.  54. 
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large   (investment  ratio  =  0.19).   New  energy  flows  asso- 
ciated with  the  proposed  island  development  were  estimated 
to  be  0.468  x  10^^  Kcalpp^/yr  (Table  11).   This  represents 
a  65%  increase  in  current  urban  energy  flows.   If  the  pro- 
posed tourist  development  were  to  expand  and  include  all 
of  St.  George  Island,  purchased  energy  flows  might  increase 
to  11.28  Kcalpp^/yr,   This  was  estimated  by  assuming  that 
the  proposed  800-acre  development  was  duplicated  to  cover 
all  20,000  acres  of  the  island.   The  purchased  flows  given 
in  Table  11  were  multiplied  by  a  factor  of  24.   Urban  en- 
ergy flows  in  the  county  increased  to  12.0  x  10   Kcalpnp/ 
yr  with  the  addition  of  the  full  island  development.   The 
investment  ratio  was  3.10  with  the  addition  of  the  tourist 
development.   This  is  comparable  to  other  developed  areas 
of  Florida  and  is  higher  than  the  national  average. 

To  obtain  an  estimate  of  the  localized  mix  of 
natural  and  purchased  energy  flows  in  the  proposed  tourist 
development,  an  area  around  the  site  was  defined  (Fig.  54). 
Figure  55  gives  a  map  view  of  the  proposed  800-acre  develop- 
ment.  It  was  assumed  that  a  five-mile  extension  on  either 
side  of  the  development  would  include  most  of  the  natural 
areas  having  exchange  with  the  development  including  fish- 
ing, vistas  and  microclimate  control.   Within  this  area 
natural  and  purchased  flows  were  compared  as  shown  in  Fig. 
53c.   Table  12  summarizes  energy  flows  for  this  area.   In 
this  case,  urban  inputs  were  very  large  compared  to  natural 
energies.   The  investment  ratio  was  9.67. 
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Notes  to  Table  11 


The  purchased  energy  flow  associated  with  the  devel- 
opment of  St.  George  Island  (800  acres;  7000  residents;  3000 
dwelling  units)  was  estimated  by  multiplying  the  dollar  cost 
(minus  advertising  costs)  times  25,000  Kcal/dollar.   Data 
were  from  the  Development  of  Regional  Impact  report  prepared 
by  Leisure  Properties  Limited  (1974).   Development  costs 
include  prorated  costs  of  both  vertical  and  hori zontal  devel- 
opment.  Advertising  costs  were  assumed  to  be  spent  outside 
the  county  and  were  not  included. 

Horizontal  development  (prorated  over  20  years) 

($.042  X  10^/yr)(25.000  Kcal/dollar) 

1.1  X  10^°  Kcalppg/yr 
Vertical  development  (prorated  over  20  years) 
($5.4  X  10^/yr)(25,000  Kcal/dollar) 

13.5  X  10^°  Kcalppp/yr 
Total  =  14.6  X  10^°  Kcalpp^/yr 

Annual  electricity  use  was  extrapolated  from  data 
developed  by  Brown  and  Genova  (1973)  for  Lee  County,  Florida. 
For  dwelling  unit  densities  of  3.5  D.U./acre,  electricity 
use  was  15  x  10^6  Kcal /D . U . /yr .   In  this  calculation  there 
Were  3.75  D.U./acre.   Over  the  twenty-year  period  considered 
here  there  was  an  average  of  2000  dwelling  units. 

(15  X  10^  Kcal/D.U./yr)(2000  D.U.) 

3.0  X  10^°  Kcal/yr 


Gasoline  use  was  calculated  using  1000  gals/yr  as 
the  average  family  consumption  (Ford  Foundation,  1974). 

(1000  gal/D.U./yr)(2000  D.U.)(3.23  x  10^  Kcal/gal) 

6.46  x 

The  energy  flow  associated  with  goods  and  services 
was  extrapolated  from  data  developed  for  Lee  County,  Florida 
by  Brown  and  Genova  (1973).   They  estimated  goods  and  ser- 
vices use  to  be  5.6  x  107  Kcal/D.U./yr  for  new  housing  at 
similar  densitites. 


10^°  Kcalppp/yr 
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(5.6  X  10^  Kcal/D.U./yr)(2000  D.U.) 


1.12  X  10 


11 


Kcalpp^/yr 


Urban  energy  flows  that  would  develop  based  on 
the  new  tax  money  from  the  island  development  were  calcu- 
lated by  multiplying  estimated  tax  receipts  (Leisure 
Properties  Limited,  1974)  by  25,000  Kcal/dollar. 

($1.52  X  10^/yr)(25,000  Kcal/dollar) 


3.8  X  10 


10 


Kcalppp/yr 


Fig.  54.   Map  showing  site  of  tourist  development  on  St. 
George  Island.   The  area  used  in  calculating 
natural  energy  input  to  the  development  (Fig. 
53c)  is  also  shown  (13,400  acres). 
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Energy  Quality. of  Water  at  Jim  Woodruff  Par 


The  methodology  for  calculating  an  energy  quality 
ratio  was  presented  in  the  Methods  section.   Given  here  is 
a  calculation  estimating  the  energy  quality  of  the  hydro- 
static head  of  water  at  Jim  Woodruff  Dam  on  the  Apalachi- 
cola  River  relative  to  fossil  fuel  quality  work. 

Fig.  56a  and  56b  are  simplified  diagrams  of  the 
process  of  converting  the  potential  energy  of  river  head 
into  electricity  at  a  hydroelectric  station.   Required  pur- 
chased subsidies  enter  the  conversion  process  at  the  top  of 
the  box  symbol  and  include  the  prorated  costs  of  construc- 
tion and  maintenance.   Land  costs  were  evaluated  in  terms 
of  diverted  metabolism  and  are  shown  as  another  energy  in- 
put in  Fig.  56b.   Degraded  heat  is  shown  at  the  bottom. 
Annual  output  of  electricity  is  shown  on  the  right.   All 
energy  flows  in  Fig.  56a  are  in  heat  kilocalories  with  no 
adjustment  to  fossil  fuel  equivalents.   Total  energy  input 
amounts  to  4.51  x  10   Kcal/yr.   Output  was  2.05  x  10^^ 
Kcal/yr  as  electrical  energy. 

In  Fig.  56b  all  flows  except  (b)  have  been  con- 
verted to  fossil  fuel  equivalents.   Gross  output  was  7.38 
X  10   Kcalpp^/yr.   After  subtracting  the  feedback  costs 
from  the  gross  output,  net  output  amounted  to  5.92  x  10^^ 
Kcalpp^/yr.   River  head  calories  are  related  to  fossil  fuel 
calories  by  the  ratio  of  net  output  to  input  (5.92  x  10^^ 
Kcal/yr  divided  by  4.14  x  10^^  Kcalppp/yr).   Thus  1.42 


Fig.  56.   Simplified  energy  circuit  language  diagrams  show- 
ing the  conversion  of  potential  energy  of  river 
head  to  electricity.   Letters  refer  to  explana- 
tory notes  in  Table  13.   (a)  Model  evaluated  in 
heat  calories,  (b)  Model  evaluated  in  fossil  fuel 
equivalent  ki 1  oca! ories . 
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5.92 


=  1.42 


INPUT  4.14 

Therefore     one     kilocolorie     of     fossil     fuel     work     is     equivalenf     to     0.7 
kilocalories    of     hydrostatic     head. 


(b) 
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Kcal  of  fossil  fuel  work  are  equal,  1n  work  ability,  to  1 
Kcal  of  river  head-type  work.   Data  and  calculations  are 
given  in  Table  13. 
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Table    13 
Documentation   of   energy    flows    given    in    Fig.    56. 


a)  Power   plant   output    in    KWH    (yearly   average) 

238,090,000  KWH 

238,090,000  KWH  x  860.5  Kcal/KWH  =  2.049  x  10^^   Kcal/yr 

A   factor   of   3.5   was    used    to   convert   electricity   to 
fossil    fuel    equivalents    (Table    1). 

(2.049  X  10^^  Kcal/yr)(3.6)  =  7.38  x  10^^   Kcalpp^/yr 

(Apalachicola  River  Basin  Reservoir  Regulation  Manual,  1972) 

b)  Calculation  of  input  energy 

net  head,  full  load  =  30.5  ft. 
turbine  discharge  (cfs) 

1  unit  operating  5,500 

2  units  operating  12,000 

3  units  operating  18,300 

Average  annual  flow  past  dam  =  21,311  cfs 

(Apalachicola  River  Basin  Reservoir  Regulation  Manual,  1972) 

Calculating  power  developed  from  this  flow  with  a  head  of 
30.5  ft. 


Pw  = 


HW 


where  H  =  height  (ft) 

W  =  weight  (lbs) 

T   =  time  interval 

Pw  =  power  (ft-lbs/sec) 


260 


Table  13  -  continued 


V  =  volume  of  water 

w  =  weight  (lbs/ft"^) 

Pw  =  HwQ 

where  Q  =  =r  =  volumetric  flow  rate  in  ft  /sec  out 
of  the  reservoir 

Pw  =  HwQ 

=  (30.5  ft)(62.3  lbs/ft^)(21,311  ft^/sec) 

=  4.049  X  10^  ft-lbs/sec 
Pw  =  4.049.10''  ft-lbs/sec  x  1.945  x  10"^  ft-l'bs/sec 

=  7.975  X  10^  Kcal/min 
Pw  =  7.875  X  10^  Kcal/min  x  60  min/hr  x  24  hr/day  x  365  days/yr 
Pw  =  4.139  X  10^^  Kcal/yr 
(Healy,  1974) 

c)  Lifetime  of  1 ock-dam-hydrof aci 1 i ties  was  estimated  to 
be  50  yrs    (Harrison,  1974). 

d)  Calculation  of  construction  costs: 

Powerhouse  $15,392,120 

Lock  8,940,191 

Spillways  and  Dikes   26.906,434 

$51 ,238,745   (1953  dollars) 

The  dam  and  hydrofaci 1 i ty  appear  to  operate  independent 
of  lock  and  because  of  this,  the  cost  of  the  lock  was 
not  included  in  this  calculation  (Snow,  1974). 

Total  cost  in  1953  dollars  =  $42,966,243 

This  figure  was  adjusted  to  1973  dollars  by  multiplying 
by  1.6.   This  was  done  so  that  the  25,000  Kcal/dollar 
ratio  could  be  used. 
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Table    13    -    continued 


Total  fossil   fuel   investment  =  $68,745,988  x  25,000  Kcal/dollar 

=  1.719  X  10^^  Kcalpp^ 

Prorated  cost  over  50  years  was  0.344  x  10      Kcalrcp/y' 

e)  Calculation  of  yearly  maintenance   costs: 

Powerhouse  $83,199/yr 

Spillway   and    Dikes  52,981/yr 

$136,180/yr 

Converting    to    Kcal/yr   using    25,000    Kcal/dollar 

$136,180/yr   x    25 ,000/dol  1  ar   =   0.034   x    10^^    KcaU^^/yr 

r  r  h 

(Snow,  1974) 

f)  Calculation  of  costs  associated  with  modified  environ- 
ments (conversion  of  forests  to  reservoir). 

Land  costs  =  (forest  loss/yr  +  river  loss/yr)  - 

(reservoir  gain/yr) 

Loss  of  forest  metabolism  =  (33.5  x  10  acres)  x 

(20  X  10^  Kcal/yr)(4.046  x  loV/acre) 

=  27.11  X  10^^   Kcal/yr 

Loss    of   river   metabolism   =    (4.0  x  10    acres)  x 

(10  X  10^  Kcal/yr) (4.046  x  loV/acre) 

=  1.62  X  10^^   Kcal/yr 

Gain    in    reservoir   metabolism   =  (37.5  x  10     acres)  x 

(4.6  X  10-^  Kcal/yr)(4.046  x  loV/acre) 


6.98  X  10^^   Kcal/yr 


Differential    cost    =    (27.11   x  10       +  1.62  x  10     ) 
(6.98  X  10^^)  =  21.75  x  10^^   Kcal/yr 
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Table  13  -  continued 


The  differential  cost  was  converted  to  fossil  fuel 
quality  by  dividing  by  20  (Table  1). 

.11 


21  .75  X  10 
20 


=  1  .09  X  10 


11 


Kcalpp^/yr 


Metabolism  estimates  are  from  Table  10, 


DISCUSSION 

This  study  considered  the  energetic  basis  of  a 
coastal  region  using  calculations,  simulation  models,  and 
areal  maps  for  purposes  of  adding  to  our  understanding  of 
the  behavior  of  large  open  systems  of  man  and  nature. 
The  study  was  also  a  part  of  a  continuing  effort  to  de- 
velop a  generalized  methodology  for  anticipating  regional 
trends  and  deciding  which  of  many  alternative  patterns  now 
facing  coastal  zone  and  other  areas  adds  to  the  vitality 
of  the  region.   The  discussion  is  presented  in  sections 
organized  by  issues  which  are  important  in  the  Franklin 
County-Apalachicola  Bay  region  and  which  may  have  broader 
applications. 

Main  Energy  Characteristics  of  the  County 

The  most  striking  characteristic  of  the  energy 
basis  of  the  county  was  the  major  role  played  by  natural 
energy  inputs.   In  1973  approximately  84%  of  all  energy 
inputs  to  the  county  came  from  natural  inputs.   The  re- 
maining 16%  was  associated  with  goods,  fuels  and  services 
for  which  payment  was  made.   The  ratio  of  urban  to  natural 
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inputs  (investment  ratio)  was  0.19.   Total  energy  input 
to  the  region  was  4.5C  x  10^^  KcaUp^/yr,  or  about  5.1 
X  10  Kcalppp/acre/yr. 

The  investment  ratio  and  the  energy  flow  per  area 
can  be  compared  to  other  systems  to  gain  perspective  on 
the  energy  basis  of  Franklin  County.  Zucchetto  (1975)  re- 
ported an  investment  ratio  of  3.90  for  the  Miami-Dade 
County  region  of  Florida.   The  ratio  for  the  South  Florida 
region  was  2.4  in  1974  while  the  overall  ratio  for  the  U.S 
in  1974  was  2.5  (Odum  and  Brown,  1975).   Franklin  County 
would  have  to  process  about  13  times  the  present  amount  of 
purchased  fuels  to  have  an  investment  ratio  similar  to  the 
national  average.   On  a  unit  area  basis,  the  Miami-Dade 
area  had  6.7  x  10^  Kcal pp^/acre/yr,  and  the  U.S.  1.7  x  10^ 
Kcalpp^/yr.   On  an  area  basis  energy  input  in  the  Franklin 
County  region  was  less  than  the  Miami-Dade  and  South  Flor- 
ida areas  but  greater  than  the  U.S.  average.   Most  energy 
inputs  to  Franklin  County  were  natural  inputs. 

Thus,  the  notion  of  Franklin  County  as  a  rural 
area  (Colberg  and  Windham,  1965)  heavily  dominated  by 
natural  energy  flows  appears  to  be  substantiated  by  quan- 
titative energy  input  calculations.   Inspection  of  maps 
from  aerial  photographs  gave  a  similar  impression.   Domi- 
nant on  the  maps  were  the  terrestrial,  estuarine,  and 
coastal  plankton  systems.   Freshwater  systems  were  less 
dominant  but  still  evident  in  the  deltaic  system  at  the 
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mouth  of  the  Apalachicola  River.   Highly  metabolic  urban 
systems  were  a  small  feature  (0.3%  of  the  total  area). 

In  establishing  the  energy  basis  of  the  study  area 
energy  flows  were  evaluated  as  they  crossed  the  boundary 
into  Franklin  County  region.   Raw  energy  flows  (represent- 
ing heat  content)  were  converted  to  fossil  fuel  equivalent 
kilocalories  to  allow  comparisons  between  flows.   On  a  fos- 
sil fuel  equivalent  basis,  the  largest  inflow  in  the  primi- 
tive pattern  was  sunlight  and  this  accounted  for  about  69% 
of  the  total  input.   Tides,  waves,  mixing  energy,  and  hydro- 
static head  were  primarily  associated  with  the  estuary  and 
amounted  to  about  1.0  x  10^^  Kcalpp^/yr  or  about  26%  of 
the  total  input.   Wind  was  also  a  major  natural  input  and 
amounted  to  about  5%  of  the  total.   Several  observations 
concerning  the  energy  input  pattern  are  offered.   First, 
estuaries  have  been  often  characterized  as  highly  dynamic 
areas  (Ketchum,  1972).   In  the  study  area  the  estuary  occu- 
pied approximately  16%  of  the  total  area.   If  sunlight  and 
wind  inputs  were  prorated  so  the  estuary  received  16%  of 
the  total  input  and  this  was  combined  with  the  estuarine 
inputs  given  above,  then  the  estuary  received  about  38% 
of  total  energy  inputs.   Secondly,  metabolism  estimates 
were  made  of  each  of  the  dominant  systems  (Table  10). 
Total  metabolism  was  estimated  to  be  about  1.5  x  10^^ 
Kcalpp^/yr.   This  figure  represents  some  portion  of  the 
work  done  by  input  energies  after  they  enter  the  system. 
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Total  ecosystem  metabolism  represents  at  least  40%  of  the 
input  energies.   Actual  work  done  by  other  inputs  was  not 
estimated  inside  the  study  area  but  some  actions  are  evi- 
dent, especially  in  the  estuary  where  shorelines  are  build- 
ing and  eroding.   The  point  here  is  that  ecosystem  metabo- 
lism may  be  the  major  factor  in  shaping  and  controlling 
landscapes  by  utilizing  input  energies. 

Estimates  of  primitive  and  present  patterns  of 
energy  input  were  calculated  for  the  Franklin  County  region 
(Tables  9  and  10)  to  quantify  the  impact  of  present  urban 
systems  on  the  full  energy  budget.   In  the  Franklin  County 
region  urban  flows  were  small  compared  to  natural  inputs. 
In  this  area  input  energies  in  the  present  pattern  are 
larger  than  those  in  the  primitive  pattern.   This  is  con- 
sistent with  Lotka's  Principal  of  Maximum  Power  (Lotka, 
1922;  Odum,  1971a),  which  states  that  successful  systems 
are  those  that  maximize  the  use  of  energy  inputs  available 
to  them.   Since  the  present  urban  area  seems  to  be  a  per- 
sistent part  of  the  landscape,  we  would  expect  that  total 
energy  flow  would  now  be  greater,  even  though  some  portion 
of  the  original  system  has  been  displaced. 

Given  in  Fig.  57  is  a  plot  of  energy  flow  per  area 
for  five  major  ecosystem  types  versus  the  percent  of  the 
total  study  area  occupied  by  each  ecosystem  type.   The 
largest  spatial  system  (coastal  plankton)  had  the  smallest 
energy  flow  per  unit  area  per  year.   Community  metabolism 


Fig.  57.   Plot  of  energy  flow  per  area  per  year  versus  per- 
cent of  total  study  area  in  each  ecosystem  type. 
Urban  energy  flow  per  area  is  in  KcalFFE/y^-   Eco- 
system energy  flows  are  in  Kcal/acre/yr  measured 
at  the  gross  photosynthesis  level.   Total  area  was 
8.97  X  10^  acres. 
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was  used  as  an  index  of  work  done  in  the  system.   The  urban 
system  which  was  small  in  area  had  the  largest  energy  flow 
per  area  per  year.   The  dashed  line  connecting  the  points 
has  the  general  shape  of  an  exponential  curve.   Odum  (1974a) 
has  suggested  that  exponential  energy  power  spectra  may  be 
developed  in  all  systems  because  of  the  requirement  to 
maximize  total  system  power.   In  maximizing  total  power, 
systems  develop  feedbacks  from  storages  to  capture  incoming 
energy.   Special  feedbacks  are  developed  that  are  of  high 
quality  to  provide  services  that  lower  quality  storages 
are  not  capable  of  delivering.   In  generating  the  high 
quality  storages  energy  is  lost  and  the  total  quantity  of 
high  quality  energy  is  less.   In  this  way  exponential  power 
spectra  are  developed. 

Curves  similar  to  the  one  in  Fig.  57  have  been  ob- 
served in  other  systems  including  food  webs  as  is  suggested 
in  data  reported  by  McKellar  (1975).   The  exact  significance 
of  this  distribution  on  a  regional  scale  is  not  clear  at 
this  point  but  it  might  be  speculated  that  all  systems  are 
characterized  by  such  relationships,  and  that  regional  sys- 
tems are  distinguished  by  such  patterns.   The  figure  is 
presented  here  as  one  way  to  view  a  relationship  between 
energy  flow/unit  area  in  major  ecosystems  and  the  areal 
extent  of  the  ecosystems  in  a  regional  pattern. 
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Relationship  of  Economy  to  Natural  Systems 

As  shown  in  the  consideration  of  the  energy  basis 
of  the  county,  natural  systems  provide  the  main  energy  sup- 
ports to  the  county.   In  Franklin  County  most  of  the  money 
earned  comes  from  fishery  products  which  at  present  are 
taken  from  wild  stocks.   The  total  value  of  the  oyster  har- 
vest was  calculated  to  be  $6.5  x  10^  (Colberg  and  Windham, 
1965)  in  1970  and  is  worth  more  at  the  present  time.   Dock- 
side  value  of  the  remaining  fishery  products  was  about  $1.3 
X  10  .   The  full  value  to  the  county  was  somewhat  more  as 
this  figure  does  not  include  value  added  in  processing. 
Tourism  in  the  county  appeared  to  be  mainly  related  to  rec- 
reational opportunities  including  fishing  in  estuarine, 
freshwater,  and  offshore  areas,  hunting,  and  relaxation. 
The  main  attractions  for  these  are  provided  by  natural  eco- 
systems with  little  or  no  cost.   Total  tourist  income  was 
estimated  to  be  $1.3  x  10^/yr  in  1970.   This  may  have  been 
an  underestimate  as  the  money  spent  by  visitors  at  the 
annual  seafood  festival  was  not  included.   Fishery  and 
tourism  incomes  were  $9.1  x  10^  in  1970,  about  60%  of  the 
total  county  income  of  about  $16  x  10  . 

Feedback  Required  of  Urban  Systems 


Figure  3  shows  a  simplified  model  of  Franklin 
County  and  Apalachicola  Bay.   The  diagram  suggests  main 
feedback  relationships  operating  between  urban  and  natural 
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systems.   These  included  the  addition  of  nitrogen  and  coli- 
forms  to  the  bay  from  the  urban  area  and  harvest  effort 
which  may  have  the  effect  of  rearranging  oyster  bar  mater- 
ials and  in  this  way  stimulating  new  growth.   On  a  larger 
scale,  exports  of  fishery  products  to  areas  outside  the 
county  have  a  return  flow  as  nutrients,  detritus,  and  col- 
iforms  from  the  river.   Only  some  unevaluated  portion  of 
these  materials  are  from  urban  areas.   Also,  there  is  a 
state-supported  program  of  shell  planting  in  the  bay  for 
purposes  of  developing  new  oyster  reefs  or  improving  the 
condition  of  existing  reefs  (Whitfield,  1974).   Connecting 
flows  from  the  estuary  include  fishery  harvest,  contribution- 
to  image  which  attracts  tourists,  residents,  and  capital, 
and  assimilation  of  nitrogen  and  col i forms  in  sewage  ef- 
f 1 uent. 

To  test  one  aspect  of  the  feedback  reinforcement 
theory  the  sensitivity  of  oyster  stocks  to  increased  har- 
vest effort  was  estimated.   As  shown  in  Figs.  45  and  46, 
increasing  harvest  effort  by  factors  of  2  and  4  times 
normal  depressed  the  stock  of  adult  oysters.   If  the  model 
is  correct,  it  would  seem  that  additional  urban  feedbacks 
to  the  estuary  and  the  oyster  reef  system  would  be  required 
to  increase  the  yield.   Rockwood  (1973)  has  proposed  several 
methods  for  improving  fishery  yields.   Whitfield  (1973)  has 
reported  results  of  a  shell  planting  program  in  Apalachi- 
cola  Bay.   Both  sets  of  suggestions  represent  potential 
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feedbacks  from  the  urban  sector  to  the  producing  system. 
Each  suggestion  should  be  evaluated  to  test  the  net  yield- 
ing characteristics  of  the  proposal.   This  has  been  done 
for  the  oyster  industry  in  this  paper. 

Energy  Value  of  Natural  Pulses 


Over  the  course  of  the  year  there  is  a  gradual 
rise  and  fall  in  sunlight,  temperature,  river  flow,  pro- 
ductivity and  other  factors  in  Apalachicola  Bay.   It  has 
been  hypothesized  elsewhere  (Odum,  1971b)  that  estuarine 
communities  have  adapted  periods  of  reproduction  and  growth 
to  coincide  with  periods  of  peak  production  in  the  estuaries 
Estuaries  characteristically  have  populations  of  fish  and 
shellfish  moving  in  and  out  on  a  seasonal  program  that 
seems  related  to  estuarine  production.   Such  migrations 
have  been  recorded  in  Texas  estuaries  by  Hellier  (1962), 
Copeland  (1965),  and  Odum  (1967),  and  in  Apalachicola  Bay 
by  Livingston  (1975)(Fig.  18). 

There  may  be  special  energy  value  related  to  puls- 
ing seasonal  programs  in  that  populations  that  were  dis- 
persed over  wide  ocean  areas  become  concentrated  in  estua- 
rine areas  and  are  available  for  harvest.   Systems  with 
pulsing  inputs  may  be  characterized  by  also  having  pulsing 
net  yields,  some  of  which  are  available  to  urban  systems. 
The  seasonal  pulse  in  water  and  materials  from  the  Apa- 
lachicola River  may  be  a  major  factor  organizing  this 
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characteristic.   Because  of  the  seasonal  pulse  in  river 
flow  bay,  salinity  and  m^ny  other  environmental  factors 
change  seasonally.   Relatively  few  species  are  adapted  to 
the  full  range  of  conditions  found  in  the  estuary  and  as 
a  result  many  species  are  present  for  a  short  period  of 
time  when  conditions  are  favorable  and  are  then  replaced 
by  others.   A  few  species  such  as  oysters  and  some  clams 
have  adapted  to  estuarine  conditions  and  are  often  present 
in  great  numbers.   The  river  pulse  which  changes  environ- 
mental condition  may  favor  a  system  with  fewer  components, 
less  spatial  diversity,  and  large  numbers  of  a  few  species, 
The  special  value  of  pulses,  such  as  river  flow,  may  be 
in  organizing  ecosystems  with  low  diversity  and  large 
stocks  of  only  several  species. 

This  concept  was  incorporated  into  the  regional 
model  (Fig.  28a).   Model  responses  were  tested  to  changes 
in  both  the  amount  of  freshwater  entering  the  bay  and  the 
amplitude  of  the  seasonal  pulse.   Results  are  shown  in 
Fig.  40.   In  this  simulation,  pulsing  river  flows  tended 
to  maintain  oyster  stocks.   River  flows  much  higher  and 
lower  than  normal  and  without  seasonal  pulses  caused  oy- 
ster stocks  to  decline.   Menzel  et  al.  (1957)  have  shown 
that  predator  stocks  increased  on  one  reef  in  Apalachicola 
Bay  when  the  salinity  pattern  stabilized  for  several  years. 
Adult  oysters  disappeared  from  the  reef  during  this  time. 
While  no  quantitative  value  has  been  given  to  the  pulse 
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of  the  Apalachicola  River  in  the  study  area,  the  ability 
of  the  urban  area  to  process  purchased  goods,  fuels,  and 
services  and  to  maintain  the  established  estuarine  oyster 
fishery  was  shown  in  this  model  to  be  related  to  river  pulse 


Alternative  Factors  Controlling  Patterns  in  the  Region 

This  section  brings  together  calculations  and  simu- 
lation model  results  concerning  county  trends  resulting  from 
proposed  or  existing  modifications.   Included  are  consider- 
ations of  nutrient  and  organic  matter  loading  with  oxygen 
depletions,  salinity  control  in  the  bay  with  passes  and  im- 
pact on  the  county  energy  budget  of  adding  urban  housing 
developments  to  St.  George  Island. 

Nutrient  and  Organic  Matter  Loading  with  Oxygen  Depletions 

The  concentration  and  seasonal  cycles  of  nitrate 
and  ammonia  were  shown  in  Fig.  17  (Estabrook,  1973).   Peak 
nitrate  concentrations  occurred  in  February  and  were  gen- 
erally lower  than  those  for  other  river-dominated  estuaries 
including  the  upper  and  mid-Chesapeake  Bay  (Carpenter  et 
aj_.  ,  1969),  Patuxent  River,  Maryland  ( Fl  emer  et  aj .  ,  1970) 
and  the  Ythan  estuary  (Leach,  1971).   In  the  Potomac  River, 
nitrate  plus  nitrate  concentrations  were  often  5  to  10 
times  the  peak  value  observed  in  Apalachicola  Bay.   The 
Potomac  receives  sewage  wastes  from  the  Washington,  D.C. 
area.   Peak  values  were  similar  to  those  observed  for 
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North  Carolina  estuaries  (Thayer,  1969).   The  seasonal 
cycle  in  nitrate  concentration  was  similar  to  those  ob- 
served in  other  river-dominated  estuaries  (Carpenter  et  al . , 
1969)  and  temperate  coastal  areas  with  early  spring  over- 
turn (Martin,  1965).   Ammonia  concentrations  were  similar 
to  those  reported  for  the  lower  portion  of  the  Chesapeake 
Bay  (Patten  et_  aj^.  ,  1963)  and  lower  than  those  reported  for 
the  Patuxent  River,  Maryland  (Flemer  ejt  aj[.  ,  1970),  which 
has  nutrient  additions  from  a  sewage  treatment  plant. 

Total  organic  carbon  was  measured  over  several  days 
in  the  fall  of  1973.   Data  were  not  available  on  a  seasonal 
basis  for  organic  carbon.   Organic  matter  concentrations  in 
the  bay  were  higher  than  those  in  the  river  apparently  be- 
cause phytoplankton  stocks  were  more  developed  in  the  bay. 
The  values  observed  were  lower  than  those  observed  in 
Georgia  marsh  creeks  (Odum  and  de  la  Cruz,  1967),  the 
upper  Chesapeake  Bay  (Flemer  and  Biggs,  1971),  and  Chinco- 
teague  Bay,  Maryland,  which  does  not  have  a  strong  river 
influence  (Boynton,  1974).   Stocks  of  organic  carbon  were 
about  5  times  higher  than  those  reported  for  shelf  waters 
of  the  Gulf  of  Mexico  (Fredericks  and  Sacket,  1968).   In 
general,  nitrogen  and  organic  carbon  levels  were  not  high 
compared  to  other  coastal  areas  with  river  inputs.   Total 
input  may  be  much  greater  than  indicated  by  the  one  set  of 
measurements  presented  here. 
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Calculations  used  in  models  indicated  that  most 

nitrogen  used  in  phytopl ankton  photosynthesis  originally 

entered  the  bay  as  organic  nitrogen  bound  in  particulate 

or  dissolved  organic  matter.   The  nitrogen  from  this  source, 

when  recycled,  constituted  about  85%  of  the  total  input  to 

the  nitrogen  storage  in  June  (compare  J,,  +  J^^  and  J^o  in 

14     lb       Id 

Table  8).   Inputs  from  urban  areas  constituted  only  about 
1%  of  the  inorganic  nitrogen  entering  the  bay. 

The  diurnal  activity  model  was  evaluated  with  late 
summer  data  from  a  station  in  the  medium  salinity  plankton 
system  which  includes  a  major  portion  of  the  bay  and  may  be 
representative  of  average  conditions  in  the  bay.   Summer 
data  were  used  because  temperature  and  salinity  were  high 
and  flushing  and  the  saturation  level  of  oxygen  were  at  a 
seasonal  low.   With  high  temperature,  oxygen  use  in  respira- 
tion may  have  been  near  a  seasonal  high.   Thus  it  may  be 
that  the  late  summer  period  is  most  likely  to  experience 
low  oxygen  levels.   Oxygen  measurements  taken  over  24-hr 
periods  during  the  summer  of  1973  never  fell  below  3.9  ppm 
even  in  the  East  Bay  where  respiration  was  dominant.   In 
general  most  stations  had  oxygen  levels  above  saturation 
during  the  day  and  slightly  below  saturation  at  night  or 
just  before  dawn.   Simulation  model  output  showed  normal 
diurnal  oxygen  variations  of  about  2.8  ppm  for  the  water 
column.   Oxygen  levels  in  the  model,  as  in  the  field,  did 
not  fall  below  saturation  even  in  the  predrawn  period. 
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The  model  indicated  the  possibility  of  nitrogen 
being  a  factor  limiting  photosynthesis  (see  gross  photo- 
synthesis plot  for  cases  1  and  2  in  Fig.  23).   Nitrogen 
additions  may  stimulate  photosynthesis  some,  increasing 
daily  oxygen  pulses  as  shown  in  Fig.  24.   However,  this 
model  was  not  very  sensitive  to  nitrogen  additions,  and 
additions  up  to  5  times  the  normal  rate  did  not  cause  oxy- 
gen levels  to  drop  below  saturation.   Estabrook  (1973) 
found  that  phosphate  may  have  been  limiting  phytopl ankton 
production  during  periods  of  calm  weather  in  late  summer. 
However,  he  suggested  that  phytopl ankton  photosynthesis  was 
primarily  controlled  by  temperature.   While  temperature  may 
stimulate  photosynthesis  some  directly  (Williams  and  Mur- 
doch, 1966)  it  also  may  be  promoting  the  recycle  of  bound 
nutrients  via  respiration.   The  large  diurnal  changes  in 
nitrogen  seen  in  all  simulations  suggest  recycle  to  be  a 
major  source  of  nitrogen  during  the  summer.   Iverson  (1975) 
has  suggested  that  there  may  be  a  substantial  stock  of 
nitrogen  in  the  bay  sediments.   If  this  is  true,  remineral- 
ization  of  this  material  would  be  another  source  of  inor- 
ganic nitrogen  available  for  phytopl ankton  photosynthesis. 
A  large  sediment  nitrogen  storage  would  also  act  as  a 
stabilizing  factor  in  the  model.   Diurnal  changes  in  in- 
organic nitrogen  would  have  been  reduced  if  this  storage 
had  been  included. 
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Oxygen  concentrations  fell  below  saturation  where 
detritus  input  was  increased  and  when  reaeration  was  re- 
stricted.  When  reaeration  was  reduced  to  10%  of  the  normal 
level  and  detritus  input  increased  to  4  times  the  normal 
input  rate  (equivalent  to  river  concentrations  of  about 
20  g  C/m  )  predawn  oxygen  levels  were  at  about  3  ppm. 
These  conditions  might  be  expected  for  areas  with  higher 
than  normal  organic  matter  inputs  and  restricted  circula- 
tion.  If  organic  inputs  increased  and  circulation  became 
restricted,  some  portions  of  the  bay  might  experience  se- 
vere oxygen  depletions  before  low  oxygen  levels  were  char- 
acteristic of  the  bay  in  general.   If  organic  matter  inputs 
to  the  bay  began  to  increase  the  East  Bay  area  may  show 
signs  of  oxygen  depletion  first  because  the  area  is  some- 
what restricted  with  poorer  circulation  and  reaeration. 
Diurnal  metabolism  measurements  indicated  that  respiration 
was  already  dominant  in  this  area. 

Salinity  Control  with  Dams  and  Passes 

The  changing  salinity  regime  in  Apalachicola  Bay 
is  an  important  factor  in  maintaining  estuarine  charac- 
teristics.  Seasonally  the  average  bay  salinity  changed 
from  about  8  ppt  in  late  winter  to  about  21   ppt  in  late 
summer.   Indications  are  that  the  seasonal  change  in  salin- 
ity due  to  river  pulse  favors  oyster  communities  over  more 
diverse  communities  with  less  yield  found  in  areas  with 
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stable  salinities.   Several  environmental  modifications 
have  been  made  or  are  planned  that  may  affect  the  bay  salin- 
ity pattern. 

In  1965  Jim  Woodruff  Dam  was  constructed  at  Chatta- 
hoochee, Florida,  on  the  Apalachicola  River.   Lake  Seminol« 
was  formed  behind  the  dam.   Monthly  average  river  flow  for 
the  period  1929-1972  was  given  in  Fig.  13.   No  obvious 
changes  in  river  flow  are  evident  from  this  plot.   Both  be- 
fore and  after  dam  construction  there  were  years  of  high 
and  low  flow.   Figure  14b  shows  monthly  river  flow  averaged 
for  a  15-year  period  before  and  after  dam  construction. 
Average  peak  flow  prior  to  dam  construction  was  slightly 
lower  than  after  the  dam  was  built.   Also,  before  the  dam 
was  built  average  river  flow  declined  to  low  flow  more 
slowly  than  after  the  dam  was  built  and  then  rose  more 
quickly.   Statistical  tests  of  the  means  for  average  monthly 
flows  before  and  after  dam  construction  (February,  March, 
and  November,  when  differences  between  average  flows  were 
the  largest)  showed  no  statistical  difference  at  the  0.01 
level  . 

Information  developed  in  this  study  and  information 
from  Estabrook  (1973)  and  Livingston  (1975)  suggest  that 
there  may  have  been  some  changes  in  the  salinity  pattern 
in  the  bay  related  to  the  construction  of  Bob  Sikes  Pass 
(Fig.  1).   Dawson  (1955a)  conducted  a  hydrographic  survey 
prior  to  pass  construction  and  reported  that  the  bay  was 
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well  mixed  with  salinity  stratification  occurring  only  oc- 
casionally in  navigation  channels  and  near  inlets.   In  sur- 
veys done  recently,  salinity  stratification  was  a  general 
feature.   Even  during  periods  of  high  river  flow  there  was 
stratification  in  the  outer  portions  of  the  bay. 

Figures  8  and  9  indicate  high  salinity  water  was 
associated  with  Bob  Sikes  Pass.   To  estimate  if  there  have 
been  changes  in  the  bay  salinity  after  construction  of  the 
pass,  area  weighted  estimates  of  total  salt  in  the  bay  were 
calculated  and  plotted  against  average  river  flow.   Figure 
58a  shows  data  from  Dawson  (1955a)  plotted  in  this  way  for 
a  period  just  prior  to  the  construction  of  the  pass.   Figure 
58b  shows  data  collected  in  1973  plotted  in  the  same  way. 
Both  plots  indicate  that  total  salt  in  the  bay  was  inversely 
correlated  with  river  flow.   Comparison  of  the  two  plots 
suggests  that  during  periods  of  low  river  flow  there  was 
little  difference  before  and  after  the  shrimp  boat  pass 
was  cut.   However,  the  plots  suggest  that  during  high  river 
flow  periods  there  was  more  salt  in  the  bay  after  construc- 
tion of  the  pass.   If  more  salt  was  in  the  bay  this  sug- 
gests that  the  degree  of  mixing  of  bay  and  river  water  has 
decreased  and  that  river  water  is  reaching  the  Gulf  of 
Mexico  faster  than  before  the  pass  was  cut.   It  seems 
possible  that  the  shrimp  pass  has  caused  this  general  in- 
crease in  salinity.   More  data  should  be  collected  to  sub- 
stantiate fully  the  suggestion  made  here. 


Fig.  58.   Plot  of  total  salt  versus  river  flow. (a)  Data 
from  Dawson  (1955a)  taken  just  prior  to  con- 
struction of  Bob  Sikes  Pass,   (b)  Data  from 
Estabrook  (1973)  and  this  study.   Average  river 
flow  was  calculated  using  average  daily  river 
flow  for  a  two-week  period  prior  to  salinity 
measurements . 
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Salinity  control  has  implications  for  the  oyster 
industry.   It  appears  that  oyster  communities  do  best  in 
areas  that  have  seasonally  changing  salinities  (Menzel 
et  al.  ,  1957).   One  reason  for  this  is  that  many  oyster 
predators  cannot  withstand  low  salinity  water.   Each  year 
when  salinity  drops  these  species  are  either  killed,  migrate 
to  other  areas  or  become  dormant  and  oyster  predation  de- 
creases.  Additionally,  the  influx  of  freshwater  kills  many 
fouling  organisms  that  are  attached  to  oyster  shell.   After 
these  organisms  die  cultch  is  clean  and  available  for  new 
oyster  spat  (Gunter,  1953).   Other  mechanisms  are  also  prob- 
ably active  that  favor  oyster  dominance  in  estuarine  areas 
with  changing  salinities. 

Data  from  Menzel  and  Cake  (1969)(Fig.  59)  show  a 
fairly  clear  relationship  between  environmental  change  (as 
indicated  by  salinity  range)  and  diversity.   Low  diversity 
was  associated  with  large  yearly  salinity  changes.   High 
diversity  was  associated  with  areas  with  stable  salinities, 
in  this  case  areas  near  the  ocean.   Oyster  bars  in  Apala- 
chicola  Bay  were  associated  with  areas  of  moderate  salinity 
change  (15-20  ppt).   Two  stations  were  deleted  from  this 
plot.   Both  had  diversities  lower  than  expected.   They 
were  located  at  the  river  mouth  where  many  factors  such 
as  sedimentation  and  dredging  may  be  controlling  diversity. 

Shown  in  Fig.  60  is  a  long-term  plot  of  the  lowest 
yearly  salinity  recorded  at  two  oyster  bars  in  Apalachicola 


Fig.  59.   Plot  of  maximum  salinity  range  per  year  at  13 
stations  in  Apalachicola  Bay  versus  yearly  av- 
erage species  diversity  recorded  at  that  sta- 
tion.  Data  were  from  Menzel  and  Cake  (1969). 
Diversity  was  reported  as   number  of  species 


number  of  individ- 
ual s 


Numbers  on  the  diagram  correspond  to  station 
numbers  reported  by  Menzel  and  Cake  (1969). 
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Bay.   Station  II  is  near  the  John  Gourie  Bridge.   Station 
III  is  near  the  eastern  end  of  St.  Vincent  Island.   Data 
and  station  numbers  are  from  Menzel  ejt  af[.  (1966).   There 
is  some  indication  of  increasing  salinity  at  station  III 
especially  after  construction  of  Jim  Woodruff  Dam  and  the 
boat  pass.   Station  II  does  not  show  this  trend.   The  lower 
graph  in  this  figure  shows  monthly  salinity  values  at  these 
same  oyster  bars.   Station  II  was  a  \/ery   productive  reef, 
had  few  predators,  and  was  subjected  to  a  well-defined  sea- 
sonal dip  in  salinity.   Station  III  was  depleted  with  only 
small  spat  present.   Predators  were  common.   Salinities  at 
this  reef  did  not  drop  as  low  as  at  Station  II  and  were 
often  above  30  ppt.   Menzel  ejt  aj^.  (1966)  noted  that  their 
study  was  done  during  a  dry  period  and  this  may  be  in  part 
the  reason  for  the  lack  of  adult  oysters  on  Reef  II.   This 
reed  had  previously  been  quite  productive.   Several  local 
oystermen  said  that  oysters  were  scarce  in  Reef  III  during 
the  1974-75  season. 

-^      The  data  suggest  that  salinities  are  higher  and 
have  less  seasonal  change  now  than  before  the  pass  was  con- 
structed.  Since  stable  salinity  regimes  do  not  favor  oyster 
production  (as  suggested  in  Fig.  59)  a  quantitative  assess- 
ment of  the  boat  pass  and  its  relation  to  salinity  and 
stratification  may  be  advisable,  especially  if  some  addi- 
tional loss  of  freshwater  associated  with  a  proposed  navi- 
gation dam  on  the  Apalachicola  River  can  be  anticipated. 
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If  the  pass  proves  to  be  a  major  factor  in  changing  salini 
ty  and  mixing,  it  may  be  possible  to  modify  the  pass 
slightly  so  that  freshwater  is  retained  and  mixed  in  the 
bay  while  the  recreational  and  commercial  use  of  the  pass 
can  be  continued. 


Urban  Housing  Developments 

The  proposed  island  housing  development  has  stimu- 
lated local  interest  in  the  consequences  of  attempting  to 
have  a  fishing  and  recreation  area  coexist  in  close  prox- 
imity to  one  another.   Sewage  disposal  schemes  have  been 
proposed  by  the  developer  as  well  as  a  monitoring  program 
to  detect  changes  in  water  quality  associated  with  the  de- 
velopment.  The  impact  of  the  new  development  in  terms  of 
energy  flow  at  the  county  level  was  given  in  Table  11.   On 
a  yearly  basis  the  new  development  represented  a   65%  in- 
crease in  purchased  energy  flows  in  the  county.   Compared 
to  energy  flows  in  other  parts  of  the  state  and  nation, 
Franklin  County  would  still  be  less  developed  with  an  in- 
vestment ratio  of  0.37.   This  ratio  may  be  low  compared  to 
the  national  average  reported  in  Odum  and  Brown  (1975)  due 
to  calculation  differences.   The  national  ratio  had  natural 
inputs  calculated  by  dividing  annual  sunlight  by  a  factor 
of  2,000.   Here  all  natural  energies  were  summed  including 
sunlight  and  thus  the  divisor  in  the  investment  ratio  was 
larger.   Using  Odum's  method  the  investment  ratio  becomes 
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0.45  which  still  indicates  low  levels  of  development  rela- 
tive to  the  national  average. 

In  Fig.  53b  a  calculation  was  made  comparing  urban 
and  natural  energy  inputs  assuming  that  the  proposed  de- 
velopment was  expanded  to  include  the  entire  island.   In 
this  case  the  investment  ratio  became  3.1  suggesting  that 
this  level  of  development  would  overdevelop  the  county  com- 
pared to  the  rest  of  Florida  and  the  nation. 

In  Fig.  53c  a  third  evaluation  of  the  development 
was  done  using  an  area  around  the  development  rather  than 
the  whole  county.   If  the  energy  flows  associated  with  the 
proposed  development  (Fig.  55)  are  compared  to  natural 
energy  flows  in  the  full  county,  the  investment  ratio  in- 
creases to  0.37.   The  implication  from  this  may  be  that  a 
development  of  the  proposed  size  would  not  overdevelop  the 
county.   This  calculation  was  done  to  see  if  the  density  of 
the  development  was  high  relative  to  surrounding  areas. 
Could  other  similar  developments  be  added  to  this  one  and 
not  overdevelop  the  local  area?  As  shown  in  Fig.  53c  the 
investment  ratio  was  quite  high  (9.63),  well  above  the 
national  average.   If  construction  costs  were  not  included, 
the  ratio  was  still  high  (>  5).   The  ratio  suggests  that  at 
the  densities  proposed,  high-qual ity  purchased  energies  will 
not  find  sufficient  natural  energies  for  matching.   There 
may  be  an  overdependence  on  purchased  energy  flows. 
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Several  simulation  experiments  with  the  regional 
model  were  also  aimed  at  estimating  the  impact  of  large 
developments  based  on  external  energy  sources.   The  pro- 
posed island  development  is  of  this  type  because  its  suc- 
cess depends  largely  on  excess  energy  being  attracted  from 
outside  the  county. 

Figure  44  shows  the  impact  of  an  attempt  to  develop 
the  county  which  fails.   In  this  run  there  was  a  temporary 
surge  in  building  (as  represented  by  the  structure,  capital, 
and  image  storages).   However,  it  was  hypothesized  in  this 
experiment  that  external  energies  (or  capital  with  which 
to  direct  energies  to  Franklin  County)  were  not  available 
for  completion  of  the  project.   The  simulation  shows  a 
"boom  and  bust"  pattern  with  no  lasting  effects.   In  this 
simulation  the  oyster  fishery  did  not  close  because  of 
water  quality  changes.   Based  on  current  information  the 
proposed  island  development  would  have  a  yearly  investment 
of  about  $9.9  x  10  /yr  (about  20  times  the  current  invest- 
ment rate)  and  thus  would  be  equivalent  to  something  be- 
tween cases  2  and  3  in  Fig.  44.   Since  all  prices  in  this 
simulation  were  constant,  the  effect  may  be  more  like  case 
2  than  case  3.   In  fact,  if  energy  limitations  in  the  coun- 
ty as  a  whole  are  as  severe  as  some  suggest,  there  may  be 
no  change  at  all  because  the  development  will  never  start. 

Several  economists  doing  work  in  the  Franklin 
County  area  (Colberg  and  Windham,  1965;  Rockwood,  1973) 
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have  stressed  the  fact  that  the  region  has  a  limited  eco- 
nomic base  and  is  susceptible  to  economic  problems  should 
anything  interfere  with  the  fishery.   According  to  pro- 
ponents of  the  island  development  this  project  would  pro- 
vide additional  employment  and  relieve  the  county  dependence 
on  fisheries.   Several  simulation  experiments  were  done  to 
investigate  this  problem. 

In  the  simulation  shown  in  Fig.  43  tourist  money 
to  the  county  was  allowed  to  expand.   As  the  simulation 
shows  the  effect  of  the  seasonal  nature  of  the  oyster  fish- 
ery lessened  as  the  tourist  input  increases  (case  4).   In 
this  simulation  there  was  closure  of  the  fishery.   The 
impact  of  fishery  closure  was  not  large  because  when  closure 
occurred  the  oyster  industry  was  a  small  part  of  county 
income.   Growth  continued.   However,  in  this  simulation, 
prices  were  constant  and  the  income  associated  with  tourism 
was  about  2  times  as  large  as  investment  associated  with 
the  proposed  development.   In  Fig.  48  the  inflation  rate 
for  both  purchases  and  sales  was  at  20%  per  year.   In  case 
(b)  the  investment  rate  was  similar  to  the  one  for  the  pro- 
posed development.   In  this  case  the  dependence  on  the  in- 
dustry was  still  evident  as  shown  by  the  response  of  the 
capital  storage.   The  inflation  rate  of  20%  per  year 
leveled  growth  and  the  decline  was  severe  because  the  oys- 
ter industry  was  closed  (indicated  by  the  spikes  in  the 
capital  storage)  because  of  high  coliform  counts.   In  this 
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case  investment  did  not  help  the  county  economy  over  the 
long  run  because  of  external  energy  shortages  and  because 
the  oyster  industry  was  closed. 

A  third  possible  outcome  was  explored  in  the  simu- 
lation shown  in  Fig.  47.   In  this  case  there  was  an  infla- 
tion rate  of  10%  per  year  and  investment  increases  to  10 
and  20  times  the  normal  rate.   In  cases  (a)  and  (b)  there 
was  some  growth  and  then  a  gradual  decline  because  of  con- 
tinued inflation.   Investments  and  prices  for  oysters  and 
other  sales  would  have  to  increase  to  offset  the  effects  of 
inflation  if  growth  were  to  continue.   Case  (c)  shows  model 
response  under  investment  conditions  similar  to  those  of 
the  proposed  development  and  inflation  rates  for  sales  and 
purchases  of  10%  per  year.   It  was  hypothesized,  however, 
that  due  to  the  development  the  oyster  industry  was  lost 
for  one  reason  or  another.   The  simulation  suggests  that  a 
development  of  the  size  proposed  is  not  large  enough  to 
cover  losses  should  it  interfere  with  the  oyster  fishery. 
As  shown  in  Fig.  47  county  structure  goes  into  a  decline 
upon  the  loss  of  the  oyster  fishing  even  though  investment 
in  the  development  continued. 

The  decline  in  county  structure  would  have  been  more 
severe  if  development  stopped  because  of  poor  water  quality 
conditions  in  the  bay.   It  seems  wise  in  this  case  to  devise 
mechanisms  for  protecting  both  the  oyster  fishery  and  the 
development  so  as  to  assure  coexistence  of  both. 
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Factors  of  Energy  Investment  Return 

The  energy  return  from  energy  investment  is  dis- 
cussed in  this  section  including  the  oyster  industry  and 
the  conversion  of  river  head  to  electricity. 

Oyster  Industry 

Calculations  showing  the  energy  basis  of  the  oyster 
fishery  were  shown  in  Fig.  52,   The  energy  basis,  expressed 
in  fossil  fuel  equivalents,  was  mainly  natural.   Harvesting 
was  most  of  the  purchased  energy  input  (about  12%).   Energy 
in  sewage  from  the  town  and  planting  of  oyster  cultch  was 
relatively  small.   On  an  energetic  basis,  the  present  oys- 
ter industry  is  not  highly  dependent  on  purchased  fuels 
within  the  county.   Hence,  the  oyster  industry  may  be  less 
susceptible  to  rising  fuel  prices. 

The  matching  of  natural  and  purchased  energies  in 
the  oyster  industry  was  traced  from  input  of  natural  ener- 
gies to  final  consumption.   At  dockside  (P)  the  investment 
ratio  was  0.18  (2.2/12).   When  oysters  left  the  county  the 
ratio  was  1 . 35  ( 1 6. 3/1 2 . 0)  and  when  consumed  the  ratio  was 
2.15  (25.8/12).   The  calculation  suggests  that  the  natural 
energies  acting  in  support  of  the  oyster  industry  eventually 
attract  fossil  fuels  in  the  same  ratio  as  the  average  of 
other  processes  in  the  country.   This  may  indicate  that 
more  f uel -i ntensi ve  oyster  fishing  will  not  be  competitive. 
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Additionally,  if  the  ratio  of  purchased  to  natural  energies 
continues  to  drop,  the  oyster  processing  chain  may  shorten 
with  more  use  of  oysters  in  the  county. 

Most  oyster  reefs  in  Apalachicola  Bay  occur  in  the 
areas  occupied  by  the  medium  salinity  plankton  ecosystem 
(Fig.  7).   This  area  covers  80,000  acres.   Rockwood  (1973) 
suggested  that  this  area  could  support  many  more  oyster 
reefs  than  at  present.   Energy  used  in  maintaining  the 
present  reefs  (Fig.  52a)  was  compared  to  total  sunlight  in- 
put to  the  medium  salinity  plankton  system.   Total  sunlight 
was  used  as  an  estimate  of  total  energy  available  for  oyster 
growth  (5  x  10^^  Kcal/hr  of  sunlight).   The  calculated  energy 
use  by  oysters  was  about  2.5  x  10^^  Kcal/yr,   Comparison 
of  these  numbers  suggests  that  the  bay  may  be  able  to  sup- 
port twice  the  present  population  of  oysters.   Oyster  bars 
are  supported  by  large  areas  of  the  estuary  through  water 
movement.   This  points  for  the  need  for  full  estuarine 
management  when  considering  the  management  of  oyster  reefs. 
If  must  be  noted  that  this  calculation  assumed  that  no  other 
energy  subsidies  were  involved  in  oyster  production. 
Higher  densities  could  be  achieved  in  the  bay  only  by  add- 
ing more  purchased  energies  such  as  supplemental  feeding. 
Intensive  aquaculture  is  an  example.   Per  area  the  output 
of  oyster  meat  from  Apalachicola  Bay  was  about  8  x  10^ 
Kcal/acre/yr.   The  harvest  came  from  about  8,000  acres  of 
oyster  beds  in  the  bay.   The  yield  for  pasture  beef  in 
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Florida  was  20  x  10   Kcal/acre/yr.   In  beef  production, 
about  30%  of  the  input  energy  was  from  purchased  fuels 
(De  Bellevue,  1975).   Meat  production  from  estuarine  areas 
compares  favorably  with  alternatives. 

Simulation  results  shown  in  Figs.  45  and  46  sug- 
gest special  energy  inputs  to  the  bay  would  be  required  if 
oyster  yields  were  to  increase  much  beyond  the  present 
yield.   In  these  simulations  an  increase  in  harvest  effort 
by  a  factor  of  2  caused  large  declines  in  the  standing 
stock  of  oysters.   If  additional  inputs  were  developed 
using  purchased  energies,  output  may  increase  but  so  also 
would  dependence  on  fossil  fuels.   The  cultch  planting  pro- 
gram may  be  an  example  of  a  mechanism  for  increasing  yields 
in  a  fashion  that  is  not  heavily  dependent  on  purchased 
fuel s . 

The  oyster  fishery  represents  a  large  net  yield- 
ing process  for  the  county.   On  a  fossil  fuel  basis,  total 

harvest  effort  was  1.64  x  10   Kcalpp^/yr  in  1970.   Oyster 

11  fi 

sales  represented  1.63  x  10   Kcalppc/yr  ($6.5  x  10  x 

25,000  Kcal/$).   This  represents  a  10  to  1  return  on  the 

investment  as  represented  by  harvest  effort. 

Energy  Costs  of  Basin  Formation 


In  the  calculation  given  in  Fig.  56  the  energy 
quality  of  river  head  came  cut  to  be  1.42  times  as  concen- 
trated as  fossil  fuel  equivalent  calories.   This  indicates 
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that  the  potential  energy  of  elevated  water  is  of  high 
quality.   Young  ejt  al-  (1974)  calculated  a  slightly  higher 
ratio  using  information  for  larger  dams.   High  quality  of 
elevated  water  has  been  utilized  for  many  years  by  the 
construction  of  hydroelectric  dams  on  many  major  river 
systems  around  the  world.   A  suitable  basin  is  a  prerequi- 
site for  most  hydroelectric  facilities.   In  the  creation  of 
the  reservoir  at  Jim  Woodruff  Dam  a  considerable  area  of 
flat  land  was  flooded.   Some  terrestrial  land  was  converted 
to  a  use  with  less  metabolism  as  a  reservoir.   This  contribu- 
tion was  75%  of  the  total  cost  and  may  be  similar  to  geo- 
logical costs  in  other  river  systems  spent  on  creating  more 
relief  and  storage  capacity. 


Energetic  Basis  of  Regional  Coastal  Image 

What  factors  attract  money  and  people  to  the  Frank- 
lin County  region?   In  interviews,  some  suggested  it  was 
the  natural  resources  of  the  area.   Others  felt  it  was  re- 
lated to  the  existing  structure  in  the  area.   Several  com- 
binations of  rates  and  storages  were  suggested.   It  was 
hypothesized  that  there  is  in  every  region  a  factor  which 
represents  the  ability  of  an  area  to  attract  outside  in- 
vestments and  thus  amplify  that  region's  ability  to  do 
work.   We  chose  to  call  this  factor  image  and  to  make  it 
proportional  to  energy  flows.   Image  had  inputs  from  the 
major  work  processes  occurring  in  the  county.   These  included 
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inputs  from  estuarine  and  terrestrial  metabolism  and  urban 
activities  as  represented  by  the  major  input  to  structure 
(Fig.  28a).   Each  input  to  image  was  adjusted  so  that  the 
energy  flow  along  the  pathway  was  in  fossil  fuel  equivalent 
ki localories.   Energy  flows  of  different  quality  were 
weighted  equally  in  their  ability  to  do  work. 

If  the  attractiveness  of  a  region  is  based  on  the 
energy  flows  avai 1 abl e,  the  coastal  zone  would  seem  to  have 
a  naturally  high  image  because  there  are  special  energy 
flows  associated  with  coastal  areas  not  generally  found 
elsewhere.   These  include  waves  and  tidal  energy,  potential 
and  kinetic  energy  of  river  flow,  coastal  currents,  mixing 
energy  available  from  the  concentration  difference  between 
river  water  and  ocean  water,  steady  winds  and  perhaps  others 
not  recognized  here.   Energy  inputs  to  the  estuarine  area 
in  Franklin  County  were  more  than  twice  what  they  were  in 
other  areas  of  the  county.   If  the  theory  is  valid,  these 
special  energy  flows  will  attract  more  external  investments. 
Over  50%  of  the  population  of  the  U.S.  lives  in  counties 
bordering  on  the  ocean  or  Great  Lakes  (Ketchum,  1972). 
Most  urban  centers  in  Florida  occur  in  coastal  areas. 

In  this  dissertation,  the  image  of  the  region  was 
based  on  both  urban  and  natural  energy  flows.   Odum  and 
Brown  (1975)  suggested  that  the  ultimate  ability  of  a  re- 
gion to  attract  high-quality  energies  depends  on  the  amount 
of  lower-quality  natural  resident  energies  available  for 
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interaction  with  high-quality  flows.   High-quality  work 
can  be  done  by  high-quality  flows  only  when  there  are  low- 
quality  energies  with  which  to  interact.   If  these  are  not 
available,  the  high-quality  flows  can  not  be  fully  used  in 
a  high-quality  fashion  and  the  system  produces  less  to 
exchange  for  more  purchased  fuels. 


World  Conditions  of  Energy  Prices  and 
Future  Trends  in  Coastal  Franklin  County 


The  conditions  of  world  fuel  prices  as  they  may 
affect  future  trends  in  Franklin  County  were  explored 
using  the  regional  simulation  model  shown  in  Fig.  28a. 
Sharp  growth  trends  were  observed  only  when  external  ener- 
gies were  not  limiting.   These  familiar  trends  were  shown 
in  Figs.  42  and  43. 

The  most  dominant  trend  in  the  regional  model  was 
that  growth  in  the  urban  sector  leveled  and  declined  slowly 
when  prices  were  allowed  to  inflate  at  rates  similar  to 
those  in  the  nation.   This  pattern  was  consistent  under 
a  variety  of  conditions.   In  Fig.  41  prices  for  oyster 
sales  and  prices  for  goods,  fuels  and  services  were  both 
inflating  at  10%  per  year.   It  was  hypothesized  that  some 
modification  could  be  made  in  the  bay  with  little  cost 
(perhaps  similar  to  the  cultch  planting  program)  that  would 
increase  the  recruitment  of  adult  oysters.   As  shown,  even 
increasing  the  recruitment  rate  by  a  factor  of  five  was 
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not  sufficient  to  offset  declines  in  the  urban  sector.   In- 
creasing the  harvest  effort  (Figs.  45  and  46)  did  not  off- 
set declines  because  oyster  stocks  became  limiting.   De- 
clines were  more  severe  when  the  inflation  rate  of  oysters 
was  less  than  for  goods,  fuels  and  services  (Fig.  46b). 
There  were  only  small  declines  when  the   investment  rate 
was  increased  to  ten  times  the  current  rate  and  some  growth 
when  the  rate  was  increased  to  20  times  the  current  rate 
(Fig.  47b).   Declines  and  "boom  and  bust"  trends  were  sharper 
with  large  investment  rates  and  20%  per  year  inflation  (Fig. 
48).   In  all  cases,  declining  energy  limited  long-term 
growth.   In  Fig.  43a  (case  4)  the  oyster  industry  was  lost 
due  to  high  coliform  counts.   With  rich  energies  available 
oyster  income  was  only  a  small  part  of  county  income  and 
loss  of  the  industry  had  little  effect  on  the  growth  pattern 
When  external  energies  were  limiting  loss  of  the  oyster  in- 
dustry caused  a  reversal  in  growth  trends  (Fig.  47c)  and 
a  decline  in  structure. 

Growth  based  on  the  oyster  industry  had  a  stable 
pattern  with  leveling  as  shown  in  Fig.  49,  a  characteris- 
tic of  all  runs  based  on  internal  resources.   Leveling  was 
initiated  when  oyster  stocks  declined. 

When  prices  were  allowed  to  inflate  the  yearly 
change  in  capital  became  more  pronounced.   This  resulted 
because  both  inputs  and  outputs  from  the  capital  storage 
increased  while  the  stock  remained  about  the  same.   Hence 
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the  stock  was  sensitive  to  interruptions  in  input  such  as 
that  caused  by  the  seasonal  oyster  harvest.   The  pulse  be- 
came larger  when  harvest  effort  was  increased  or  when  in- 
creased investment  created  new  structure  which  required  ex- 
penditures for  maintenance. 

This  model  suggests  that  when  external  energies 
are  limiting,  efforts  to  continue  growth  lead  to  increased 
turnover  of  capital  and  sensitivity  to  interruptions.   At- 
tempts at  growth  based  on  external  investment  were  shown 
in  Fig.  48.   In  case  (a)  declines  were  smooth  and  gradual. 
In  cases  (b)  and  (c)  there  was  a  boom  and  bust  cycle  with 
large  yearly  changes  in  capital  and  in  case  (c)  loss  of  the 
oyster  industry.   The  end  result  in  all  three  cases  was 
about  the  same.   The  model  suggests  that  when  external 
energies  become  limiting,  as  they  may  be  now,  declines 
in  structure  may  result.   Large  investments  aimed  at  growth 
were  not  long- 1 asting ,  created  large  yearly  changes  in 
capital,  and  in  one  case  caused  the  loss  of  the  previously 
successful  oyster  industry. 

While  this  model  predicts  leveling  or  declines 
in  the  county  due  to  external  energy  shortages,  there  still 
may  be  some  significant  growth  in  the  area  in  the  future. 
This  growth,  if  it  occurs,  may  be  based  on  the  redistribu- 
tion of  purchased  energies  from  presently  developed  urban 
areas  in  other  parts  of  the  country.   During  initial  periods 
of  growth,  based  on  large  energy  sources,  there  was  a 


302 


premium  on  getting  assets  generated  with  which  to  feed  back 
to  the  energy  source.   Those  areas  with  the  most  early 
assets  won  when  in  competition  with  other  systems  with 
smaller  assets.   The  result  was  spotty  development  with  un- 
equal distribution  of  fossil  fuels.   However,  when  energy 
sources  become  limiting,  as  they  may  be  at  present,  this 
advantage  declines  both  because  further  efforts  to  process 
more  fuels  are  not  successful  and  because  maintenance  costs 
are  high.   The  work  done  in  these  areas  per  unit  of  fuel 
expended  may  decline.   Rural  areas  with  more  natural  ener- 
gies available  for  use  with  fossil  fuels  may  become  more 
competitive  and  attract  high-quality  energies  away  from 
present  urban  areas.   High-quality  energies  may  become  more 
evenly  spread  across  the  landscape.   It  seems  possible  that 
this  may  happen  in  Franklin  County  because  of  the  abundance 
of  natural  energies  available. 


Fisheries  and  Urban  Economy 


The  future  coexistence  of  coastal  fisheries  and 
urban  development  was  a  central  issue  in  this  paper.   Sim- 
ulation models  and  calculation  results  suggest  several 
trends  and  strategies  that  may  have  use  in  Franklin  County 
or  other  similar  areas. 

■■^.    Energy  calculations  suggested  that  the  support  of 
each  acre  of  oyster  reef  was  derived  from  about  five  acres 
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of  estuary.   The  present  ratio  of  oyster  reef  area  to  es- 
tuary in  Apalachicola  Bay  is  about  one  to  ten.   It  may  be 
possible  to  expand  the  natural  oyster  stock  by  a  factor  of 
two  without  the  addition  of  large  purchased  energies  sub- 
sidies. 

The  evaluation  of  the  oyster  industry  suggested  that 
about  37%  of  the  value  of  oysters  was  exported  from  the 
county  and  matched  with  fossil  fuels  elsewhere.   If  the 
oysters  were  more  fully  processed  in  the  county  more  high- 
quality  energy  may  be  attracted  to  the  county  and  add  to 
the  energy  basis  of  the  area.   The  low  ratio  of  purchased 
to  natural  energies  suggests  that  there  would  be  a  good  re- 
turn on  additional  inputs  of  purchased  energies. 

Simulation  results  from  the  regional  model  (Fig. 
28a  and  Fig.  47)  and  investment  ratio  calculations  suggest 
that  addition  of  more  purchased  energies  would  bring  the 
county  economy  into  a  more  competitive  position  relative  to 
the  rest  of  the  nation.   In  this  sense  development  would  be 
good.   However,  simulation  results  also  suggest  (Fig.  47c) 
that  new  developments  that  affect  water  quality,  and  hence 
the  oyster  industry,  should  be  kept  well  insulated  from  the 
estuary.   If  this  can  be  done,  coexistence  may  be  quite 
possible.   The  disposal  or  recycle  of  some  wastes  from  ur- 
ban developments  might  be  done  in  the  abundant  wetlands  in 
the  county.   This  type  of  coupling  would  make  purchased  fuels 
available  for  other  uses,  and  may  provide  a  needed  buffer 
between  estuarine  fisheries  and  urban  activities. 


APPENDICES 


APPENDIX   A 
SYMBOLS    USED    IN    MODEL    DEVELOPMENT 


The    symbols    used    in    diagrams    of   models    are    those    of 
the   enery   circuit   language    developed    by   H.    T.    Odum.      Each 
symbol    has    both    a    verbal    meaning    and    an   exact  mathematical 
equivalent   which    can    be    found    in    Odum    (1971a,    1972a). 

Forcing  Function.     An  external 
source  of  energy  with  or  without 
materials  whose  driving  forces  are 
independent  of  model   behavior. 
Program  can  be  constant,  sinosoidal, 
etc.  and  is  controlled  from  outside 
the  model . 

Flow  Limited  Forcing  Function.     An 
external   source  of  energy  with  or 
without  materials  whose  input  can 
be  a  limiting  factor  due  to  inter- 
actions within  the  model. 
kJ. 


^  =  ^0-  kJ.X;  J  =  k(^-^)(X) 
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Forcing  Function.  An  external 
source  of  energy  with  or  without 
materials  whose  input  is  determined 
by  some  variable  within  the  model 
(X).  Inflow  can  only  be  limited 
by  the  variable  with  which  the 
forcing  function  interacts. 

Storage  Module.  Represents  a  stor- 
age of  energy  of  materials  within  a 
system  where  a  quantity  is  stored  as 
the  balance  of  inflov;s  and  outflows 
(^  =  J  -  kQ)  and  where  outflow  in- 
cludes depreciation. 

Self-Maintaining  Consumer  Model.  A 
group  module  which  represents  a  con- 
sumer unit  including  a  combination  of 
a  storage  module  and  at  least  one 
multiplier  where  energy  stored  in  one 
or  more  places  in  the  module  is  fed 
back  to  do  work  on  processing  input 
energy  to  that  unit;  response  is 
autocatalytic  if  the  above  features 
are  included.  The  group  symbol  is  of- 
ten used  to  organize  model  components. 
When  used  in  this  way,  it  does  not  im- 
ply additional  pathways  beyond  those 
actually  shown. 
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Production  and  Regeneration  Module. 
A  group  module  representing  an  inter- 
active production  process  and  stor- 
age. Normally  used  to  depict  green 
plant  photosynthesis.  On  a  region- 
al scale  the  module  represents  the 
production  and  consumption  of  entire 
ecosystems  (P/R).  Details  of  rela- 
tionships in  a  particular  model  are 
shown  within  the  group  symbol. 


>f 


Pathway.  Shows  a  flow  of  energy 
with  or  without  materials  which  is 
proportional  to  a  quantity  in  stor- 
age or  external  sources  at  each  end 
(J  =  k(Q^  -Q^)).  The  heat  sink  re- 
presents energy  losses  due  to  fric- 
tional  forces  and  backforce  along 
the  pathway. 


Adding  Junction.  Shows  the  inter- 
section of  two  pathways  capable  of 
adding.  Arrow  indicates  direction 
of  flow  and  absence  of  any  backforce. 


Money  Pathvjay.  Dashed  line  indicates 
a  flow  of  money  with  arrow  indicating 
direction. 
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Force  from  a  Flow  Symbol.  Flow 
rate  of  one  pathway  (J  )  delivers 
a  force  X  that  is  proportional  to 
the  sensed  flow  and  derives  its 
energy  from  it. 

Price  Transactor.  Symbol  indicates 
an  economic  transaction  with  price  (P) 
the  ratio  of  money  flow  to  energy 
flow  (J2/J1).  Price  may  be  constant 
or  may  vary  in  a  variety  of  ways. 
Heat  sink  indicates  the  energy  cost 
of  maintaining  transactions. 

Multiplicative  Workgate.  Symbol  in- 
dicates intersection  of  two  pathways 
coupled  to  produce  an  outflow  propor- 
tional to  the  product  of  the  forces 
driving  both  flows.  General  re- 
sponse is  a  limiting  factor  type 
(J  =  kN^N2). 

Drag  Action  Workgate.  Symbol  indi- 
cates an  intersection  where  an  in- 
crease in  one  flow  has  a  retarding 
effect  on  the  output  flow  (J  = 
kN^(l  -  kNp)).  Sensor  symbol  indi- 
cates there  is  no  appreciable  loss 
from  Np  in  this  interaction. 
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Reversing  Multiplier.     Symbol   indi- 
cates an  intersection  where  flows 
may  reverse  direction  when  the  down- 
stream force  is  greater  than  the  up- 
stream force  [J  =  kN2(N^   -  N3)].     Np 
is  the  common  driving  force  variable. 


V 


Group  Logic.     Group  symbol   for  on-off 
processes  controlling  flow  J.     De- 
tails are  drawn  within  such  as  the 
Logic  Comparator  shown  next.     Heat 
sink  indicates  the  cost  of  maintain- 
ing logic  control . 

Logic  Comparator.  A  logic  module 
which  is  on  or  off  depending  on  a 
critical   threshold  (T). 


APPENDIX  B 
ECOSYSTEM  CLASSIFICATION  SYSTEM 

A  classification  of  ecosystem  types  was  developed 
for  Franklin  County  and  used  in  the  land-use  maps  in  Figs. 
6  and  7.   Categories  were  named  so  as  to  recognize  dominant 
energy  flows  or  vegetation  types. 

Urban  Systems 

Scraped  Land.   Lands  cleared  and  prepared  for  de- 
velopment. 

Light  Residential  .   Residential  areas  with  densi- 
ties of  four  units  per  acre  or  less. 

Medium-Dense  Residential.   Residential  areas  with 
densities  of  four  units  per  acre  or  more. 

Commercial  and  Industrial.   Commercial  and  indus- 
trial areas  within  the  city. 

Airport.   System  designation  used  here  to  identify 
position  of  local  airports.   Airports  in  Franklin  County 
not  used  for  regular  commercial  traffic. 
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Coastal  Systems 


Oligohaline  System.   System  existing  at  the  river 
mouth  area  of  estuaries  (normal  salinity  0.4-10  ppt). 
Characterized  by  large  seasonal  salinity  changes  due  to 
river  pulses,  turbid  water,  and  low  diversity.   Biomass  is 
often  high  being  composed  of  mollusks,  and  migrating  inver- 
tebrate and  fish  subsystems  (Copeland  ejt  aj[.  ,  1974). 

Grass  Flats.   Estuarine  system  occurring  in  shallow 
waters  with  moderate  currents,  fairly  stable  salinities  and 
low  turbidity.   In  turbid  bays,  grass  flats  occur  on  shallow 
sand  bars.   Characteristic  plant  species  include  Diplanthera 
and  Ruppia  in  intertidal  areas,  and  Thalassia ,  Syringodium 
and  H a  1 0 p h i 1  a  below  water  levels.   Zostera  is  dominant  in 
more  temperate  areas.   Animal  species  are  varied  and  often 
depend  on  the  grass  flat  for  part  or  all  of  their  life 
cycle  (Phillips,  1974). 

Salt  Marsh.   Estuarine  system  most  often  character- 
ized by  dense  stands  of  Sparti  na  and/or  Juncus  with  large 
seasonal  organic  matter  exports  to  adjacent  estuarine 
areas.   Marshes  develop  where  there  are  areas  of  inter- 
tidal soft  sediments  not  too  strongly  stressed  by  ocean 
waves  or  winter  ice.   Large  numbers  of  estuarine  species 
utilize  marsh  creeks  for  feeding,  protection  and  as  nur- 
sery areas  (Cooper,  1974). 

Oyster  Reef.   Estuarine  subsystem  found  where  cur- 
rents and  tidal  action  bring  suspended  material  that  serves 
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as  food  and  eliminates  waste  from  the  reef  community.   Bio- 
mass  of  oysters  and  associated  organisms  is  high,   diver- 
sity moderate  and  metabolism  high  and  dependent  on  sup- 
plies of  imported  organics.   Most  commercial  oysters  bars 
occur  in  the  medium  salinity  plankton  system  zone  or  inter- 
tidally  in  higher  salinity  areas  (Chestnut,  1974). 

Sand  Bar.   Intertidal  or  subtidal  sand  areas  formed 
by  current  action.   Characterized  by  shifting  sands,  lack 
of  rooted  vegetation  and  sparse  benthic  populations. 

Medium  Salinity  Plankton  System.   System  character- 
ized by  salinities  between  5-18  ppt  and  often  located  in 
the  middle  broad  portion  of  east  coast  estuaries.   Seasonal 
cycles  of  nutrients,  temperature,  phytopl ankton  and  zoo- 
plankton  production  are  evident  in  northern  and  temperate 
areas.   The  food  web  is  based  on  phytopl ankton  production 
and  supports  commercial  quantities  of  fish,  shrimp,  oysters 
and  clams.   These  estuaries  often  have  partial  stratifica- 
tion with  wedges  of  dense  salt  water  along  the  bottom. 
Anaerobic  conditions  do  not  generally  occur  on  the  bottom 
although  oxygen  is  often  lower  than  at  the  surface  (Bellis, 
1974). 

Coastal  Plankton  System.   The  coastal  system  is 
located  between  the  estuarine  zone  and  the  open  sea,  ex- 
tending to  a  depth  of  approximately  300  ft.   Phytopl ankton 
production  is  the  main  photosynthetic  food  source.   Res- 
piration and  nutrient  recycle  is  accomplished  both  by 
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zooplankton  and  diversified  benthic  communities.  The 
coastal  plankton  system  supports  commercial  fisheries 
(Odum  ejt  al.  ,  1974a). 

High  Velocity  Channel.   Estuarine  ecosystem  charac- 
terised by  high  velocity  currents  (2-20  miles/hr)  and  having 
reef-like  communities  of  encrusting  organisms  along  the 
sides  of  the  channel.   Bottoms  are  often  swept  clean.   In 
areas  where  waters  are  clear,  dense  stands  of  attached  macro- 
phytic  algae  are  present.   High  metabolic  activity  is  facil- 
itated by  the  strong  water  flow  delivering  organisms  and 
nutrients  and  flushing  wastes.   Where  high  velocity  ecosys- 
tems occur  near  open  ocean  areas  with  high  salinities,  di- 
versity is  moderately  high  (Odum  ejt  al.  ,  1974a). 

High  Energy  Beach.   Characteristic  ecosystem  occur- 
ring along  sandy  shores  receiving  strong  wave  action.   Well- 
sorted  sands  and  zones  with  distinct  chemical  and  physical 
characteristics  support  a  specialized  biota  of  sand  dwell- 
ing (psammon)  animals.   Different  portions  of  the  epipsammon 
live  on  the  sand  and  include  birds,  fish,  etc.;  the  endo- 
psammon  burrow  into  the  sand  and  include  snails,  bivalves 
and  Crustacea;  mesopsammon  live  between  sand  grains  and 
include  many  invertebrate  classes.   The  beach  is  primarily 
heterotrophic,  acting  as  a  large  filter  extracting  both 
particulate  and  dissolved  organics  and  nutrients  delivered 
by  ocean  waves  (Riedl  and  McMahan,  1974).   In  this  classi- 
fication the  high  energy  beach  extends  out  to  the  20  ft 
contour  or  500  yards  offshore,  whichever  comes  first. 
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Coastal  Dune.   Includes  area  extending  landward 
from  the  extreme  high  t-!de  line  to  the  beginning  of  pine- 
lands.   The  vegetation  is  adapted  to  salt  spray  and  occas- 
sional inundation.   Characteristic  species  include  sea 
purslane,  bermuda  grass,  seaside-pennywort,  sea  ort,  sea- 
side evening  primrose,  saw  palmetto  and  others.   Cabbage- 
palms  and  hardwoods  may  be  present  on  the  leeward  side  of 
old  established  dunes  (Kurz,  1942). 

Boat  Basins  and  Harbors.   These  new  estuarine  sub- 
systems are  often  characterized  by  depth  exceeding  natural 
flushing  capacity  resulting  in  deposition  of  sediments  and 
organic  matter.   Phytopl ankton  production  is  high  due  to 
high  nutrient  concentrations  but  the  euphotic  zone  is 
limited  and  respiration  generally  exceeds  photosynthesis. 
Oxygen  stress  on  the  bottom  can  be  severe  and  species  di- 
versity is  often  low. 

Dredged  Navigation  Channel.   These  new  estuarine 
subsystems  are  characterized  by  being  deeper  than  surround- 
ing areas,  steep-sided,  and  often  stratified  with  respect 
to  temperature,  salinity,  oxygen,  and  other  factors.   At 
times  oxygen  concentrations  may  be  very  low. 

Freshwater  Systems 


Shallow  Woodland  Pond.   Freshwater  ecosystems  char- 
acterized by  open  water  with  slow  exchange   and  seasonal 
stratification.   In  ponds,  the  littoral  zone  is  large 
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compared  to  the  limnetic  and  profundal  zones  and  is  the 
main  production  area.   Ponds  were  characteristic  of  the 
study  region  and  were  probably  warm  monomictic  with  cir- 
culation in  the. winter  (Odum,  1971).   Some  ponds  were  sur- 
rounded by  cypress,  and  many  have  rooted  aquatic  plants, 
dark  water  and  seasonal  water  level  fluctuation.   Included 
here  are  a  few  brackish  water  ponds  on  the  barrier  islands. 

Rivers.   Freshwater  ecosystem  characterized  by 
through-flowing  waters,  and  general  absence  of  stratifica- 
tion.  Attached  algae  and  diatoms  may  be  important  producers 
but  imported  organics  are  generally  the  primary  food  source. 
Rivers  characteristically  export  organics,  recycle  nutri- 
ents, and  support  seasonal  fish  and  insect  populations 
(Odum,  1971).   Two  types  of  rivers  appear  in  the  study 
area  and  include  low  pH,  blackwater  rivers  carrying  swamp 
drainage,  and  turbid  piedmont  rivers  carrying  upland  silts 
and  clays . 


Terrestrial  Systems 


Pine  Forest.   This  type  includes  all  planted  pine 
(slash)  stands  in  the  county  in  which  the  trees  are  large 
enough  to  be  recognized  from  aerial  photographs.   Those 
areas  recently  planted  in  which  trees  are  not  recognizable 
from  the  photographs  have  been  included  in  the  cleared/ 
planted  type.   Also  included  here  are  scattered  stands  of 
native  longleaf,  slash,  and  sand  pine  which  could  not  be 
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clearly  distinguished  from  planted  pines  on  aerial  photo- 
graphs. 

Cleared/Planted  Lands.   This  type  includes  lands  that 
have  been  cut  for  timber  and  recently  replanted  in  pine  or 
allowed  to  naturally  reestablish  vegetative  cover.   Ground 
cover  is  characterized  by  young  planted  pines  (not  recog- 
nizable from  aerial  photographs),  shrubs,  grasses,  and 
scattered  mature  pines  and  saw  palmetto.   No  attempt  was 
made  to  distinguish  between  wet  and  dry  habitats,   or  from 
natural  savannahs  as  recognized  by  Clewell  (1971). 

River  Swamp.   Climax  hardwood  system  found  in  sea- 
sonally flooded  areas  with  poor  drainage  (along  creeks, 
rivers,  sloughs,  and  depressions)  characterized  by  presence 
of  Fraxinum  carol ini ana  (Carolina  ash),  Taxiodium  distichum 
(bald  cypress),  Sabal  palmetto  (Cabbage  palmetto),  and 
Ulmus  f loridana  (white  or  American  elm).   Acer  rubrum  (red 
maple),  Nyssa  sylvatica  (Black  tupelo),  and  Liquidambar 
styracif 1 ua  (Sweetgum)  are  major  species  but  also  occur  in 
bay  swamps.   Quercus  nigra  (water  oak),  ^.  laurifol ia 
(laurel  oak),  Carya  aquatica  (water  hickory),  and  Magnol ia 
virginiana  (sweet  bay)  may  also  occur.   Cephalanthus  occi- 
dental is  (Cottonbush)  ,  Ilex  cassine  (cassena).  My  r  i  c  a  c  e  r  i - 
fera  (wax  myrtle)  and  C  a  r  p  i  n  u  s  carol ini ana  (American  horn- 
beam) are  the  major  understory  trees.   These  swamps  are 
more  fertile,  less  acidic,  and  have  more  flooding  than  bay 
swamps.   Tree  dominance  is  more  variable  than  in  bay  swamps. 
Litter  accumulation  is  slight.   Soils  are  alluvial. 
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Bay  Swamp.   A  broad  classification  which  includes 
titi,  bay,  cypress,  and  blackgum  dominated  swamps.   Soils 
are  very  poorly  drained  highly  organic  sands  of  the  Plummer- 
Rutledge  Association.   Clewell  (1971)  gives  a  full  descrip- 
tion of  each  swamp  type.  Also  included  here  are  savannahs 
characterized  by  frequent  fire,  poorly-drained  Plummer- 
Rutledge  soils,  and  dominant  herbaceous  vegetation.   Soil 
moisture  is  higher  than  in  the  pine-palmetto  flatwoods  and 
some  bays  (Clewell ,  1971  ). 

Freshwater  Marsh.   Characteristic  vegetation  includes 
pickerel  weed  (Pontederia  lanceolata),  fireflag  (Thai i a 
geniculota  L.),  water  lilies  (Nymphacea  spp. )  ,  U  t  r  i  c  u 1  a  r  i  a 
i  nf lata  (blatterwort) ,  cattail  (Typha  spp  .  )  .  spike  rush 
(Eleocharis  cellulosa),  bullrush  (Scirpus  spp.),  and  spat- 
terdock  (Nuphar  luteum).   Seasonally  to  continually  inun- 
dated by  water  with  slow  water  drainage.   Water  levels  gen- 
erally higher  than  in  wet  prairies  (Davis,  1943). 

Pine  Scrub.   Areas  characterized  by  sand  pine 
(  P i n u s  clausa  ) ,  scrub  oaks  (Quercus  myrti  fol  ia  ,  Cj.-  chap- 
ma  n i  i ,  and  £.  vi  rgi  niana  var .  mari  ti  na)  ,  palmetto  (Sere- 
noa  repens ) ,  and  various  lichens.   Wiregrass  is  absent. 
Soils  are  excessively  drained  and  sandy  soils  of  the  St. 
Lucie-Kureb-Rimi ni  association  (Kurz,  1942).   Included  here 
are  the  Franklin  County  barrier  islands. 

High  Pineland.   Areas  characterized  by  longleaf 
pine  (Pinus  palustris)  with  understory  of  bluejack  oak 
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(Quercus  incana).  sand-post  {Q^.    margaretta) ,  turkey  oak 
(^.  laevis)  and  sand  live-oak  (Q.  gemineta).   Wiregrass 
(Aristida  stricta)  is  abundant  with  other  herbaceous  spe- 
cies present.   Fire  is  frequent.   Soils  are  deep,  dry  sands 
of  the  Lakeland  group.   Also  called  sandhill,  turkey  oak 
ridge,  or  scrub  oak  (Clewell,  1971). 

Pine-Palmetto  Flatwoods.   Area  characterized  by 
longleaf  (P.  palustris)  and  slash  pine  (P.  el  1 iottii ) . 
Ground  cover  of  wiregrass,  saw  palmetto,  gallberry  and 
others.   Fires  are  frequent  in  the  natural  situation.   Soils 
are  sandy  with  an  underlying  hardpan  several  feet  below  the 
surface  and  are  predominantly  of  the  Leon  group  (Clewell, 
1971). 


APPENDIX  C 

SUPPLEMENTARY  INFORMATION  ON 
THE  DIURNAL  ACTIVITY  MODEL 


Table  C-1.   Data  used  in  simulation  of  the  diurnal  activity 
model  . 


Table  C-2.   Scaled  differential  equations  and  pot  settings 
used  in  simulating  the  diurnal  activity  model. 

Table  C-3.   Maximum  values  for  state  variables  and  forcing 
functions  used  in  diurnal  activity  model. 

Fig.  C-1.   Analog  patching  diagram  of  the  diurnal  activity 
model  shown  in  Fig.  19. 
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Table  C-2 


Scaled  differential  equations  and  pot  settings 
used  in  simulating  the  diurnal  activity  model. 


Phytoplankton  biomass 


0.05(10) 


0.3 


-  0.06 
^4 


k.  JOj 


0.04 
^7 


0.004 


Consumer  biomass 


=  0.011 
•^10 


100 


0.002 
^11 


15 


Dissolved   oxygen 
6. 


=  0.001 
k 


0.002 
k 


-  0.022 
k. 


0.35 


A        ^       -• 


+     0.027 
k 


13 
30 


14     ^    -" 


13      ^ 


-  0.010 
k„ 


T 
100 

30 

- 

'"31 
30 

\' 

30 

'"3" 
30 

""4" 
100 

0.021 
k 


23 
30 


+  0.168 
k2 


6       ^    -■ 


700 

0.30 

100 


Detritus 


100 


0.004 
k. 


0.003 
ko 


0.4 


+  0.004 
^7 


30 

100 

0.002 
•^19 


100 


-  0.006 

^7 


30 

100 
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Table   C-2   -   continued 


Ni  trogen 


0.3 


=  0.004 


+  0.109 


hsS 


N 


0.013 


+  0.122 
'^12*^4 


100 


0.001 
■^20 


0.3 


0.014 

hehi 


0.083(10) 


700 

0.3 

Source  limited 
sunlight  inflow 


700 


^  [s- 

0.330  .700. 

0.330  +  0.793 
P      0 

.0.30J 
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Table   C-3 

Maximum   values    for   state    variables    and    forc- 
ing   functions    used    in   diurnal    activity   model 
Scaled    potentiometer   settings    are   given    for 
each   state   variable   and   forcing   function. 
Potentiometer   settings   for  the   sine  wave   gen- 
erator  are   given   at   the   bottom. 


Maximum  Value  Scaled  Pot  Maximum  Value  Scaled  Pot 

of  State  Variable  Setting  of  Forcing  Function        Setting 

^1-.  n  o-,n  r    A    -,  0.001 


Cl^]  0.210  l^-] 


[^f]  0.547  [^]  0.25 


[3§]  0.547  [q^]  0.004 


[^]  0.197  [q-^]  0.004 


[^]  0.20  [%}  0.485 


sine  wave  settings 

I.e.  =  1 .00 
0)  =  0.262 
u)  =  0.262 


Fig.  C-1.   Analog  patching  diagram  of  the  diurnal  activity 
model  shown  in  Fig.  19.   Equation  coefficients 
are  identified  by  numbered  potentiometers. 
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APPENDIX  D 


SUPPLEMENTARY  INFORMATION  ON  THE  REGIONAL  SIMULATION  MODEL 
OF  FRANKLIN  COUNTY  AND  APALACHICOLA  BAY 


Table  D-1.   Data  used  in  the  regional  model  of  Franklin 
County  and  Apalachicola  Bay. 

Table  D-2.  Scaled  differential  equations  used  in  simula- 
tion of  the  regional  model  of  Franklin  County 
and  Apalachicola  Bay,  Florida. 

Table  D-3.   Scaled  state  variables  and  forcing  functions 

and  scaled  initial  condition  pot  settings  used 
in  simulating  the  regional  model  of  Franklin 
County  and  Apalachicola  Bay. 

Fig.  D-1.   Analog  patching  diagram  for  the  regional  model 
of  Franklin  County  and  Apalachicola  Bay, 
Florida. 


Fig.  D-2.   Analog  patching  diagram  used  to  generate  sun- 
light, river  flow  and  logic  control  of  oyster 
fishing  season  in  regional  simulation  model  of 
Franklin  County  and  Apalachicola  Bay,  Florida. 
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Table    D-3 


Scaled    state    variables    and    forcing    functions    and    scaled 
initial    condition    pot    settings    used    in    simulating   the 
regional    model    of   Frank.Tin    County   and   Apalachicola   Bay. 


Maximum  Value  of 
State  Variable 

Scaled 

Initial  Condition 

Pot  Setting 

Maximum  Value  of 
Forcing  Functions 

Scaled  Pot 
Setting 

Ql 

0.62 

^0 

See  Appendix, 
Table  D-2  (re- 
fer to  diagram 

15x10^ 

2.5x10^^ 

of  sine  wave) 

^2 

0.39 

^1 
2x10^ 

See  Appendix 

ISxlO'' 

Table  D-2 

^3 

0.20 

h 

a 

10x10^^ 

10x10^ 

^4 

0.077 

1.0 

10x10^ 

^5 

0.141 

^4 
20 

10x10^^ 

^6 

0.106 

^5 

7x10^ 

3x10^ 

Qy 

0.10 

^6 
1x10^ 

75x10^ 

^8 

0.50 

h 

1.04x10^ 

10x10^° 

% 

0.106 

^7- 

10x10^0 

10x10^ 

Forcing  functions  with  no  scaled  pot  setting  were  combined  with  trans- 
fer coefficient  pot  setting.     Combined  values  can  be  found  in  Appendix  Table  D-2. 


370 


Table    D-3   -    continued 


Maximum  Value  of 
State  Variable 

Scaled 

Initial   Condition 

Pot  Setting 

Max 
Fore 

imum  Value  of 
ing  Functions 

Scaled  Pot 
Setting 

^10 

0.110 

^8 

5x10^ 

1x10^ 

^11 

0.174 

^9 

5x10^ 

5x10^ 

Ql2 
2x10^ 

0.153 

^12 
1x10^ 

^13 
15x10^^ 

0.367 

^13 
1x10^ 

^20 
3x10^ 

^r 

0.50 

2.5x10^^ 

Pi 

0.1 

4x10"^ 

'z 

0.10 

4x10"^ 

^3 

0.10 

4x10"^ 

P4 

0.10 

0.473 

L 

0.548 

15x10^^ 

Fig.  D-1 .   Analog  patching  diagram  for  the  regional  model 

of  Franklin  County  and  Apalachicola  Bay,  Florida 
Equation  coefficients  are  identified  by  numbered 
potentiometers.   Patching  for  oscillating  forc- 
ing functions  and  logic  control  is  given  in  Fig. 
D-2. 
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APPENDIX  E 

FULL-SIZED  MAPS  OF  PRIMITIVE 
AND  PRESENT  LAND-USE  PATTERNS 


These  maps  are  on  file  in  the  library  of  the  Department 
of  Environmental  Engineering  Sciences,  University  of 
Florida.   Reduced  versions  of  the  maps  are  given  as  Figs 
6  and  7  in  the  main  text. 
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